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The  capability  of  applying  spatially  variable  inputs  to  field  crops  has  triggered 
new  opportunities  in  crop  management  (Site-Specific  Management,  SSM).  The  success 
of  SSM  depends  on  the  ability  to  relate  the  status  of  a crop  system  with  its  physical 
performance.  One  of  the  current  challenges  of  SSM  is  the  implementation  of  a tool  for 
that  purpose.  The  objective  of  this  dissertation  is  to  develop  and  validate  procedures  to 
predict  the  spatial  pattern  of  yield.  Two  procedures  deserve  special  attention:  artificial 
neural  networks  and  crop  simulation  models.  In  addition  we  aim  to  understand  the 
mechanisms  that  lead  to  spatial  yield  variation  under  soil-water  limiting  conditions. 

A field  experiment  was  carried  out  in  a 3.6  ha  com  farm  field  in  Michigan,  USA, 
for  two  years.  In  1997,  a dry  year,  plant  population  and  plant-available  soil-water,  as 
affected  by  effective  soil  depth,  soil-water  holding  limits  and  seasonal  amount  of  rainfall, 
accounted  for  as  much  as  82%  of  the  total  final  grain  yield  spatial  variability  (range  of 
6.6  - 12.0  t/ha).  In  1998,  with  average  seasonal  rainfall  amounts  and  a more  uniform 


xvm 


stand,  plant-available  soil-water  alone  was  able  to  explain  68%  of  final  grain  yield  spatial 
variability  (range  of  9.8  - 12.6  t/ha). 

Effective  soil  depth  was  more  strongly  related  to  topography  than  sand  and  clay 
content.  The  topographic  attributes  that  were  most  related  to  these  soil  properties  were 
downstream  flow  length  and  distance  to  flow  network.  Corn  growth  and  yield  were 
consistently  related  to  topography  in  both  a dry  and  a normal  rainfall  year.  These 
relations  were  stronger  for  the  dry  year,  especially  for  downstream  flow  length  and 
distance  to  flow  network. 

The  neural  network  model  that  included  all  variables  in  the  study  produced 
smaller  predictive  errors  (average  of  91 8 kg/ha).  The  agronomic  and  topographic  neural 
networks  had  RMSEP  over  1000  kg/ha.  CERES-Maize  accurately  predicted  the  spatial 
variability  of  corn  growth  and  yield  for  two  years  of  independent  data  when  accurate 
inputs  were  provided  for  soil  properties  and  plant  population  (average  RMSEP  of  501 
kg/ha). 
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CHAPTER  1 
INTRODUCTION 


What  Is  Site-Specific  Management? 

The  technological  ability  to  apply  agricultural  inputs  to  field  crops  in  a spatially 
variable  and  economically  viable  manner  has  triggered  new  opportunities  in  crop 
management.  This  new  area  of  endeavor  has  received  several  designations  and 
definitions.  Some  of  the  most  common  designations  are  precision  agriculture,  precision 
farming,  prescription  farming,  site-specific  crop  (and/or  soil)  management  and  spatially 
variable  crop  (and/or  soil)  management.  Pierce  and  Sadler  (1997)  defined  it  as  "the 
management  of  soil  and  crops  according  to  localized  conditions  within  a field"(p.xvii). 
Another  more  complete  definition  was  given  by  Robert  et  al.  (1995):  "an  information  and 
technology  based  agricultural  management  system  to  identify,  analyze,  and  manage  site- 
soil  spatial  and  temporal  variability  within  fields  for  optimum  profitability,  sustainability, 
and  protection  of  environment"(p.vii). 

Although  some  of  the  technological  tools  used  in  site-specific  management  (SSM) 
are  new,  the  overall  idea  is  not.  Farmers  and  agronomists  have  been  aware  of  the  spatial 
variability  of  soil  properties  and  crop  performance  for  long  a time  (Crawford  et  al., 

1997).  Most  farmers  know  their  best  and  worst  areas  within  each  field,  the  'wet'  and  the 
'dry'  spots  most  commonly.  Agronomists  have  composed  their  soil  samples  from 
different  locations  in  each  field  to  get  some  representation  of  existing  spatial  variability. 
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Until  recently,  however,  the  ability  to  move  from  field-averaged  to  more  precise  spatially 
variable  field  characterizations  and  input  recommendations  has  not  been  technically  and 
economically  feasible. 

Technology  in  Site-Specific  Management 

The  essence  of  SSM  is  the  increase  of  precision  with  which  usual  crop  and  soil 
management  practices  are  performed  in  order  to  apply  optimal  levels  of  inputs  in  all  field 
locations.  This  implies  spatially  variable  nutrient  application,  pesticide  application,  plant 
density,  irrigation,  tillage,  etc,  according  to  specific  characteristics  of  each  point  in  the 
field.  To  accomplish  these  tasks,  several  technological  breakthroughs  such  as  yield 
mapping,  variable  rate  technology  (VRT)  and  real-time  sensors  had  to  occur.  In  addition, 
some  existing  technologies  had  to  become  economically  viable.  Examples  of  these  are 
global  positioning  system  (GPS),  geographic  information  systems  (GIS)  and  remote 
sensing. 

Each  of  these  technologies  plays  a specific  role  in  SSM.  Yield  mapping,  real- 
time sensors  and  remote  sensing  contribute  to  a more  precise  characterization  of  soil 
properties  and  crop  performance  within  fields,  i.e.,  they  assist  variability  measurement. 
On  the  other  hand  VRT  permits  a more  precise  implementation  of  spatially  variable  input 
recommendation,  i.e.,  assists  recommendation  implementation.  GIS/GPS  works  as  the 
geo-referenced  spatial  database  manager  to  combine  data  with  agronomic  knowledge  to 
produce  spatially  variable  input  recommendations,  i.e.,  assists  knowledge  application  and 


recommendation  development. 
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The  Challenge  of  Site-Specific  Management 

The  knowledge  application  and  recommendation  development  step  constitutes  one 
of  the  biggest  challenges  of  SSM  (Schueller,  1996;  Acock  and  Pachepsky,  1997).  A 
strong  link  between  the  diverse  layers  of  data  describing  existing  variability  and  the 
consequent  spatially  variable  management  recommendation  is  currently  missing.  The 
most  attractive  procedure  is  to  apply  the  same  rules  of  whole-field  management  to  every 
single  location  in  the  field  despite  the  scale  dependence  of  those  rules  and  their 
appropriateness  (Sudduth  et  al.,  1996).  In  addition,  these  traditional  methods  do  not 
consider  layers  of  data  that  were  only  made  available  recently  (examples  of  remote 
sensing  data,  terrain  characteristics,  etc).  As  a result,  persistence  in  using  whole-field 
based  recommendation  rules  in  SSM  will  lead  to  inefficient  use  of  inputs,  a misuse  of 
data  resources  and,  ultimately,  the  failure  of  SSM  to  prove  itself  economically  viable. 

The  successful  implementation  of  SSM  depends  on  the  ability  to  relate  the  current 
status  of  a crop  system  and  its  final  physical  performance  (grain  yield  in  most  cases)  for 
every  location  of  a field.  In  some  cases  under  environmental  constraints,  it  will  also 
depend  on  how  well  conflicting  goals  can  be  managed.  For  all  cases,  a procedure  is 
needed  to  predict  the  spatial  distribution  of  the  variable(s)  of  interest  (grain  yield  and/or, 
for  example,  nitrate  leaching).  The  desirable  method  should  (1)  predict  the  spatial 
distribution  of  yield  as  early  in  the  season  as  possible  so  that  possible  repairing 
agronomic  actions  could  take  place;  (2)  use  the  minimal  number  of  variables  as  possible 
to  make  it  easy  to  use;  use  easily  available  data  to  reduce  cost;  (3)  be  portable  in  space 
and  time  to  be  usable  in  new  fields  and  in  different  years;  (4)  be  sensitive  to  the  most 
important  yield  limiting  factors  to  be  able  to  explain  variability;  (5)  incorporate 
uncertainty  in  inputs  to  estimate  confidence  limits  for  the  variables  of  interest. 
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Only  when  such  a tool  is  available  can  spatially  variable  management 
recommendations  be  done  efficiently  and  correctly. 

Objective  of  the  Dissertation 

The  specific  objective  of  this  dissertation  is  to  develop  and  validate  over  space 
and  time,  procedures  to  predict  the  spatial  pattern  of  grain  yield  in  a soil-water  limiting 
environment  with  as  many  of  the  characteristics  specified  above  as  possible.  The  focus 
on  the  water  limiting  conditions  results  from  the  fact  that  most  intensively  farmed  fields 
are  managed  to  limit  yield  reductions  due  to  nutrients,  diseases,  pests  and  weeds.  An 
additional  goal  of  the  dissertation  is  to  understand  the  mechanisms  that  lead  to  spatial 
yield  variation  under  soil-water  limiting  conditions  by  describing  and  quantifying  the 
spatial  variability  of  the  responsible  variables. 

Two  procedures  to  describe  and  predict  the  spatial  pattern  of  grain  yield  deserve 
special  attention  in  this  work:  artificial  neural  networks  and  crop  simulation  models.  The 
two  procedures  are  inherently  different  in  nature.  An  artificial  neural  network  is  a totally 
empirical  tool  whereas  crop  simulation  models  are  based  on  the  simulation  of 
physiological  and  physical  processes. 

Outline  of  the  Dissertation 

This  dissertation  assembles  five  distinct  pieces  of  research,  each  of  which  is 
presented  in  a separate  chapter.  All  chapters  are  organized  independently  so  that  they 
could  be  understood  without  knowledge  of  the  preceding  ones.  The  five  studies  are 
based  on  an  experiment  carried  out  in  a corn  farm  field  in  Michigan,  USA,  during  the 


summers  of  1997  and  1998. 
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In  Chapters  2 and  3 results  of  exploratory  analysis  of  the  data  are  presented.  The 
intent  was  to  identify,  describe  and  quantify  the  spatial  variability  of  the  relevant 
variables  causing  spatial  yield  variability.  In  the  following  three  chapters  (chapters  4,  5 
and  6)  two  distinct  tools  are  investigated  for  their  ability  to  describe  and  predict  the 
spatial  pattern  of  grain  yield. 

In  chapter  2,  the  relationships  among  seasonal  rainfall,  spatial  variability  of  soil 
profile  characteristics  and  yield  components  are  analyzed  in  detail.  Maps  of  each  of  the 
variables  are  presented  and  related  to  each  other  according  to  the  seasonal  rainfall  pattern. 
Conclusions  are  drawn  focusing  in  the  most  important  variables  related  to  spatial  yield 
variability  and  the  mechanisms  that  lead  to  it. 

In  chapter  3,  an  additional  variable  is  included  in  the  study:  topography.  The 
interest  in  this  new  variable  arrives  from  its  close  relation  to  important  hydrological 
processes  such  as  surface  and  subsurface  lateral  water  flows.  In  addition,  topography  is  a 
relatively  easy  variable  to  obtain  and  is  related  to  soil  properties  that  in  turn  are  decisive 
in  grain  yield  formation.  The  maps  of  the  several  attributes  used  to  describe  the 
topography  of  the  studied  field  are  presented,  and  the  relationships  to  soil  profile  and 
crop  variables  studied.  The  conclusions  focus  the  relationships  between  the  most 
important  topographic  attributes  and  spatial  yield  variability. 

In  chapter  4 three  artificial  neural  networks  models  are  implemented,  trained  and 
tested  for  their  ability  to  explain  and  predict  final  grain  yield  variability  across  the  field. 
The  distinction  among  models  was  due  to  the  types  of  variables  used.  The  maps  of 
predicted  final  grain  yield  are  presented  and  compared  to  the  maps  of  observed  final  grain 
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yield  in  order  to  conclude  which  was  the  most  accurate  artificial  neural  network  model 
for  spatial  yield  variability  prediction. 

In  chapter  5 a crop  simulation  model  is  used  to  describe  and  predict  growth,  final 
grain  yield  and  soil-water  content  spatial  variability  within  the  field.  The  maps  of 
predicted  final  grain  yield  are  presented  and  compared  to  the  observed  final  grain  yield. 
Conclusions  are  drawn  in  terms  of  the  overall  ability  of  the  crop  simulation  model  to 
describe  and  predict  accurately  the  observed  spatial  yield  variability. 

In  chapter  6 a procedure  to  estimate  site-specific  soil  characteristics  required  by 
crop  simulation  models  is  described  and  tested.  The  procedure  is  based  on  the  inverse 
modeling  process  using  specific  simulated  and  observed  variables  in  the  objective 
function.  Comparison  is  made  for  two  different  objective  function  variables:  the  more 
theoretically  correct  soil-water  content  time-series  and  the  more  practical  final  grain 
yield. 

Finally,  chapter  7 summarizes  the  dissertation  and  the  conclusions. 


CHAPTER  2 

SPATIAL  VARIABILITY  OF  CORN  GROWTH  AND  YIELD  DETERMINED  BY 
SOIL-WATER  AND  PLANT  POPULATION 

Introduction 

An  essential  aspect  in  site-specific  management  (SSM)  is  the  understanding  of  the 
most  important  factors  limiting  final  yield  (Mulla  and  Schepers,  1997;  Sudduth  et  al., 
1996;  Verhagen  et  ah,  1995).  Only  with  that  knowledge  can  informed  and  economically 
based  decisions  take  place.  Early  attempts  to  explain  within-field  spatial  yield  variability 
focused  on  macronutrients,  mainly  nitrogen  and  phosphorus,  and  pH  (e.g.  Everett  and 
Pierce,  1996;  Khakural  et  ah  1996;  Sudduth  et  ah,  1996;  Mallarino  et  ah,  1996;  Nolin  et 
ah,  1996).  These  studies  showed  that  those  fertility  factors  were  not  the  main  cause  for 
the  spatial  pattern  found  in  yield  maps.  They  also  concluded  that  there  was  a common- 
sense  reason  for  the  finding:  farmers  that  implemented  SSM  were  already  applying  levels 
of  those  inputs  that  were  not  limiting  yield,  and  therefore,  those  factors  were  not 
responsible  for  the  variability.  Furthermore,  due  to  its  relative  low  cost,  economical 
gains  from  spatially  variable  applications  of  those  factors  were  not  as  beneficial  as 
originally  conceived,  particularly  if  environmental  issues  were  not  to  be  considered. 

After  these  early  attempts,  research  shifted  to  other  yield  limiting  factors:  plant- 
available  soil-water,  plant  population,  weeds,  pests  and  diseases.  This  study  focuses  on 
the  first  two  aspects. 

Soil-water  can  be  a limiting  factor  to  plant  growth  and  yield  when  available  in  too 

low  or  too  high  amounts.  At  the  macroscopic  level,  soil-water  deficit  can  have  distinct 

effect  on  yield  depending  on  the  severity  and  the  timing  of  its  occurrence  in  the  crop 
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season.  Early  in  the  season,  during  vegetative  growth,  soil-water  deficit  leads  to  low  leaf 
area  index,  low  intercepted  radiation  and  low  plant  biomass  growth.  Ultimately,  this  will 
reduce  yield  by  reducing  assimilate  production  during  the  grain  filling  phase  due  to  low 
radiation  interception.  Later  in  the  season,  soil-water  deficit  can  limit  yield  by  its  effect 
on  important  yield  components  such  as  number  of  grains  per  ear  and  unit  grain  weight. 
The  dynamic  aspect  of  soil-water  deficit  impacts  on  grain  yield  highlights  the  importance 
of  a spatial-temporal  approach  when  analyzing  grain  yield  spatial  variability.  Soil-water 
excess,  either  transient  or  permanent,  can  also  limit  crop  growth  and  yield.  Its  damaging 
effect  is  due  to  lack  of  oxygen  to  the  root  system,  thus  reducing  growth  and  accelerating 
senescence.  Similar  to  soil-water  deficit,  the  timing  of  occurrence  of  excess  water 
determines  the  degree  of  damage  and  yield  loss. 

Soil-water  Availability  to  Plants 

For  corn  grown  under  non-irrigated  conditions,  plant-available  soil-water 
(PASW)  is  a common  sense  factor  for  analyzing  yield  variability.  In  these  systems, 
PASW  is  mostly  influenced  by  effective  soil  depth,  rooting  depth,  soil-water  holding 
limits  (that  determine  soil-water  holding  capacity),  position  in  the  landscape  (which 
relates  to  soil  depth  and  soil-water  holding  limits  in  conjunction  with  the  parent  material) 
and  seasonal  rainfall  pattern  and  amount. 

Sudduth  et  al.  (1995)  obtained  significant  correlation  coefficients  of  0.62  between 
corn  yield  and  depth  to  claypan  in  a field  with  effective  soil  depth  varying  from  10  to  70 
cm.  Morton  et  al.  (1999)  studied  the  effect  of  spatial  variability  of  plant-available  soil- 
water  by  deriving  linear  regression  models  between  grain  yield  and  evapotranspiration, 
deep  percolation  and  water  stress  (measured  as  percent  of  water  between  field  capacity 
and  permanent  wilting  point).  By  combining  all  variables  in  a multiple  regression 
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approach,  they  were  able  to  explain  83%  of  total  final  grain  yield  variability.  Logsdon  et 
al  (1999),  conducting  research  in  a high  water  table  soil  in  Iowa,  only  found  significant 
correlation  between  water  holding  capacity  and  corn  yield  for  field  areas  with  higher 
elevations.  The  overall  field  coefficient  of  determination  for  the  relation  between  the  two 
variables  was  13%.  Other  studies  reported  the  occurrence  of  spatial  variability  of  soil- 
water  holding  capacity  or  water  content  but  without  direct  assessment  of  its  impact  on 
crop  yields  (e.g.  Fitzjohn  et  al.,  1998;  Mulla,  1998;  Torner  and  Anderson,  1995;  Moore  et 
al.  1988;  Hanna  et  al.  1982). 

From  simulation  studies  on  wheat  in  South  Eastern  Australia,  using  CERES- 
Wheat  simulation  model  (Hoogenboom  et  al.,  1994),  Moore  and  Tyndale-Biscoe  (1999) 
concluded  that  the  strategic  importance  of  profile  soil-water  holding  capacity  was  much 
higher  than  spatially  variable  nitrogen  management  in  determining  yield  spatial 
variability. 

Paz  et  al.  (1997)  used  the  CROPGRO-Soybean  simulation  model  (Hoogenboom 
et  al.,  1994)  to  test  the  hypothesis  that  soil-water  related  stresses  were  responsible  for  the 
spatial  variability  of  grain  yield  of  a soybean  field  in  Iowa,  USA.  They  concluded  that 
variable  rooting  and  water  table  depths  explained  69%  of  the  total  variability.  The 
additional  variability  was  hypothesized  to  be  due  to  spatially  variable  plant  population, 
disease  incidence  and  rainfall/infiltration  ratio  across  the  field. 

None  of  the  above  studies  accounted  for  systematic  and  detailed  measurement  and 
analysis  of  the  effects  of  within-field  differential  PASW  on  corn  growth,  yield  and  yield 
components.  Only  such  detailed  studies  can  provide  an  understanding  of  the  spatial  and 
temporal  components  of  the  effect  of  PASW  on  final  grain  yield. 
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Plant  Population 

For  crops  with  limited  compensation  effect  (non-tillering  species),  such  as 
temperate  com,  the  effect  of  low  plant  population  (PP)  on  yield  is  dramatic.  Potential  PP 
is  determined  by  seeding  rate,  which  is  reduced,  then,  due  to  poor  germination  and/or 
emergence  thereby  producing  plant  skips.  Poor  germination  and/or  emergence  can  be 
due  to  variable  sowing  depth,  soil-water  stresses  (too  wet  or  too  dry),  existence  of 
physical  obstacles  at  the  soil  surface  (stones,  clods  or  crop  residue  patches),  pests  and 
diseases  and  herbicides  residue  (Pommel  and  Bonhomme,  1 998).  All  these  factors  are 
enhanced  by  a delay  in  germination/emergence  that  can  occur  due  to  low  soil 
temperature.  Another  source  of  PP  variability  is  the  existence  of  plant  spacing  variability 
and  "doubles".  These  sources  of  variability  are  mainly  due  to  machine/operator  factors. 
The  determination  of  the  spatial  pattern  of  PP  is  a costly  and  time-consuming  operation. 
For  this  reason,  the  PP  at  harvest  is  not  usually  considered  in  efforts  to  explain  final  grain 
yield  spatial  variability.  The  situation  might  change  in  the  short  term  with  combine- 
mounted  devices  for  measuring  corn  plant  spacing,  and,  therefore,  PP  (Easton,  1996; 
Plattnet  and  Hummel,  1996;  Birrell  and  Sudduth,  1995).  To  our  knowledge,  no  studies 
exist  on  the  quantification  of  within-field  spatial  variability  of  PP  in  corn  and  its 
contribution  to  final  yield  spatial  variability. 

Objective 

The  objective  of  this  paper  was  to  describe  and  quantify  the  effects  of  plant 
population,  effective  soil  depth,  soil-water  holding  limits  (soil-water  holding  capacity) 
and  seasonal  pattern  of  rainfall  on  the  spatial  pattern  of  corn  growth,  yield  and  yield 
components  (grain  number  per  ear,  grain  number  per  unit  area  and  unit  grain  weight). 

The  general  hypothesis  is  that  plant-available  soil-water  and  plant  population  vary 
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considerably  in  space  and  account  for  most  of  the  within-field  spatial  variability  found  in 
a non-irrigated  cornfield.  The  present  results  constitute  part  of  an  analysis  of  a broader 
effort  to  understand,  describe  and  predict  spatial,  within-field,  corn  yield  patterns. 

Materials  and  Methods 

Field  Experiment 

A field  experiment  was  carried  out  in  a 3.6  ha  corn  (Zea  Mays  L.)  farm  field 
located  in  Michigan,  USA  (42°  50’N,  83°  58’ W)  during  the  crop  seasons  of  1997  and 
1998.  The  field  experiment  had  no  imposed  levels  of  any  kind  of  inputs.  The  experiment 
consisted  of  observing  variables  of  interest  in  a uniformly  managed  field.  The  farmer 
managed  the  field  exactly  in  the  same  way  as  he  had  been  doing  for  over  10  years.  The 
field  was  chosen  from  among  others  for  its  characteristic  spatial  variability  of  final  grain 
yield  according  to  the  farmer. 

The  uniform  management  consisted  of  a no-tillage,  non-irrigated  corn-soybean 
rotation  with  two  nitrogen  applications  and  one  pre-emergence  herbicide  application. 

The  only  interference  in  the  "normal"  management  of  the  field  was  to  break  the  crop 
rotation  and  plant  corn  for  the  two  consecutive  years  of  the  experiment.  The  farmer 
planted  the  crop  on  April  29  (DOY  1 19)  in  1997  and  on  May  1 1 (DOY  131)  in  1998. 

The  sites  for  observing  soil  and  crop  variables  (total  of  43)  were  spatially 
distributed  in  a grid  of  approximately  30.4  m (100  ft).  The  grid  was  considered  a good 
representation  of  all  combinations  of  soil  and  crop  conditions  important  in  the  study 
(Figure  1). 

Effective  soil  depth  (simply  designated  as  soil  depth)  was  determined  by  the  depth 
of  an  existing  hard  clay  layer,  determined  through  direct  observation  using  a probe 
(hydraulic  auger).  Soil  particle  size  analysis  was  carried  out  for  each  15-cm  increment 
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layer  down  to  the  constraining  clay  layer  using  the  standard  hydrometer  method  (Gee  and 
Bauder,  1986).  Average  clay  and  sand  contents  (simply  designated  as  clay  and  sand)  of 
the  top  seven  layers  (0-105  cm)  were  used  as  direct  indicator  variables  of  the  soil-water 
holding  capacity,  i.e.,  the  difference  between  field  capacity  and  permanent  wilting  point 
of  each  layer. 

Plant  dry  weight  per  unit  area  and  leaf  area  index  (LAI)  were  determined  at 
anthesis  (DOY  21 1 in  1997  and  207  in  1998)  for  three  individual  plants  in  each  of  the  43 
sites.  The  limited  size  of  the  sample  was  due  to  farmer  constraints.  The  area  of  the 
leaves  of  each  individual  plant  was  measured  using  a LICOR  leaf  area  meter  instrument. 

Grain  yield  and  its  components,  namely  grain  number  per  ear  (G#EAR),  grain 
number  per  unit  area  (G#AREA)  and  unit  grain  weight  (UGW)  as  well  as  PP  were 
determined  at  maturity  (DOY  289  in  1997  and  288  in  1998)  on  a 3 x 2 m2  plot  at  all  43 
sites.  PP  was  determined  by  measuring  individual  plant  spacing.  Since  every  plant  had 
only  one  ear,  the  yield  component  'number  of  ears  per  plant'  was  the  same  as  PP  and  thus 
was  not  considered  separately  in  the  analysis.  No  sub-sampling  technique  was  used.  All 
ears  and  grains  were  used  to  calculate  G#EAR,  G#AREA  and  UGW.  The  number  and 
dry  weight  of  the  grains  of  each  individual  ear  were  determined  after  drying  at  70  C until 
constant  weight  (6  to  8 days). 

A tipping  bucket  precipitation  sensor  (1-mm  accuracy)  that  was  located  200 
meters  from  the  border  of  the  field  recorded  rainfall  data  daily.  The  reset  time  for  daily 
amounts  was  done  at  midnight.  Maximum  and  minimum  daily  temperature  and  solar 
radiation  were  also  recorded  at  the  same  site. 

According  to  the  soil  survey  of  Shiawassee  County  (Threlkeld  and  Feenstra, 

1974)  the  soils  in  the  field  belong  to  a Miami-Conover-Brookston  Association 
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characterized  as  well  drained  to  poorly  drained,  loamy  soils  on  till  plains  and  moraines. 
Three  series  cover  97.3%  of  the  field's  area:  Breckenridge  (Bt)  with  46.2%  of  the  area, 
Macomb  (MaA)  with  39.2%  and  Metamora  (MsA)  with  12.0%.  The  remaining  area  is 
occupied  by  the  series  Conover  (CtA)  and  Wasepi  (WeA)  with  0.9%  and  1.8%  of  the 
total  area.  A map  of  the  soils  in  the  field  is  shown  in  Figure  1 . More  complete 
classification  of  the  soil  series  present  in  the  field  is  given  in  Table  1 . 

The  soils  in  the  field  have  a profile  drainage  constraint  caused  by  a clay  layer 
located  at  variable  depths.  It  is  usually  shallower  in  zones  of  higher  elevation  and  deeper 
in  depressions.  Typically,  this  layer  causes  a perched  water  table  that  is  present  in  the 
profile  at  the  beginning  of  spring,  due  to  snow  melt,  and  at  the  end  of  the  summer  in  high 
rainfall  years  (summer  storms).  In  the  middle  of  the  season,  the  water  table  disappears 
once  depleted  by  roots.  Therefore,  the  water  table  may  play  a very  important  role  in 
supplying  water  during  critical  yield  forming  phases:  flowering  and  grain  filling.  In  years 
with  higher  rainfall,  the  water  table  is  present  throughout  the  season  at  higher  and  more 
variable  depth  causing  excess  water  stress  to  plants. 


CTA  WEA  BT 


Figure  1 . Map  showing  the  location  of  the  43  sampling  points  used  in  the  study  and  the 
soil  series  designated  by  2 or  3-letter  codes  according  to  the  order  two  soil  survey  at 
1 :20000  scale  (Threlkeld  and  Feenstra,  1974).  (Bt  = Breckenridge;  MaA  = Macomb; 
MsA  = Metamora;  WeA  = Wasepi;  and  CtA  = Conover). 
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Table  1 . Classification  of  soil  series  present  at  the  experimental  field. 


Series 

Symbol 

Family 

Subgroup 

Order 

Breckenridge 

Bt 

Coarse-loamy,  mixed,  nonacid, 
frigid 

Mollic  Haplaquepts 

Inceptisols 

Macomb 

MaA 

Fine-loam,  mixed,  mesic 

Udollic  Ochraqualfs 

Alfisols 

Metamora 

MsA 

Fine-loam,  mixed,  mesic 

Udollic  Ochraqualfs 

Alfisols 

Conover 

CtA 

Fine-loam,  mixed,  mesic 

Udollic  Ochraqualfs 

Alfisols 

Wasepi 

WeA 

Coarse-loamy,  mixed,  mesic 

Aquollic  Hapludalfs 

Alfisols 

Source:  Threlkeld  and  Feenstra,  1974. 


Statistical  Analysis 

Maps  showing  the  spatial  pattern  of  all  variables  were  developed.  Semivariance 
analysis  and  kriging  were  used  for  interpolation  proposes  (Issacks  and  Srivastava,  1 989). 
The  variograms  were  derived  by  plotting  the  semivariance  of  the  data  as  a function  of 
vector  distance  before  being  modeled.  The  software  Variowin  (Pannatier,  1996)  and 
Surfer  (Golden  Software,  CO)  were  used  for  the  analyses. 

The  relationships  between  independent  and  dependent  variables  were  studied 
using  correlation  analysis  for  simple  effects  and  forward-stepwise  multiple  regression 
analysis  for  multiple  combined  effects.  To  test  the  significance  of  the  effects  of  the 
independent  variables,  as  well  as  their  interactions,  on  the  dependent  variables,  analysis 
of  variance  was  used  assuming  an  unbalanced  fixed  factor  factorial  design.  For  this 
purpose,  each  independent  variable  was  categorized  into  two  classes  based  on  the 


median. 


Cumulative  rainfall  (mm) 
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Figure  2.  Cumulative  rainfall  for  1997  and  1998  years  and  for  the  corn  growing  seasons 
in  these  two  years  measured  200  m away  from  the  field  border. 
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Results  and  Discussion 

Independent  Variables 

Figure  2 shows  cumulative  rainfall  for  1997  and  1998  and  for  the  corn  growing 
seasons  in  these  years.  1 997  was  clearly  drier  than  1 998  during  the  crop  season  and  the 
whole  year.  The  final  values  for  1997  and  1998  were  529.1  and  849.1mm  for  the  year 
and  225.9  mm  and  436.1  mm  (difference  was  210.2  mm)  for  the  crop  season.  This 
difference  made  these  two  years  satisfactory  to  study  the  effects  of  seasonal  pattern  of 
rainfall  on  grain  yield  and  its  spatial  pattern.  In  1997  there  were  periods  of  50,  27,  14  and 
1 1 days  with  no  rainfall  while  in  1998  the  longest  dry  periods  were  8,  7 and  6 days. 

From  a daily  weather  data  series  (1979-1988)  of  the  same  location,  the  average  seasonal 
rainfall  amount  is  459  mm  with  minimum  and  maximum  values  of  313  and  701  mm. 
Considering  these  statistics,  1 998  can  be  considered  an  average  year  in  terms  of  seasonal 
rainfall  whereas  1997  was  a dry  one  (the  driest  in  the  series). 

The  longest  period  of  rainfall  absence  in  1 997  followed  planting.  This  dry  period 
was  accompanied  by  low  air  temperature  (Figure  3 above)  causing  crop  emergence  to  be 
delayed.  In  1997  emergence  occurred  30  days  after  planting  versus  only  8 days  in  1998. 
Consequently,  PP  at  maturity  was  more  variable  in  1997  than  in  1998. 


Daily  solar  radiation  (MJ/m2)  Daily  mean  temperature  (C) 
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Figure  3.  Daily  mean  air  temperature  (above)  and  solar  radiation  during  1997  and  1998 
com  growing  seasons. 
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Table  2 summarizes  the  descriptive  statistics  of  all  the  independent  variables 
considered  in  the  study.  Figure  4 shows  the  maps  of  sand  and  clay  in  the  field.  All 
combinations  of  point  values  of  sand  and  clay  classified  the  textures  as  loam,  sandy-loam 
and  silt-loam  types.  The  values  of  sand  and  clay  varied  between  34.3-67.7%  and  12.8- 
20.9%,  respectively.  The  values  of  those  two  variables  were  negatively  (r2  = 0.78) 
correlated  for  the  43  sites  used  in  the  study  indicating  that  only  one  of  these  two  variables 
can  be  used  as  an  indicator  for  soil-water  holding  capacity. 

The  map  of  soil  depth  is  shown  in  Figure  5.  The  point  values  varied  between  105 
and  150  cm.  Most  of  the  sites  had  a soil  depth  greater  than  or  equal  to  120  cm  with 
almost  even  distribution  in  the  120,  135  and  150  cm  depths  (16,  14  and  1 1 sites 
respectively).  Only  2 sites  had  a soil  depth  of  105  cm. 

According  to  the  pedotransfer  function  developed  by  Saxton  et  al.  (1986)  and 
validated  by  Kern  (1995),  for  the  textures  at  all  sites  and  all  depths,  sand  (clay)  is 
negatively  (positively)  correlated  with  soil-water  holding  capacity  (varying  from  0.080  to 
0. 1 70  cm3  cm'3).  The  values  of  soil-water  holding  capacity  given  by  the  same 
pedotransfer  function  for  all  the  43  profiles  averaged  164  mm  with  extreme  values  of  105 
and  200  mm.  The  main  driver  for  the  spatial  variability  found  in  the  textures  and  soil 
depth,  and  consequently  in  soil-water  holding  capacity,  is  topography.  The  inclusion  of 
topographic  attributes  as  independent  variables  is  studied  in  chapter  3. 
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Table  2 - Mean  values  and  range  of  variation  for  independent  and  dependent  variables 
in  the  study.  


Unit 

Mean 

Minimum 

Maximum 

Standard 

deviation 

Sand  content 

% 

46.2 

34.3 

67.7 

8.0 

Clay  content 

% 

17.6 

12.8 

20.9 

2.1 

Effective  soil  depth 

cm 

135 

105 

150 

13.3 

97 

#/m2 

8.4 

5.7 

10.3 

1.2 

Plant  population 

98 

#/m2 

9.0 

7.4 

10.3 

0.7 

97 

Kg/ha 

9571 

6620 

11963 

1324 

Grain  yield 

98 

Kg/ha 

11331 

9834 

12629 

760.0 

97 

- 

365 

274 

501 

54.6 

Grain  number  per  ear 

98 

- 

463 

387 

553 

42.3 

97 

mg 

303.4 

266 

329 

20 

Unit  grain  weight 

98 

mg 

271.2 

238 

286 

12.0 

97 

ft/m2 

3007 

2614 

3593 

219.1 

Grain  number  per  unit  area 

98 

ft/m2 

4172 

3631 

4698 

271.2 

97 

tt/m2 

8.4 

5.7 

10.3 

1.2 

Plant  population 

98 

ft/m2 

9.0 

7.4 

10.3 

0.7 

97 

- 

3.9 

2.6 

5.4 

0.74 

LAI  at  anthesis 

98 

- 

5.1 

3.7 

6.4 

0.66 
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Figure  4.  Maps  of  the  average  sand  (above)  and  clay  (below)  content  (%)  for  the  top  105 
cm  of  soil. 


Figure  5.  Map  of  the  effective  soil  depth  (cm). 
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As  mentioned  previously,  PP  was  more  variable  in  1997  than  in  1998.  The 
maximum  PP  values  were  approximately  the  same  for  both  years  (10.3  #/m'2)  whereas 
minimum  values  were  much  smaller  for  1997  (5.4  #/m'2  in  1997  and  7.4  #/m"2  in  1998). 
The  map  of  PP  for  both  years  is  shown  in  Figure  6.  The  maps  show  little  similarity 
between  years.  The  zone  of  greater  PP  reduction  in  1997  corresponds  to  an  area  of 
higher  sand  and  shallower  soil  profiles.  This  suggests  that  the  delayed  emergence  was 
not  the  only  cause  for  lower  PP  but  a combined  effect  with  low  soil-water  content. 
Dependent  Variables 

Table  2 summarizes  the  descriptive  statistics  all  of  the  dependent  variables 
considered  in  the  study.  Figure  7 shows  the  maps  of  grain  yield  for  1997  and  1998  and 
Figure  8 shows  the  respective  area-based  cumulative  frequency  for  yield.  Grain  yield 
was  lower  in  1997  than  in  1998  in  the  entire  field  (mean  of  9571  vs.  11331  kg/ha, 

F=161 .99***).  The  spatial  pattern  was  similar  between  years.  The  difference  between 
years  varied  from  around  3000  kg/ha  in  lower  yielding  areas  to  about  500  kg/ha  in  higher 
yielding  areas.  The  lower  yielding  zones  corresponded  to  areas  of  shallower  soil  profiles 
with  higher  sand  in  both  years.  In  1997,  lower  yielding  zones  also  corresponded  to  areas 
of  lower  PP.  These  observations  are  consistent,  in  part,  with  the  analysis  of  variance 
where  significant  differences  were  found  for  the  simple  effects  of  soil  depth  (F= 
125.41***)  and  for  the  interaction  between  year  and  the  effect  of  PP  (F=8.73**).  No 
other  effects  were  significant,  including  that  of  sand.  The  coefficients  of  determination 
between  the  independent  variables  and  grain  yield  are  shown  in  Table  3 and  Table  4. 
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Figure  6.  Map  of  plant  population  (#/m2)  for  1997  (above)  and  1998  (below)  com  growth 


seasons. 


Figure  7.  Map  of  com  grain  yield  (kg/ha)  for  1997  (above)  and  1998  (below)  growing 
seasons. 
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Figure  8.  Area-based  cumulative  frequency  of  com  grain  yield  for  1997  and  1998 
growing  season. 


Table  3.  Coefficients  of  determination  between  independent  and  dependent  variables  for 
1997.  Non  significant  relationships  are  marked  with  n.s.  The  significant  correlation  are 
classified  according  to  the  level  of  significance  (*  p<0.05,  **  p<0.01,  ***  p<0.001). 


Sand 

Clay 

Soil  depth 

Plant 

Population 

Grain  yield 

(-)  0.46*** 

(+)  0.48*** 

(+)  0.79*** 

(+)  0.56*** 

Grain  number  per  ear 

(+)  0.23** 

(-)  0.20** 

(.)  o.49*** 

(-)  0.73*** 

Unit  grain  weight 

(-)  0.50*** 

(+)  0.41*** 

(+)  0.74*** 

(+)  0.49*** 

Grain  number  per  unit 
area 

(-)  0.12* 

(+)0.11* 

n.s. 

(+)  0.12* 

Plant  weight  at 
anthesis  per  unit  area 

(-)  0.65*** 

(+)  0.56*** 

(+)  0.43*** 

(+)  0.34** 

LAI  at  anthesis 

(-)  0.33*** 

(+)  0.32** 

(+)  0.28** 

(+)  0.54*** 
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Table  4.  Coefficients  of  determination  between  independent  and  dependent  variables  for 
1998.  Non  significant  relationships  are  marked  with  n.s.  The  significant  correlations  are 
classified  according  to  the  level  of  significance  (*  p<0.05,  **  p<0.01,  ***  p<0.001). 


Sand 

Clay 

Soil  depth 

Plant 

Population 

Grain  yield 

(-)  0.41*** 

(+)  0.36** 

(+)  0.56*** 

n.s. 

Grain  number  per  ear 

n.s. 

n.s 

n.s. 

(-)  0.36** 

Unit  grain  weight 

(-)  0.55*** 

(+)  o.44*** 

(+)  0.59*** 

n.s. 

Grain  number  per  unit 
area 

(-)  0.18* 

(+)  0.21** 

n.s. 

n.s. 

Plant  weight  at  anthesis 
per  unit  area 

(-)  0.28** 

(+)  0.27** 

(+)  0.13* 

(+)  0.27** 

LAI  at  anthesis 

n.s. 

n.s. 

n.s. 

(+)  0.12* 

Figure  9 shows  the  maps  of  G#EAR  for  1997  and  1998.  G#EAR  had  a lower 
mean  in  1997  than  in  1998  (365  vs.  463,  F=196.36***).  The  spatial  pattern  was  very 
distinct  for  both  years.  Areas  with  higher  values  of  G#EAR  showed  some  degree  of 
match  between  years  but  there  was  no  relationship  for  most  of  the  field.  In  these  areas 
with  higher  G#EAR,  average  differences  between  years  of  up  to  150  grains  per  ear  were 
found.  The  analysis  of  variance  showed  no  significant  effect  of  sand  and  two  significant 
interactions:  between  PP  and  year  (F=T9.52***)  and  between  soil  depth  and  year 
(F=5.90*).  In  1997,  high  G#EAR  zones  corresponded  to  areas  of  lower  PP  and  shallower 
profiles.  This  indicates  the  existence  of  some  degree  of  compensation  between  PP  and 
G#EAR  in  1997  when  rainfall  was  more  limiting.  In  addition,  shallower  profiles  have  a 
perched  water  table  nearer  the  surface  where  roots  can  reach  it  more  rapidly.  In  1998, 
with  apparently  no  water  stress  effects  (at  least  up  to  grain  number  setting  stage),  G#EAR 
was  not  correlated  with  any  variable  included  in  the  study.  The  coefficients  of 
determination  between  these  variables  and  G#EAR  are  shown  in  Table  3 and  Table  4. 
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Figure  9.  Map  of  grain  number  per  ear  for  1997  (above)  and  1998  (below)  com  growing 
seasons. 

Figure  10  shows  the  maps  of  UGW  for  1997  and  1998.  UGW  had  a higher  mean 
in  1997  than  in  1998  (303.4  vs.  271.2  mg,  F=344.45***).  The  difference  was  due  to  the 
compensation  between  number  of  grain  and  weight.  The  map  had  some  similarities 
between  both  years.  The  differences  between  years  were  higher  for  zones  with  high 
UGW,  reaching  values  of  40  mg.  The  analysis  of  variance  showed  significant  simple 
effects  of  soil  depth  (236.89***)  and  sand  (4.22*)  and  interaction  between  year  and  the 
effect  of  PP  (6.32*).  No  other  effects  were  significant.  Higher  UGW  zones 
corresponded  to  areas  of  deeper  profiles  with  lower  sand.  In  1997,  high  UGW  zones  also 
corresponded  to  areas  with  higher  PP.  The  relation  between  PP  and  UGW  is  explained 
by  the  intermediate  relationships  with  G#EAR.  The  coefficients  of  determination 
between  these  variables  and  UGW  are  shown  in  Table  3 and  Table  4. 
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Figure  10.  Map  of  unit  grain  weight  for  1997  (above)  and  1998  (below)  com  growing 
seasons. 

Figure  1 1 shows  the  maps  of  G#AREA.  G#AREA  had  a higher  mean  in  1998 
than  in  1997  (3007  vs.  4172  m'2,  F=6 10.42***).  The  map  shows  no  similarity  between 
the  patterns  of  variation.  This  is  not  surprising  since  this  yield  component  depends  on 
two  other  variables  (PP  and  G#EAR),  which,  themselves,  do  not  show  consistency 
among  years.  The  compensation  effect  between  PP  and  G#EAR  has  led  to  the  lower 
degree  of  variability  of  G#AREA  within  each  year.  As  a result,  UGW  played  a more 
important  role  in  defining  the  spatial  pattern  of  final  grain  yield  (especially  in  1 997). 
G#AREA  was  relatively  more  important  in  defining  the  magnitude  of  grain  yield  than  its 
spatial  pattern.  Only  the  year  had  a significant  effect  on  G#AREA.  The  coefficients  of 
determination  between  these  variables  and  G#AREA  are  shown  in  Table  3 and  Table  4. 
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Figure  11.  Map  of  grain  number  per  unit  area  (#/m2)  for  1997  (above)  and  1998  (below) 
com  growing  seasons. 

Vegetative  growth  was  also  reduced  in  1997;  the  amount  of  plant  weight  at 
anthesis  (1997  average  of  5290  kg/ha  vs.  1998  average  of  7233  kg/ha,  F=181.61***)  and 
leaf  area  index  at  anthesis  (1997  average  of  3.95  vs.  1988  average  of  5.09,  F=70.92***) 
were  smaller.  Neither  map  showed  great  similarities  between  years  (Figure  12  and 
Figure  13).  The  differences  between  years  were  higher  for  areas  with  less  plant  weight 
and  LAI  (almost  less  2000  kg/ha  of  plant  weight  and  1 .5  LAI  units).  The  analysis  of 
variance  showed  significant  interactions  between  year  and  soil  depth  (F=31.75***),  sand 
(F=5.86*)  and  PP  (F=6.53*)  for  plant  weight  at  anthesis.  For  leaf  area  index,  only  the 
interaction  between  year  and  PP  was  significant  (F=6.33*).  This  means  that  in  1997, 
zones  with  higher  plant  weight  corresponded  to  deeper  profiles,  higher  PP  and  lower 
sand.  For  the  same  year,  zones  with  higher  LAI  corresponded  to  areas  with  higher  PP. 
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Figure  12.  Map  of  leaf  area  index  at  anthesis  for  1997  (above)  and  1998  (below)  com 
growing  seasons. 


These  results  emphasize  the  complexity  of  the  yield  forming  process,  its 
interactions  with  the  environment,  and  the  resulting  spatial  yield  variability.  In  1 997,  the 
driest  year  in  1 0 years,  a period  of  uncharacteristically  low  mean  temperature  and  no 
rainfall  delayed  seed  emergence.  These  events  led  to  high  spatial  variability  of  PP  that 
had  lower  values  in  zones  of  shallower  soils  and  higher  sand  content  (terrain  tops).  In 
1998,  PP  was  much  more  homogenous.  Up  to  flowering,  plant  growth  per  unit  area  was 
reduced  in  1 997  when  compared  to  1 998  due  to  lower  PP  and,  possibly,  some  water 
stress  before  roots  reached  deeper  layers.  This  effect  was  also  crucial  in  the 
determination  of  G#EAR  in  1997.  Lower  PP  in  shallower  profiles  led  to  higher  grain 
number  per  ears  and  vice-versa.  Shallower  profiles  had  the  hard-clay  layer  closer  to  the 
surface  and,  therefore,  any  perched  water  table  present  is  depleted  sooner  by  roots.  In 
1 998,  apparently  no  serious  limitation  was  imposed  on  the  number  of  grains  per  ear.  Due 
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to  the  compensation  effect  between  G#EAR  and  PP,  the  G#AREA  did  not  show  a lot  of 
spatial  variability.  Nevertheless,  it  varied  greatly  from  one  year  to  the  other,  resulting  in 
totally  different  magnitudes  for  final  grain  yield. 


Figure  13.  Map  of  the  plant  weight  at  anthesis  (kg/ha)  for  1997  (above)  and  1998 
(below)  com  growing  seasons. 

The  UGW  maps  showed  great  spatial  variability  for  both  years.  In  general,  UGW 
was  lower  in  1998  than  in  1997.  This  was  due  to  the  compensation  between  number  of 
grain  and  weight.  The  spatial  variability  within  each  year  was  probably  more  related  to 
soil  depth  variation.  Perched  water  tables  close  to  the  surface  were  responsible,  in  part, 
for  the  higher  G#EAR  of  shallower  profile  zones.  For  the  UGW  the  result  was  reversed: 
deeper  profiles  had  their  perched  water  tables  depleted  by  roots  later  in  the  season  (during 
grain  filling)  and,  therefore,  were  able  to  provide  a non-water  stress  environment  during 
grain  filling  phase.  The  result  was  higher  UGW  for  deeper  profile  zones  and  vice-versa. 
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Finally,  the  yield  map  was  the  integration  of  all  these  relations.  1998  had  higher  final 
grain  yield  due  mainly  to  the  higher  G#AREA.  1 997  had  higher  spatial  variability  of 
final  grain  yield  due  mainly  to  more  variable  UGW. 

Independent  Variables  Effect:  Regression  Analysis 

Multiple  regression  linear  equations  were  first  obtained  for  1997  and  1998 
independently.  The  purpose  was  to  test  the  combined  linear  effects  of  the  different 
independent  variables  on  final  grain  yield  and  to  assess  how  much  variability  those 
variables  could  explain  yield  for  the  year  not  used  in  the  regression  analysis.  Because  of 
its  correlation  to  sand,  clay  was  not  selected  for  any  of  the  models.  For  1997,  all  3 other 
independent  variables  were  included:  soil  depth,  sand  and  PP.  This  model  resulted  in  a 
very  good  root  mean  square  error  of  fit  of  470  kg/ha  and  coefficient  of  determination  of 
82%  (F=98.26***).  This  means  that  these  3 variables  explained  82%  of  the  total  spatial 
variability  found  in  the  1 997  grain  yield  map.  In  relative  terms,  soil  depth  was  the  most 
important  variable  (measured  by  the  magnitude  of  the  standardized  regression 
coefficients),  followed  by  PP  and  sand.  In  1998,  PP  did  not  contribute  significantly  to 
explain  grain  yield  variability.  Only  sand  and  soil  depth  were  included  in  the  regression, 
which  was  able  to  account  for  68%  (F=45.27***)  of  total  variability.  As  for  the  1997 
model,  in  1998,  soil  depth  contributed  more  to  the  regression  than  did  sand.  The  root 
mean  square  error  of  fit  was  similar  to  that  of  1 997  (43 1 kg/ha).  The  final  fits  for  these 
two  linear  models  are  shown  in  Figure  14.  Accounting  for  what  has  been  said  previously, 
these  results  make  sense.  PP  was  not  a limiting  factor  in  1998  in  contrast  to  1997; 
therefore,  it  was  not  incorporated  in  the  1998  regression  model.  1997  was  a drier  year 
than  1998;  therefore,  the  portion  of  total  grain  yield  variability  explained  by  plant- 
available  soil-water  indicators  (soil  depth  and  sand)  was  higher  in  1997. 
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Figure  15  shows  the  predicted  versus  observed  grain  yield  when  the  single-year 
model  was  used  to  predict  the  other  year's  grain  yield,  i.e.,  1997  regression  model  used  to 
estimate  1 998  grain  yields  and  vice-versa.  Although  the  individual  year  models  were 
able  to  explain  the  variability  found  in  final  grain  yield,  they  performed  poorly  when 
tested  for  prediction  in  an  independent  data  set.  This  result,  again,  stresses  the  peculiarity 
of  each  individual  year's  data  set  in  terms  of  limiting  factors  and  the  inability  of  these 
multiple  linear  regression  models  to  account  for  the  complex  interactions  present  in  a 
yield  formation  process. 


Observed  Grain  Yield  (Kg/ha) 


Figure  14.  Regressed  versus  observed  grain  yield  using  individual  year  regression 
curves. 
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Figure  15.  Predicted  versus  observed  grain  yield  using  crossed  regressions  curves 
(predict  grain  yield  of  1998  with  a regression  curve  obtained  from  1997  data  and  vice- 
versa). 


Conclusions 

Plant-available  soil-water  as  affected  by  effective  soil  depth,  soil-water  holding 
limits  and  seasonal  amount  of  rainfall,  was  the  single  most  important  factor  responsible 
for  the  spatial  pattern  of  corn  growth,  yield  and  yield  components  in  the  studied  non- 
irrigated  field.  In  the  year  of  unfavorable  conditions  for  seed  germination/emergence  and 
limited  seasonal  rainfall,  plant  population  also  played  a very  important  role  in  the  spatial 
yield  variability.  These  two  factors  accounted  for  as  much  as  82%  of  the  total  final  grain 
yield  spatial  variability  in  1997  (range  of  6.6  - 12.0  t/ha).  In  the  year  of  more  uniform 
stand  and  average  seasonal  rainfall  amounts,  plant-available  soil-water  alone  was  able  to 
explain  68%  of  final  grain  yield  spatial  variability  (range  of  9.8  - 12.6  t/ha).  In  both 
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situations  effective  soil  depth  was  more  important  than  sand  content,  or  soil-water 
holding  capacity,  in  causing  spatial  variability  of  grain  yield.  Unit  grain  weight  was  more 
important  in  defining  the  spatial  pattern  of  yield  variability  whereas  the  number  of  grains 
per  unit  area  was  more  important  in  defining  the  magnitude  of  total  field  yield.  Multiple 
linear  regression  models  cannot  assess  the  complex  interactions  among  the  different 
environmental  forces  and  crop  growth  and  yield  in  such  systems.  Other  more  complex 
prediction  tools  will  have  to  be  used. 


CHAPTER  3 

WITHIN  FIELD  SPATIAL  VARIABILITY  OF  CORN  YIELD  AND  SOIL 
PROPERTIES  IN  RELATION  TO  TOPOGRAPHIC  TERRAIN  ATTRIBUTES 


Introduction 

In  an  effort  to  characterize  and  understand  existing  spatial  variability  of  crop  and 
soil  variables,  SSM  researchers  have  been  looking  into  new  sources  of  data.  The  two 
most  important  examples  of  these  new  sources  are  probably  remote  sensing  and 
topography.  Both  sources  derive  from  relatively  old  disciplines  that  only  recently  have 
been  used  to  characterize  and  understand  the  spatial  variability  of  crop  and  soil  variables. 
This  study  focuses  on  topographic  data.  Topography  is  described  by  several  indices 
designated  as  topographic  terrain  attributes. 

Topographic  Terrain  Attributes 

Topographic  terrain  attributes  describe  the  topographic  characteristics  of  each 
location  based  on  its  own  properties  or  those  of  their  neighbors.  The  interest  is  not  only 
for  sampling  scheme  delineation  (Frazen  et  al.  1999),  but  also  for  prediction  of  spatial 
grain  yield  pattern  and,  therefore,  to  help  in  delineating  areas  with  similar  yield  potential 
(MacMillan  et  al.  1999).  The  use  of  topographic  terrain  attributes  in  trying  to  explain  and 
predict  spatial  yield  pattern  is  based  on  the  assumption  that  they  reflect  the  effects  of 
lateral  water  movement  on  site-specific  soil  properties,  as  well  as  on  soil-water  content. 

In  fact,  any  landscape  with  an  elevation  gradient  will  have  a lateral  component  added  to 
the  vertical  water  movement.  The  lateral  water  movement  can  transport  soil  particles 
resulting  in  deeper  soils  with  finer  textures  at  lower  elevation  areas.  If  no  soil  particles 
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are  transported  along,  the  lateral  water  movement  will  simply  create  a lateral  gradient  in 
soil-water  content.  Hall  and  Olson  (1991)  provide  a detailed  analysis  of  the  qualitative 
effects  of  landscape  morphology  on  soil  physical  and  chemical  properties  and  soil  water 
movement.  The  effect  of  lateral  water  movement,  at  the  surface  or  sub-surface,  on  site- 
specific  soil  properties,  such  as  clay  and  sand  contents  and  soil  depth,  occurs  over  a 
geological  time  scale.  In  contrast,  its  effect  on  the  soil-water  can  take  place  in  a matter  of 
hours,  days  or  weeks.  The  relative  importance  of  each  component  is  landscape-specific. 

Topographic  terrain  attributes  have  been  related  to  several  soil  properties  such  as 
organic  matter  content,  A-horizon  thickness,  degree  of  profile  development,  pH,  etc,  as 
well  as  soil  erosion  (Moore  et  al.,  1993;  Bell  et  al.  1995).  Other  studies  relate  them  to 
plant-available  soil-water  and  final  yield.  Denholm  et  al.  (1993)  showed  evidence  that 
corn  yield  was  negatively  correlated  with  "slope  position  and  elevation"  but  positively 
correlated  with  "percent  clay  in  the  Ap  horizon."  Their  data  set  was  obtained  from  a 
three-year  continuous-corn  field  with  loamy  sand  and  sandy  textured  soil  in  Ontario, 
Canada  where  the  topographic  relief  was  about  lm.  Nolin  et  al.  (1996),  who  also  found 
significant  correlation  between  altitude  and  crop  and  soil  related  variables,  provided 
quantification.  In  a field  with  an  elevation  relief  of  less  than  5 m,  they  obtained 
correlation  coefficients  of -0.41,  -0.30,  0.37  and  -0.25  between  altitude  and  thickness  of 
Ap  horizon,  clay  percent,  sand  percent  and  com  yield,  respectively.  Torner  et  al. 

(1995a,b)  were  able  to  explain  72%  of  soil-water  holding  capacity  of  the  top  1 .7  m of  soil 
with  elevation,  slope  and  plan  curvature  on  a sandy  textured  field  with  a topographic 
relief  of  about  3 m.  In  the  same  study,  the  thickness  of  the  top  layer  was  also  well 
correlated  with  elevation.  Timlin  et  al.  (1999)  studied  the  effects  of  weather-topography 
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interactions  on  corn  yield  on  well-drained  glacial  tills  with  topographic  relief  of  15  m. 
They  found  significant  correlation  between  yield  and  plan  curvature  (higher  yield  in 
concave  areas)  for  normal  and  dry  years  whereas  no  correlation  was  found  in  wetter 
years. 

Topographic  terrain  attributes  can  be  classified  as  primary  and  secondary 
according  to  its  derivation  (Bell  et  ah,  1995).  Primary  attributes  are  directly  derived  from 
elevation  data  and  the  secondary  attributes  are  derived  from  primary  attributes.  Other 
possible  classifications  might  be  between  the  attributes  that  are  of  broad  scope  and  the 
ones  that  are  hydrology-based  or  between  topographic  attributes  that  are  neighborhood- 
based  or  flow-based.  The  neighborhood-based  approach  assumes  an  area  of  homogeneity 
outside  the  calculation  neighborhood.  In  contrast,  flow-based  attributes  of  a single  cell 
account  for  the  characteristics  of  all  cells  flowing  to  it.  Moore  et  al.  (1993)  presented  a 
complete  table  of  the  most  important  primary  topographic  attributes  with  relevance  to 
hydrological  and  soil  applications.  Slope,  specific  catchment  area  and  plan  curvature 
were  identified  as  topographic  attributes  having  direct  effect  on  site-specific  soil-water 
content.  On  the  other  hand,  slope,  upslope  slope,  profile  and  plan  curvatures  were  related 
to  runoff  or  subsurface  flow  rate  and  volume. 

Objective 

The  objective  of  this  paper  was  to  investigate  the  relationships  among  topographic 
terrain  attributes  and  the  spatial  pattern  of  corn  growth,  yield  and  yield  components  (grain 
number  per  ear,  grain  number  per  unit  area,  unit  grain  weight,  plant  population)  and  soil 
properties  related  to  plant-available  soil-water  (effective  soil  depth,  sand  and  clay 
contents).  The  general  hypothesis  is  that  the  spatial  pattern  of  soil  properties  related  to 
plant-available  soil-water  and  crop  growth,  yield  and  yield  components  of  corn  grown 
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under  non-irrigated  conditions  are  related  to  certain  topographic  terrain  attributes.  In 
addition,  we  aimed  to  understand  how  topography  affects  crop  growth  and  yield  for  this 
particular  field.  We  hypothesized  that  the  effects  can  be  through  the  short-term  effect  on 
soil-water  content  (lateral  water  movement)  or  through  its  long-term  (geological  time) 
effects  on  site-specific  soil  properties,  that  in  turn,  determines  differences  in  plant- 
available  soil-water  patterns.  This  quest  is  considered  relevant  because  of  its  implication 
in  terms  of  modeling  efforts.  The  present  results  constitute  part  of  an  analysis  of  a 
broader  effort  to  understand,  describe  and  predict  the  spatial,  within-field,  corn  yield 
patterns. 


Materials  and  Methods 

Field  Experiment 

This  study  was  based  on  the  field  experiment  fully  described  in  chapter  2.  The 
experiment  was  carried  out  in  a corn  farm  field  located  in  Michigan,  USA  during  the  crop 
seasons  of  1997  and  1998.  It  consisted  of  observing  variables  of  interest  in  a uniformly 
managed  field.  The  sites  for  observing  soil  and  crop  variables  (total  of  43)  were  spatially 
distributed  in  a grid  of  approximately  30.4-m  (Figure  1). 

Effective  soil  depth  (simply  designated  as  soil  depth)  and  particle  size  analysis 
(carried  out  for  each  15 -cm  increment  layer  up  to  the  constraining  clay  layer)  were 
determined  prior  to  planting  in  1997.  Plant  weight  and  leaf  area  index  were  determined 
at  anthesis  (DOY  21 1 in  1997  and  207  in  1998).  Grain  yield  and  its  components,  namely 
plant  population  (PP),  grain  number  per  ear  (G#EAR),  grain  number  per  unit  area 
(G#AREA)  and  unit  grain  weight  (UGW)  were  determined  at  maturity  (DOY  289  in 
1997  and  288  in  1998).  Rainfall  data  were  recorded  daily. 
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Topographic  Attributes 

We  used  three  different  types  of  topographic  terrain  attributes  as  defined  in  Moore 
et  al.  (1993):  topographic-based  (altitude,  slope,  plan  curvature,  profile  curvature,  local 
relief  and  surface  roughness),  hydrology-based  (flow  direction,  specific  catchment  area, 
distance  to  network,  downstream  and  upstream  flow  length)  and  secondary  attributes 
(wetness  index  and  stream  power  index).  Elevation  data  consisted  of  a 1 m-resolution 
digital  elevation  model  (DEM)  of  the  field  obtained  by  F.J.  Pierce  and  T.G.  Mueller 
(personal  communication,  1997)  using  a precise  deferential  GPS.  The  raw  DEM  was 
checked  for  the  existence  of  sinks  and  none  were  found.  The  definitions  of  the 
topographic  terrain  attributes  used  in  the  study  are  presented  in  Table  5. 

FDIR  was  calculated  assuming  a unidirectional  model.  If  two  directions  were 
designated  as  steepest  descent,  the  final  flow  direction  was  calculated  as  the  intermediate 
between  the  two.  Physically,  topographic  curvature  indicates  if  the  slope  is  convex  or 
concave.  PRC  and  PC  are  two  complementary  ways  of  defining  curvature.  The  profile 
curvature  affects  the  acceleration  (concave  slope  and  curvature  smaller  than  zero)  and 
deceleration  of  flow  (convex  slope  and  curvature  greater  than  zero).  The  plan  curvature 
influences  convergence  (concave  slope  and  curvature  smaller  than  zero)  and  divergence 
(convex  slope  and  curvature  greater  than  zero)  of  flow.  The  catchment  area  of  1000  m2 
used  to  define  a stream  network  and,  consequently,  DTN  was  considered  to  be 
appropriate  for  the  general  topography  characteristics  of  the  field. 
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Table  5.  Definitions  of  the  topographic  terrain  attributes  used  in  the  study  to  characterize 
the  field  topography. 


Attribute 

Acronym 

Definition 

Altitude 

ELEV 

Elevation  (m) 

Slope 

SLP 

Maximum  relative  rate  of  change  in  elevation  from 
each  cell  to  its  eight  neighbors  (%) 

Plan  Curvature 

PLC 

Rate  of  change  of  slope  across  the  flow  path 

Profile  Curvature 

PRC 

Rate  of  change  of  slope  along  the  flow  path 

Local  relief 

LR 

Range  of  elevation  variation  within  the  eight  neighbor 
cells  (m) 

Surface 

Roughness 

SR 

Variance  of  elevation  of  the  eight  neighbor  cells  (m2) 

Flow  Direction 

FDIR 

Direction  of  steepest  descent  from  each  cell 

Specific 

Catchment  Area 

SCA 

Upslope  (catchment)  area  per  unit  cell  width  (m2/m) 

Distance  to 
network 

DTN 

Horizontal  distance  to  a flow  network  defined  by  the 
cells  that  had  an  catchment  area  of  at  least  1 000  m2 
(m) 

Downstream  Flow 
Length 

DFL 

Flow  path  distance  to  the  corresponding  flow  outlet 
(m) 

Upstream  Flow 
Length 

UFL 

Longest  flow  path  distance  contributing  to  the 
specified  cell  (m) 

Wetness  Index 

WI 

Logarithm  of  ratio  between  SCA  and  SLP 

Stream  Power 
Index 

SPI 

Logarithm  of  product  between  SCA  and  SLP 
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Statistical  Analysis 

Maps  showing  the  spatial  pattern  of  variability  of  crop  and  soil  variables  were 
developed.  Semivariance  analysis  and  kriging  were  used  for  interpolation  proposes 
(Issacks  and  Srivastava,  1989).  The  variograms  were  derived  by  plotting  the 
semivariance  of  the  data  as  function  of  vector  distance  before  being  modeled.  The 
software  Variowin  (Pannatier,  1996)  and  Surfer  (Golden  Software,  CO)  were  used  for  the 
analyses. 

The  relation  between  independent  (topographic  attributes)  and  dependent 
variables  (soil  and  crop  variables)  was  studied  using  correlation  analysis  for  individual 
relations  and  multiple  linear  regression  analysis  for  combined  effects.  Before  proceeding 
with  correlation  analysis  between  independent  and  dependent  variables,  the  correlation 
coefficients  between  topographic  attributes  were  calculated  to  eliminate  redundant 
variables. 

The  rationale  for  understanding  how  topography  affects  crop  growth  and  yield  for 
this  particular  field  was  based  on  a correlation  analysis.  If  topographic  attributes  are 
better  correlated  to  soil  variables  than  to  crop  variables  and  crop  variables  are  better 
correlated  to  soil  variables  than  to  topographic  attributes,  then  there  is  evidence  that 
topography  impact  on  crop  variables  is  through  its  long  term  (geological  time)  effects  on 
site-specific  soil  properties,  which  in  turn,  determines  deferential  plant-available  soil- 
water  patterns.  On  the  other  hand,  if  topographic  attributes  are  better  correlated  to  crop 
variables  than  to  soil  variables  and  crop  variables  are  better  correlated  to  topographic 
attributes  than  to  soil  variables,  then  there  is  evidence  that  topography  impact  on  crop 
variables  is  through  short-term  effect  on  soil-water  content  (lateral  water  movement). 
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Results  and  Discussion 

Weather 

Rainfall  in  1997  was  lower  than  in  1998  during  the  crop  season  and  whole  year 
(Figure  2).  The  final  values  for  1997  and  1998  were  849.1  and  529.1  mm  for  the  year 
and  225.9  mm  and  436.1  mm  (difference  was  210.2  mm)  for  the  crop  season.  This 
difference  made  these  two  years  satisfactory  to  study  the  interaction  between  the  effects 
of  seasonal  pattern  of  rainfall  and  the  effects  of  topographic  attributes  on  crop  growth  and 
yield.  From  a weather  series  (1979-1988)  of  the  same  location,  the  average  seasonal 
rainfall  amount  was  459  mm  with  minimum  and  maximum  values  of  313  and  701  mm. 
Considering  these  statistics,  1998  can  be  considered  an  average  year  in  terms  of  seasonal 
rainfall  whereas  1997  was  a dry  one  (the  driest  in  the  series). 

The  longest  period  of  rainfall  absence  (50  days)  in  1997  followed  planting.  This 
dry  period  was  accompanied  by  low  air  temperature  causing  emergence  to  be  delayed 
(Figure  3).  In  1997  emergence  occurred  30  days  after  planting  versus  only  8 days  in 
1998.  Consequently,  PP  at  maturity  was  more  variable  in  1997  than  in  1998. 

Independent  Variables:  Topographic  Attributes 

Figure  16  through  Figure  25  display  the  map  of  each  of  the  topographic  attributes 
used  in  the  study.  Figure  26  shows  the  stream  network  used  to  calculate  DTN,  DFL  and 
UFL  topographic  attributes  as  well  as  the  field  drainage  basins.  The  elevation  relief  of 
the  field  was  2.65  m.  This  gradient  was  well  distributed  in  the  area  with  an  average  slope 
of  1 .3%.  Only  three  small  areas  (about  0.7%  of  the  field  area)  had  slopes  higher  than 
4.5%.  These  three  areas  corresponded  to  field  outlet  locations.  Only  one  of  them  (in  the 
middle  of  the  field)  belonged  to  a major  drainage  basin  outlet  with  1 .25  ha  contributing  to 
its  flow.  Most  of  the  field  (64.2%)  had  low  slopes  (smaller  than  1 .5%). 
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Figure  16.  Map  of  the  farm  field  used  in  the  experiment  (total  area  of  3.6  ha)  displaying 
the  elevation  (ELEV,  m)  at  each  location  (1-m  resolution).  Lines  represent  contours  of 
designated  elevation. 


Figure  1 7.  Map  of  the  farm  field  used  in  the  experiment  (total  area  of  3.6  ha)  displaying 
the  slope  (SLP,  %)  at  each  location  (1-m  resolution).  Zones  of  higher  slope  correspond  to 
field  outlets  for  water  flow. 
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Figure  1 8.  Map  of  the  farm  field  used  in  the  experiment  (total  area  of  3.6  ha)  displaying 
the  plan  curvature  (PLC)  at  each  location  (1-m  resolution).  See  Table  5 for  attribute 
definition. 
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Figure  19.  Map  of  the  farm  field  used  in  the  experiment  (total  area  of  3.6  ha)  displaying 
the  profile  curvature  (PRC)  at  each  location  (1-m  resolution).  See  Table  5 for  attribute 
definition. 
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Figure  20.  Map  of  the  farm  field  used  in  the  experiment  (total  of  3.6  ha)  displaying  the 
specific  catchment  area  (SCA,  xlOOO  m2  m'1)  at  each  location  (1-m  resolution).  See  Table 
5 for  attribute  definition. 


Figure  2 1 . Map  of  the  farm  field  used  in  the  experiment  (total  area  of  3.6  ha)  displaying 
the  distance  to  network  (DTN,  m)  at  each  location  (1-m  resolution).  See  Table  5 for 
attribute  definition. 
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Figure  22.  Map  of  the  farm  field  used  in  the  experiment  (total  area  of  3.6  ha)  displaying 
the  downstream  flow  length  (DFL,  m)  at  each  location  (1-m  resolution).  See  Table  5 for 
attribute  definition.  Also  shown  are  the  drainage  basins. 


Figure  23.  Map  of  the  farm  field  used  in  the  experiment  (total  area  of  3.6  ha)  displaying 
the  upstream  flow  length  (UFL,  m)  at  each  location  (1-m  resolution).  See  Table  5 for 
attribute  definition. 
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Figure  24.  Map  of  the  farm  field  used  in  the  experiment  (total  area  of  3.6  ha)  displaying 
the  wetness  index  (WI)  at  each  location  (1-m  resolution).  See  Table  5 for  attribute 
definition. 
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Figure  25.  Map  of  the  farm  field  used  in  the  experiment  (total  area  of  3.6  ha)  displaying 
the  stream  power  index  (SPI)  at  each  location  (1-m  resolution).  See  Table  5 for  attribute 
definition. 
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Figure  26.  Map  of  the  farm  field  used  in  the  experiment  displaying  the  stream  network 
defined  by  the  cells  with  a catchment  area  greater  than  1000  m2  (1-m  resolution)  and  the 
drainage  basins  with  areas  larger  than  600  m2.  The  darker  color  in  the  streamline  signifies 
a greater  catchment  area  of  those  cells. 

The  field  had  four  drainage  basins  larger  than  600  m2.  One  was  already  referred 
to  with  an  area  of  1 .25  ha  and  the  other  three  had  approximately  half  of  that  area  (Figure 
26).  Due  to  the  predominantly  low  slopes,  SCA  dominates  calculation  of  both  WI  and 
SPI.  The  higher  values  of  both  indices  appear  on  the  flow  network  defined  from  SCA. 

In  other  words,  higher  upslope  contributing  area  causes  higher  wetness  index  and  stream 
power  index,  independently  of  high  or  low  predominant  slope,  respectively,  for  the  range 
of  slopes  present  in  this  field. 

Curvature  maps  show  the  distribution  of  concave  and  convex  areas  according  to 
the  flow  direction  and  perpendicular  plans.  According  to  PC,  the  zones  of  flow 
convergence  and  divergence  seem  to  equally  dominate  the  landscape.  On  the  other  hand, 
according  to  the  PRC,  zones  of  accelerated  flow  dominate  the  flow  deceleration  zones. 
DFL  resembled  the  DTN  map,  but  since  it  considers  flow  length  and  not  just  horizontal 
distance,  the  shapes  of  the  iso-lines  match  the  drainage  basins  map. 
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The  correlation  coefficients  between  topographic  attributes  are  presented  in  Table 
6.  SR,  LR  and  SLP  were  considered  redundant  since  the  correlation  coefficients  between 
them  were  0.99.  Of  these  attributes,  only  SLP  was  kept  for  the  analysis  because  it  is 
easier  to  interpret  than  SR  and  LR.  Another  highly  correlated  group  of  attributes  was 
SCA,  UFL,  SPI  and  WI.  This  was  expected  since  both  secondary  attributes  (WI  and  SPI) 
were  calculated  based  upon  SCA  and  SLP.  The  high  correlation  coefficients  between 
SCA,  SPI  and  WI  show  that  SLP  is  not  as  relevant  as  SCA  in  characterizing  this 
particular  landscape.  The  correlation  between  SCA  and  UFL  is  obvious:  lengthier 
upslope  flow  corresponds  to  larger  upslope  flow  contributing  area.  The  significant 
correlation  between  the  different  curvatures  (PRC  and  PC)  and  WI  results  from  the  fact 
that  convergent  decelerated  water  flows  tend  to  proportion  more  infiltration  and  therefore 
higher  WI. 

Dependent  Variables:  Soil  Variables 

Figure  4 shows  the  map  of  sand  and  clay  and  soil  depth  is  shown  in  Figure  5.  The 
coefficients  of  determination  between  the  topographic  attributes  and  soil  variables  are 
shown  in  Table  7.  Field  areas  with  higher  elevation,  farther  from  the  stream  network, 
with  lengthier  downstream  flow  and  smaller  wetness  index  had  higher  (lower)  sand  (clay) 
content  and  vice-versa.  Plan  curvature  had  a significant  correlation  coefficient  for  both 
sand  and  clay  contents  indicating  that  zones  of  convergent  flow  had  higher  clay  and 
lower  sand  contents.  However,  most  of  the  sites  had  plan  curvatures  between  ± 0.2. 

Only  seven  sites  had  values  outside  that  interval  and  those  were  responsible  for  the 
significant  correlation  coefficients. 


49 


Table  6.  Correlation  coefficients  between  topographic  attributes.  Topographic  attribute 
acronyms  are  defined  in  Table  5.  Non  significant  relationships  are  marked  with  n.s.  The 
significant  correlations  are  classified  according  to  the  level  of  significance  (*  p<0.05,  ** 

p<0.01,  ***  p<0.001). 
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Table  7.  Coefficients  of  determination  between  independent  and  dependent  variables. 
Topographic  attribute  acronyms  are  defined  in  Table  5.  Non  significant  relationships  are 
marked  with  n.s.  The  significant  correlations  are  classified  according  to  the  level  of 
significance  (*  p<0.05,  **  p<0.01,  ***  p<0.001). 


Sand 

Clay 

Soil  depth 

ELEV 

(+)0.30** 

(-)0.22** 

(-)0.29** 

SLP 

n.s. 

n.s. 

(-)0. 14* 

SCA 

n.s. 

n.s. 

n.s. 

DTN 

(+)0.41  *** 

(-)0.40*** 

(-)0.23** 

WI 

(-)0.28** 

(+)0.26** 

n.s. 

SPI 

n.s. 

n.s. 

n.s. 

PC 

(+)0.27** 

(-)0.23* 

(-)0.19** 

PRC 

n.s. 

n.s. 

n.s. 

DFL 

(+)0.37*** 

(-)0.40*** 

(-)0 .62*** 

UFL 

n.s. 

n.s. 

n.s. 

Field  areas  with  lower  elevation,  shorter  downstream  flow,  lower  slope  and  closer 
to  the  flow  network  had  higher  effective  soil  depths.  The  strongest  relation  in  the  study 
was  between  effective  soil  depth  and  downstream  flow  length. 

Dependent  Variables:  Crop  Variables 

Figure  9 through  Figure  13  show  the  maps  of  the  crop  variables  used  in  the  study. 
Table  8 and  Table  9 display  the  coefficients  of  determination  between  topographic 
attributes  and  crop  variables  for  1997  and  1998,  respectively. 
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Table  8.  Coefficients  of  determination  between  independent  and  dependent  variables  for 
1997  corn  growing  season.  Topographic  attribute  acronyms  are  defined  in  Table  5.  Non 
significant  relationships  are  marked  with  n.s.  The  significant  correlations  are  classified 
according  to  the  level  of  significance  (*  p<0.05,  **  p<0.01,  ***  p<0.001). 


Grain 

Yield 

Grain 

number  per 
ear 

Unit  grain 
weight 

Grain 
number 
per  unit 
area 

Plant 

Population 

Plant 
weight  at 
anthesis 
per  unit 
area 

Leaf  area 
index  at 
anthesis 

ELEV 

(-)0.30** 

(+)0.29** 

(-)0.27** 

n.s. 

(-)0.26** 

(-)0.20** 

(-)0.37** 

SLP 

n.s. 

(+)0.17* 

(-)0. 14* 

n.s. 

n.s. 

n.s. 

n.s. 

SCA 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

DTN 

(-)0.34** 

(+)0.17* 

(-)0.28** 

n.s. 

(-)0.23** 

(-)0.28** 

(-)0.49*** 

WI 

(+)0.12* 

n.s. 

(+)0.13* 

(+)0.13* 

n.s. 

(+)0.14* 

(+)0. 15* 

SPI 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

PC 

(-)0.20** 

(+)0.17** 

(-)0.19** 

n.s. 

n.s. 

(-)0.19** 

(-)0. 1 8** 

PRC 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

DFL 

(-)0.57*** 

(+)0.28** 

(-)0.45*** 

n.s. 

(-)0.21** 

(-)0.29** 

(-)0.41*** 

UFL 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 
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Table  9.  Coefficients  of  determination  between  independent  and  dependent  variables  for 
1998  com  growing  season.  Topographic  attribute  acronyms  are  defined  in  Table  5.  Non 
significant  relationships  are  marked  with  n.s.  The  significant  correlations  are  classified 
according  to  the  level  of  significance  (*  p<0.05,  **  p<0.01,  ***  p<0.001). 


Grain  Yield 

Grain 
number 
per  ear 

Unit  grain 
weight 

Grain 
number 
per  unit 
area 

Plant 

Population 

Plant 
weight  at 
anthesis 
per  unit 
area 

Leaf  area 
index  at 
anthesis 

ELEV 

(-)0.24** 

(+)0.14* 

(-)0. 14* 

n.s. 

n.s. 

(-)0. 14* 

n.s. 

SLP 

n.s. 

n.s. 

(-)0. 12* 

n.s. 

n.s. 

n.s. 

n.s. 

SCA 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

DTN 

(_)0  jo*** 

n.s. 

(-)0. 19** 

n.s. 

n.s. 

n.s. 

(-)0.10* 

WI 

(+)0.10* 

n.s. 

(+)0.10* 

n.s. 

n.s. 

n.s. 

n.s. 

SPI 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

PC 

(-)0.20** 

n.s. 

(-)0.19** 

n.s. 

n.s. 

n.s. 

n.s. 

PRC 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

DFL 

(-)0.29*** 

n.s. 

(_)0.34*** 

n.s. 

n.s. 

n.s. 

n.s. 

UFL 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 

n.s. 
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The  number  of  significant  relationships  between  topographic  attributes  and  crop 
variables  in  1997  was  higher  and  they  were  larger  in  value  when  compared  to  1998.  This 
fact  makes  sense  considering  that  1997  was  a drier  year  than  1998,  and  therefore  PASW 
was  more  limiting  to  crop  growth  in  1997. 

The  topographic  attributes  that  were  more  significantly  correlated  with  crop 
variables  were  the  same  for  both  years:  ELEV,  DTN,  WI,  PC  and  DFL.  Furthermore, 
these  are  the  topographic  attributes  that  showed  the  highest  correlation  with  soil 
properties.  DFL  was  the  topographic  attribute  that  had  the  highest  correlation  with  yield 
and  most  yield  components.  This  attribute  was  also  more  highly  correlated  to  soil 
variables,  especially  soil  depth.  Similarly,  DTN  was  the  second  most  correlated 
topographic  attribute  with  yield  as  well  as  with  soil  variables. 

In  terms  of  yield  components,  UGW  showed  the  highest  correlation  to 
topographic  attributes  in  both  years.  In  1998,  UGW  was  the  only  yield  component  that 
showed  any  correlation  to  topographic  attributes.  This  was  probably  because  this  yield 
component  was  formed  in  the  later  part  of  the  growing  season  when  stored  soil-water 
played  an  important  role  in  supplying  water  to  the  crop.  In  1997,  G#EAR  also  showed 
some  degree  of  correlation  with  the  topographic  attributes.  In  neither  of  the  years  was 
G#AREA  correlated  with  topographic  attributes.  In  1997,  PP  showed  some  degree  of 
correlation  with  DFL,  DTN  and  ELEV.  The  crop  growth  variables  were  also  more  highly 
correlated  to  topographic  attributes  in  1 997. 

These  observations  are  consistent  with  the  relationships  described  previously 
concerning  soil  properties  and  yield  and  yield  components  (chapter  2).  Effective  soil 
depth  and  clay  and  sand  contents  affected  yield  formation.  On  the  other  hand,  they  were 


54 


affected  by  topography  (especially  related  to  altitude,  distance  to  network  and 
downstream  flow  length).  These  relationships  were  affected  by  the  interactions  with 
weather,  in  particular,  rainfall.  In  the  dry  year  (1997),  the  correlations  between 
topographic,  soil  and  crop  variables  were  stronger  than  in  the  normal  rainfall  year  (1998). 

The  results  seem  to  show  that  for  dry  and  normal  rainfall  years,  the  effect  of 
topography  on  the  spatial  variability  of  com  growth  and  yield,  in  this  particular 
landscape,  is  relatively  more  important  through  the  long  term  effect  (geological  time)  on 
soil  properties  (particularly  effective  soil  depth)  than  the  short  term  effect  on  seasonal 
plant-available  soil-water.  In  wetter  years,  the  relative  importance  of  the  two  effects 
might  change  since  the  magnitude  of  the  seasonal  lateral  surface  and  sub-surface  soil- 
water  movement  changes. 


Conclusions 

Landscape  topography  impacts  soil  properties  and  corn  growth  and  yield  in  non- 
irrigated  conditions.  For  the  particular  field  studied,  effective  soil  depth  was  more 
strongly  related  to  topography  than  sand  and  clay  content.  The  topographic  attributes 
more  related  to  these  soil  properties  were  downstream  flow  length  and  distance  to  flow 
network.  Corn  growth  and  yield  were  consistently  related  to  topography  in  both  a dry 
and  a normal  rainfall  year.  These  relations  were  stronger  for  the  dry  year,  especially  for 
downstream  flow  length  and  distance  to  network.  The  yield  component  that  more 
strongly  correlated  to  topography  was  unit  grain  weight.  The  results  indicate  that  the 
effect  of  topography  on  corn  growth  and  yield  for  dry  and  normal  rainfall  years  was 
mainly  explained  through  its  effect  on  soil  properties. 


CHAPTER  4 

PREDICTING  THE  SPATIAL  PATTERN  OF  CORN  YIELD  USING  A NEURAL 

NETWORK  MODEL 

Introduction 

Once  the  key  factors  responsible  for  spatial  yield  variability  are  identified, 
predicting  that  variability  becomes  the  next  challenge.  Being  able  to  predict  the  spatial 
pattern  of  grain  yield  is  the  ultimate  proof  that  the  most  important  limiting  factors  have 
been  identified.  If  we  can  predict  yield  variability,  it  means  that  we  have  an  algorithm 
that  relates  environmental  variables  and  yield.  Such  an  algorithm  is  of  crucial  importance 
to  site-specific  management  (SSM)  because  it  allows  for  economically  and  agronomically 
feasible  recommendations  to  be  made  (Sudduth  et  al.,  1996). 

Most  of  the  algorithms  used  to  estimate  yield  from  environmental  factors  can  be 
divided  into  empirical  and  process-based  tools.  The  basic  difference  between  the  two 
approaches  is  that  process-based  algorithms  attempt  to  predict  final  grain  yield  by 
describing  the  most  important  physical/chemical/biological  processes  responsible  for 
yield  formation.  An  example  of  this  approach  is  a crop  simulation  model.  Empirical 
methods,  usually  static,  simply  describe  non-physical  relations  between  environment  and 
yield.  Examples  of  this  approach  are  statistical  methods  such  as  multiple  linear 
regression,  tree  regression  analysis  and  neural  network  models.  Both  approaches  have 
advantages  and  disadvantages. 

The  method  that  is  of  interest  in  this  paper  is  neural  network  models. 
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Neural  Network  Models 

Neural  network  models  are  especially  well  suited  for  complex  systems  where  the 
underlying  processes  are  not  well  understood  or  are  difficult  to  translate  into 
mathematical  equations  (Hashimoto,  1997;  Schultz  and  Wieland,  1997).  In  such  cases, 
the  researcher  is  exposed  to  a myriad  of  environmental  variables  that  are  supposed  to 
have  an  effect  on  the  variable(s)  of  interest.  This  is  the  particular  case  of  SSM  where  in 
most  projects,  the  researcher  has  collected  many  soil  parameters  and  weather  and  crop 
variables  and  has  the  task  of  understanding  grain  yield  variability.  In  some  cases,  linear 
regression  can  be  applied  if  the  relations  are  not  too  complex,  meaning  that  they  are  not 
non-linear  and  correlated.  Since  this  is  not  the  case  for  most  situations,  and  given  the  fact 
that  most  crop  simulation  models  have  yet  to  incorporate  some  of  the  important  yield 
limiting  variables,  neural  network  models  appear  to  be  a promising  technique. 

A good  introduction  into  neural  networks  is  given  by  Batchelor  (1998).  Neural 
network  models  have  been  developed  as  a branch  of  artificial  intelligence.  They  were 
designed  to  mimic  the  human  learning  process  by  repeatedly  exposing  the  neural  network 
with  several  representative  sets  of  input-output  relationships.  After  the  neural  network 
model  is  trained,  i.e.,  after  the  internal  parameters  have  been  optimized  to  minimize  the 
differences  between  neural  network  outputs  and  real  systems  responses,  it  can  be  used  to 
predict  the  system  response  in  a totally  new  set  of  input  conditions.  Of  course  its  ability 
to  do  this  depends  on  how  well  it  has  been  trained.  A well-trained  neural  network  model 
is  one  that  has  been  trained  with  a wide  range  of  input-output  relationships  (called 
patterns)  and,  therefore,  is  able  to  be  extrapolated  to  new  patterns. 
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Although  a trained  neural  network  captures  the  essence  of  the  meaningful  system 
relationships,  it  cannot  be  translated  into  'real  system  knowledge'  (systems  governing 
equations),  i.e.,  it  uses  a totally  empirical  approach.  This  constitutes  one  of  the 
disadvantages  of  these  models.  But  at  the  same  time  this  can  be  viewed  as  an  advantage 
if  one  is  only  interested  in  prediction.  Neural  networks  can  also  be  used  to  reveal 
meaningful  relationships  as  part  of  an  exploratory  analysis  exercise. 

The  computational  implementation  of  a neural  network  is  done  in  complete 
analogy  to  biological  neural  systems.  It  contains  three  distinct  types  of  layers:  input, 
hidden  and  output.  The  number  of  hidden  layers  can  vary.  On  the  other  hand  there  is 
only  one  input  and  output  layer.  Each  hidden  layer  has  its  own  number  of  nodes  whereas 
the  numbers  of  nodes  in  the  input  and  output  layers  are  determined  by  the  number  of 
desired  inputs  and  outputs.  Each  connection  between  nodes  (all  nodes  in  one  layer  are 
connected  to  the  nodes  in  the  following  layer)  has  a weight  that  is  multiplied  by  the  input 
signal.  Each  node  has  its  transfer  function  that  is  used  to  transform  the  sum  of  all  input 
signals.  The  weights  of  the  different  connections  in  a network  are  optimized  during  the 
learning  process. 

There  are  several  types  of  neural  networks  differentiated  by  the  number  of  nodes 
in  the  hidden  layer(s),  the  type  of  transfer  function,  the  method  of  connecting  the  nodes 
together  and  the  methods  of  computing  the  weights. 

The  most  widely  used  connection  method  is  a feedforward  architecture  because  it 
is  simple  and  able  to  generalize  relationships  very  well.  In  this  type  of  network,  data  are 
processed  from  left  to  right  through  the  layers  before  producing  the  output  values.  Data 
processing  consists  of  multiplying  all  input  signals  from  the  previous  layer  by  each 
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connection  weight  and  then  transforming  the  multiplication's  sum  using  a transfer 
function. 

The  transfer  function  is  usually  a sigmoid  curve  that  transforms  all  input  signals 
into  a numerical  range  of  0 to  1 . The  number  of  hidden  layers  and  their  nodes  depend  on 
the  complexity  of  the  problem  to  solve.  A network  with  three  hidden  layers  is  known  to 
capture  all  possible  levels  of  complexity,  but  in  most  cases  one  layer  is  enough  to  handle 
the  existing  complexity.  The  number  of  nodes  in  the  hidden  layer  is  usually  determined 
by  a sensitivity  analysis. 

The  method  to  update  the  connection  weights  is  usually  the  back-propagation 
algorithm.  In  its  most  simple  case,  the  update  is  done  by  multiplying  the  error  (sum  of 
squares  of  individual  errors)  by  a user-defined  parameter  designated  as  learning  rate  (0  to 
1 parameter).  The  weights  are  updated  from  right  to  left.  The  larger  the  learning  rate  the 
faster  the  network  will  converge  to  a solution.  In  complex  problems  (non-smooth  error 
surfaces)  this  may  result  in  local  instead  of  global  optimal  solutions,  and  a sensitivity 
analysis  is  needed  to  determine  the  best  learning  rate  for  specific  problems. 

Objective 

The  objective  of  this  paper  is  to  test  the  use  of  a neural  network  model  to  predict 
the  spatial  pattern  of  corn  yield.  An  additional  objective  is  to  study  the  importance  of  the 
different  spatial  variables  in  network  prediction  error  by  testing  several  combinations  of 
inputs.  Three  neural  network  models  were  trained  and  tested:  (1)  using  only  topographic 
variables,  (2)  using  only  agronomic  variables  and  (3)  using  both  types  of  variables. 
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Materials  and  Methods 

Neural  Network  Design.  Training  and  Testing 

Three  types  of  neural  networks  were  implemented,  trained  and  tested.  The  first 
network  (designated  as  topographic  attributes)  used  the  calculated  topographic  attributes 
and  total  season  rainfall  as  inputs.  The  second  network  (designated  as  agronomic 
variables)  used  measured  or  estimated  agronomic  variables  and  total  seasonal  rainfall  as 
inputs.  Finally,  the  third  network  (designated  as  mixed  variables)  used  all  of  the  above 
variables  as  inputs. 

All  three  neural  networks  had  a similar  structure:  feedforward  architecture, 
backpropagation  weight  adjustment  algorithm,  one  hidden  layer  with  sigmoid  transfer 
function  and  linear  transfer  functions  for  input  and  output  layers.  The  inputs  were 
standardized  (to  mean  and  standard  deviation)  before  being  processed  in  the  network. 
The  major  difference  between  the  networks  was  the  number  of  nodes  in  the  hidden  layer. 
For  each  network  a sensitivity  analysis  based  on  the  final  predictive  ability  was  carried 
out  to  find  the  optimal  number  of  nodes. 

The  training  and  validation  sets  were  obtained  from  a two-year  experiment  in  one 
field.  For  this  reason,  the  data  set  had  to  be  split.  In  order  to  provide  four  data  sets  in  a 
spatial  continuum,  the  field  was  divided  into  two  sections  along  the  transversal  direction. 
In  this  way,  two  pairs  of  training  and  testing  sets  were  obtained.  The  neural  networks 
were  trained  with  two-years  of  data  for  the  east  section  of  the  field  and  tested  on  the  west 
section  and  vice-versa.  Since  rainfall  was  obviously  a major  variable  needed  to  explain 
final  grain  yield,  the  whole  field  data  set  for  1997  and  1998  could  not  be  used  as  training 
and  testing  sets  independently. 
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The  training  and  testing  data  sets  were  evaluated  through  the  root  mean  square 
error  of  fit  (RMSEF)  and  of  prediction  (RMSEP)  among  sites,  respectively.  A field-wide 
evaluation  was  also  carried  out  through  visual  inspection  of  the  spatial  pattern  of  yield 
variation  obtained  through  kriging.  In  addition,  the  spatial  patterns  of  the  differences 
between  predicted  and  observed  yield  are  presented. 

Agronomic  Variables 

This  study  was  based  on  the  field  experiment  fully  described  in  chapter  2.  The 
experiment  was  carried  out  in  a corn  farm  field  located  in  Michigan,  USA  during  the  crop 
seasons  of  1997  and  1998.  It  consisted  of  observing  variables  of  interest  in  a uniformly 
managed  field.  The  sites  for  observing  soil  and  crop  variables  (total  of  43)  were  spatially 
distributed  in  a grid  of  approximately  30.4-m  (Figure  1). 

Effective  soil  depth  (simply  designated  as  soil  depth)  and  particle  size  analysis 
(carried  out  for  each  15-cm  increment  layer  up  to  the  constraining  clay  layer)  were 
determined  prior  to  planting  in  1997.  Grain  yield  and  plant  population  (PP)  were 
determined  at  maturity  (DOY  289  in  1997  and  288  in  1998).  Rainfall  data  was  recorded 
daily. 

Season  Rainfall 

Rainfall  data  were  recorded  daily  by  a tipping  bucket  precipitation  sensor  (1-mm 
accuracy)  that  was  located  200  meters  from  the  border  of  the  field.  The  reset  time  for 
daily  amounts  was  done  at  midnight.  Rainfall  in  1997  was  lower  than  in  1998  during  the 
crop  season  and  whole  year  (Figure  2).  The  final  values  for  1997  and  1998  were  849.1 
and  529.1  mm  for  the  year  and  225.9  mm  and  436.1  mm  (difference  was  210.2  mm)  for 
the  crop  season.  This  difference  made  these  two  years  satisfactory  to  study  the 
interaction  between  the  effects  of  seasonal  pattern  of  rainfall  and  the  effects  of 
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topographic  attributes  on  crop  growth  and  yield.  From  a weather  series  (1979-1988)  of 
the  same  location,  the  average  seasonal  rainfall  amount  was  459  mm  with  minimum  and 
maximum  values  of  313  and  701  mm.  Considering  these  statistics,  1998  can  be 
considered  an  average  year  in  terms  of  seasonal  rainfall  whereas  1997  was  a dry  one  (the 
driest  in  the  series). 

Topographic  Attributes 

We  used  the  topographic  terrain  attributes  as  defined  in  Moore  et  al.  (1993)  and 
fully  described  in  chapter  3.  The  definitions  of  the  topographic  terrain  attributes  used  in 
the  study  are  presented  in  Table  5.  Figure  16  through  Figure  25  display  the  map  of  each 
of  the  topographic  attributes  used  in  the  study.  FDIR  was  calculated  assuming  a 
unidirectional  model.  If  two  directions  were  designated  as  steepest  descent,  the  final 
flow  direction  was  calculated  as  the  intermediate  between  the  two.  Physically, 
topographic  curvature  indicates  if  the  slope  is  convex  or  concave.  PRC  and  PLC  are  two 
complementary  ways  of  defining  curvature.  The  profile  curvature  affects  the  acceleration 
(concave  slope  and  curvature  smaller  than  zero)  and  deceleration  of  flow  (convex  slope 
and  curvature  greater  than  zero).  The  plan  curvature  influences  convergence  (concave 
slope  and  curvature  smaller  than  zero)  and  divergence  (convex  slope  and  curvature 
greater  than  zero)  of  flow.  The  catchment  area  of  1 000  m2  used  to  define  a stream 
network  and,  consequently,  DTN  was  considered  to  be  appropriate  for  the  general 
topography  characteristics  of  the  field. 
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Results  and  Discussion 

Table  10  shows  the  RMSEF  for  training  sets  and  RMSEP  for  testing  sets  of  all  the 
neural  networks  models.  Figure  27  through  Figure  32  display  the  respective  predicted 
and  fitted  grain  yield  versus  observed  grain  yield  for  both  training  and  testing  sets. 

As  expected,  the  RMSEF  for  all  training  sets  were  smaller  than  the  ones  for 
testing  sets.  The  differences  between  the  RMSEF  and  RMSEP  for  training  and  testing 
sets,  respectively,  ranged  from  100  to  230  kg/ha  and  was  progressively  smaller  on 
average  from  the  topography  variable  (170  kg/ha)  to  the  mixed  variables  neural  network 
(122  kg/ha). 


Table  10.  Root  mean  square  error  (kg/ha)  for  the  different  neural  network  models.  The 
average  relative  error  (%)  is  presented  between  parentheses. 


T opography-based 

Agronomic-based 

Mixed 

Training 

Testing 

Training 

Testing 

Training 

Testing 

set 

set 

set 

set 

set 

set 

East 

1006.3 

1236.2 

869.4 

991.4 

805.2 

904.9 

section  t 

(9.6) 

(11.9) 

(8.4) 

(9.5) 

(7.7) 

(8.7) 

West 

1181.3 

1291.1 

902.4 

1076.1 

787.2 

931.9 

section  J 

(11.4) 

(12.4) 

(8.7) 

(10.4) 

(7.6) 

(9.0) 

Average 

1093.8 

1263.7 

885.9 

1033.8 

796.2 

918.4 

(10.5) 

(12.2) 

(8.6) 

(10.0) 

(7.7) 

(8.9) 

t Training  is  performed  on  east  section  and  testing  on  the  west  section. 


t Training  is  performed  on  west  section  and  testing  on  the  east  section. 


Predicted  Yield  (Kg/ha)  Fitted  Yield  (Kg/ha) 


63 


Figure  27.  Fitted  versus  observed  grain  yield  for  east-section  training  set  (above)  and 
predicted  versus  observed  grain  yield  for  east-section  testing  set  (below)  for  the 
topographic  attributes  neural  network. 


Predicted  Yield  (Kg/ha)  Fjtted  Yje|d  (Kg/ha) 
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Figure  28.  Fitted  versus  observed  grain  yield  for  west-section  training  set  (above)  and 
predicted  versus  observed  grain  yield  for  west-section  testing  set  (below)  for  the 
topographic  attributes  neural  network. 


Predicted  Yield  (Kg/ha)  Fitted  Yield  (Kg/ha) 
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Figure  29.  Fitted  versus  observed  grain  yield  for  east-section  training  set  (above)  and 
predicted  versus  observed  grain  yield  for  east-section  testing  set  for  the  agronomic 
variables  neural  network. 


Predicted  Yield  (Kg/ha)  Fitted  Yield  (Kg/ha) 
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Figure  30.  Fitted  versus  observed  grain  yield  for  west-section  training  set  (above)  and 
predicted  versus  observed  grain  yield  for  west-section  testing  set  (below)  for  the 
agronomic  variables  neural  network. 


Predicted  Yied  (Kg/ha)  Fitted  Yje|d  (Kg/ha) 
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Figure  3 1 . Fitted  versus  observed  grain  yield  for  east-section  training  set  (above)  and 
predicted  versus  observed  grain  yield  for  east-section  testing  set  (below)  for  the  mixed 
variables  neural  network. 


Predicted  Yield  (Kg/ha)  Fitted  Yield  (Kg/ha) 


68 


Figure  32.  Fitted  versus  observed  grain  yield  for  west-section  training  set  (above)  and 
predicted  versus  observed  grain  yield  for  west-section  testing  set  for  the  mixed  variables 
neural  network. 
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The  topographic  attributes  neural  network  had  the  worst  predictive  ability.  The 
RMSEP  averaged  1264  kg/ha  for  this  neural  network  model,  whereas  for  the  agronomic 
variables  network  the  same  average  was  1034  kg/ha,  a decrease  of  230  kg/ha.  The  neural 
network  that  accounted  for  all  variables  in  the  study  resulted  in  a decrease  in  RMSEP  of 
1 15.4  kg/ha  below  that  of  the  agronomic  variables  neural  network.  The  gain  in 
predictability  by  using  all  variables  did  not  increase  linearly.  Instead,  it  was  only  half  of 
the  gain  when  only  agronomic  variables  were  used.  The  possible  reason  for  this  is  the 
fact  that  the  agronomic  variables  are  correlated  with  the  topographic  variables,  especially 
sand  and  clay  contents  and  effective  soil  depth.  Including  all  variables  in  the  neural 
network  model  increased  the  predictability  of  final  grain  yield. 

Figure  33  and  Figure  34  show  the  predicted  yield  maps  for  the  three  neural 
network  models  for  1997  and  1998,  respectively.  Figure  35  and  Figure  36  show  the  maps 
of  the  differences  between  predicted  and  observed  yield  for  the  three  neural  network 
models  for  1997  and  1998,  respectively.  When  compared  to  the  observed  yield  maps 
(Figure  7),  we  can  see  that  the  predicted  yield  maps  captured  the  major  patterns  of  spatial 
variability.  Nevertheless,  gross  errors  are  easily  spotted  in  these  maps.  This  is  especially 
true  for  low  and  high  yield  production  areas  and  for  topographic  variables  and  agronomic 
variables  neural  network  models.  For  the  mixed  variables  neural  network  model,  the 
error  did  not  reach  1500  kg/ha  for  the  entire  field  in  1998,  and  exceeded  1500  kg/ha  for 
only  a small  area  of  the  field  in  1997  (Figure  35  and  Figure  36). 
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Figure  33.  Predicted  yield  maps  (kg/ha)  from  the  topographic  variables  (above), 
agronomic  variables  (middle)  and  mixed  variables  (below)  neural  networks  for  1997. 
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Figure  34.  Predicted  yield  maps  (kg/ha)  from  the  topographic  variables  (above), 
agronomic  variables  (middle)  and  mixed  variables  (below)  neural  networks  for  1998. 
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Figure  35.  Map  of  the  difference  between  predicted  and  observed  grain  yield  (kg/ha)  for 
the  topographic  variables  (above),  agronomic  variables  (middle)  and  mixed  variables 
(below)  neural  networks  for  1 997. 
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Figure  36.  Map  of  the  difference  between  predicted  and  observed  grain  yield  (kg/ha)  for 
the  topographic  variables  (above),  agronomic  variables  (middle)  and  mixed  variables 
(below)  neural  networks  for  1998. 
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Conclusions 

In  theory,  neural  network  models  are  good  candidates  for  relating  environmental 
and  management  variables  to  yield  and,  therefore,  for  predicting  yield.  We  successfully 
tested  the  prediction  ability  of  three  types  of  neural  networks  using  two  years  of  observed 
data.  In  both  years  the  neural  network  that  included  all  variables  in  the  study,  i.e., 
agronomic  variables  (plant  population,  soil  texture  and  soil  depth),  topographic  variables 
and  seasonal  rainfall,  produced  smaller  predictive  errors  (average  of  91 8 kg/ha).  The 
other  two  neural  networks  (based  only  on  agronomic  variables  or  on  topographic 
variables)  always  had  average  errors  greater  than  1 000  kg/ha.  Due  to  existing  correlation 
in  observed  variables,  the  predictive  error  did  not  decrease  linearly  with  the  amount  of 
variables  included  in  the  neural  network  models. 

Neural  network  models  allowed  for  new  sources  of  data  (topography)  to  be  used 
in  predicting  crop  performance.  The  model  that  included  all  soil,  crop  and  topographic 
variables  presented  a reasonable  predictability  of  the  spatial  variability  of  final  grain  yield 
for  practical  applications.  However,  one  serious  limitation  of  this  approach  is  the  limited 
gain  in  meaningful  system  relationships  that  could  be  translated  into  'real  system 
knowledge'  and  used  to  predict  the  spatial  variability  of  final  grain  yield  in  other 
environments  (different  weather,  soils,  fields,  etc).  Nevertheless,  if  the  future  of  the 
application  of  prediction  tools  in  SSM  is  to  involve  specific  field  calibrations,  neural 
network  models  are  promising.  Their  role  can  either  be  as  independent  tools  in  predicting 
grain  yield  from  environmental  variables  or  as  a support  to  other,  more  mechanistic 
simulation  tools,  in  providing  initial  conditions  values  or  site-specific  soil  parameters. 
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CHAPTER  5 

PREDICTING  WITHIN  FIELD  SPATIAL  VARIABILITY  OF  CORN  GROWTH  AND 
YIELD  WITH  A CROP  SIMULATION  MODEL 

Introduction 

Predicting  within-field  spatial  variability  of  grain  yield  should  be  considered  as  a 
major  objective  in  Site-Specific  Management  (SSM).  Having  a tool  that  relates  resource 
availability  to  yield  allows  for  testing  of  management  prescriptions  and  for  correct 
assessment  of  agronomic  and  economic  outcomes  (Sudduth  et  ah,  1996). 

Most  of  the  algorithms  used  to  estimate  yield  from  environmental  factors  can  be 
divided  into  empirical  and  process-based  tools.  The  basic  difference  between  the  two 
approaches  is  that  the  process-based  algorithms  attempt  to  predict  final  grain  yield  by 
describing  the  most  important  physical/chemical/biological  processes  responsible  for 
yield  formation.  An  example  of  this  approach  is  a crop  simulation  model.  Empirical 
methods,  usually  static,  simply  describe  non-physical  relations  between  environment  and 
yield.  Examples  of  this  approach  are  statistical  methods  such  as  multiple  linear 
regression  and  neural  network  models.  Both  approaches  have  advantages  and 
disadvantages. 

The  method  that  is  of  particular  focus  in  this  paper  is  crop  simulation  models. 

Crop  Simulation  Models 

A crop  simulation  model  is  a mathematical  algorithm  that  relates  genetic  and 
environmental  variables  (soil,  weather  and  management)  to  yield  and  yield  components. 
This  tool  provides  a comprehensive  approach  to  predicting  crop  growth  and  development 
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and  their  environmental  determinant  factors  and  interactions.  By  simulating  intermediate 
processes  that  lead  to  yield  formation,  crop  models  are  more  general  and  robust  than 
empirical  procedures  when  used  in  independent  conditions. 

The  use  of  crop  simulation  models  in  site-specific  management  is  not  new. 

Moore  and  Tyndale-Biscoe  (1999)  stressed  their  advantages,  when  compared  to  field 
assessment,  to  address  questions  of  economic  viability  of  different  management 
strategies:  "sufficient  data  under  a variety  of  weather  conditions  (as  well  as  soil  types  and 
genotypes)  are  not  easily  obtained  (through  field  assessment).  In  the  short  term,  field 
assessment  will  certainly  not  answer  reliably  the  questions  about  economic  benefits  that 
the  farmers  are  seeking.  In  the  long  term,  management  practices  are  changing  more 
quickly  than  the  cycle  of  time  required  to  conduct  experiments  that  assess  the  influence 
of  different  weather  conditions"(p,28).  They  used  the  CERES-Wheat  model 
(Hoogenboom  et  ah,  1994)  to  assess  the  strategic  importance  of  soil- water  holding 
capacity  and  spatially  variable  nitrogen  management  as  factors  to  determine  spatial  yield 
variability. 

Sadler  et  al.  (1993)  mentioned  the  use  of  CERES-Maize  (Hoogenboom  et  ah, 

1 994)  and  other  crop  simulation  models  to  understand  yield  constraining  factors  in  South 
Carolina  and  stressed  the  importance  of  good  soil  input  data,  in  particular  those  related  to 
rooting  habits.  They  also  mentioned  the  importance  of  soil-water  balance  simulation  in 
general  and  the  rainfall-runoff  partitioning  and  water  stress  in  particular.  In  1995,  the 
same  researchers  were  still  struggling  with  the  same  difficulties,  hoping  for  newer 
versions  of  the  model  to  address  the  problem  again  (Sadler  et  ah,  1995). 
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Fraisse  et  al.  (1999)  also  used  the  CERES-Maize  simulation  model  with  similar 
objectives  as  Sadler  et  al.  (1995),  but  for  poorly  drained  soils  with  a claypan.  The  version 
ot  the  model  used  in  the  study  (Calmon  et  al.  1999a,  1999b)  included  the  impedance  and 
hospitality  factors  to  better  describe  root  growth  in  such  adverse  conditions.  With  the 
calibrated  model,  the  authors  were  able  to  explain  measured  yield  in  seven  sites  and  one 
weather  year  with  an  extremely  good  precision  (no  numerical  value  was  given  for  the 
goodness  of  fit).  However,  the  model  was  not  able  to  simulate  the  seasonal-course  of  leaf 
area  index  very  well,  suggesting  cautious  interpretation  of  the  achieved  predictions.  In 
addition,  no  real  independent  dataset  was  used  for  validating  the  calibrated  model  since 
the  presented  results  came  from  an  on-going  project.  It  was  suggested  that  the  capability 
of  simulating  run-on  and  run-off  at  a field  level  would  greatly  improve  water  balance 
simulation  and,  therefore,  crop  yield  prediction. 

Cora  et  al.  (1999)  calibrated  and  validated  CERES-Maize  for  a Michigan  location 
with  a water  table  present  throughout  the  growing  season.  The  calibration  consisted  of 
determining  the  genetic  coefficients  for  one  site  (out  of  55)  in  one  of  the  weather  years. 
The  model  was  then  validated  for  the  other  sites  for  the  same  calibration  weather  year  and 
for  all  sites  in  another  weather  year.  For  the  calibration  weather  year,  the  model 
predicted  yield  within  2%  of  measured  values  showing  only  a slight  under-prediction  for 
yield  over  9000  kg/ha.  For  the  next  year  (a  year  with  lower  than  normal  rainfall  during 
the  grain-filling  period)  the  model  greatly  under  predicted  grain  yield  (field  average  yield 
difference  of  7760  kg/ha)  for  the  whole  range  of  measured  values.  Once  the  authors  ran 
the  model  with  no  water  stress  simulation,  the  fit  improved  greatly  bringing  the  authors  to 
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the  conclusion  that  the  model  does  not  simulate  the  water  table  depth  and  upward  water 
flux  well  for  these  soil  conditions. 

Paz  et  al.  (1997)  used  CROPGRO-Soybean  (Hoogenboom  et  al.,  1994)  to  test  the 
hypothesis  that  soil-water  related  stresses  were  responsible  for  grain  yield  spatial 
variability  of  a soybean  field  in  Iowa.  They  concluded  that  variable  rooting  and  water 
table  depths  explained  69%  of  the  total  variability.  The  additional  variability  was 
hypothesized  to  be  due  to  variable  plant  population,  disease  intensity  and  differential 
rainfall/infiltration  ratios  across  the  field.  Examples  of  using  crop  simulation  models  for 
decision  support  in  SSM  are  Boone  and  Landivar  (1999),  Braga  et  al.  (1999),  Paz  et  al. 
(1999),  Bouma  et  al  (1997),  Booltink  and  Verhagen  (1997)  and  Van  Uffelen  et  al.  (1997). 
Objective 

The  objective  of  this  paper  is  to  test  the  ability  of  a crop  simulation  model  to 
predict  the  spatial  variability  of  corn  growth  and  yield  and  yield  components.  An 
additional  objective  is  to  study  the  relative  contribution  of  each  spatially  variable  input  to 
yield  variability  and  its  predictability.  Based  on  previous  studies,  we  hypothesize  that  for 
our  field  experiment  site,  factors  affecting  plant-available  soil-water  (soil-water  holding 
limits  and  effective  soil  depth)  and  plant  population  were  the  major  causes  for  spatial 
variability  of  crop  growth  and  yield. 

Material  and  Methods 

Field  Experiment 

This  study  was  based  on  the  field  experiment  fully  described  in  chapter  2.  The 
experiment  was  carried  out  in  a corn  farm  field  located  in  Michigan,  USA  during  the  crop 
seasons  of  1997  and  1998.  It  consisted  of  observing  variables  of  interest  in  a uniformly 
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managed  field.  The  sites  for  observing  soil  and  crop  variables  (total  of  43)  were  spatially 
distributed  in  a grid  of  approximately  30.4-m  (Figure  1).  The  crop  was  planted  on  April 
29  (DOY  1 19)  in  1997  and  on  May  1 1 (DOY  131)  in  1998. 

Effective  soil  depth  (simply  designated  as  soil  depth)  and  particle  size  analysis 
(carried  out  for  each  1 5-cm  increment  layer  up  to  the  constraining  clay  layer)  were 
determined  prior  to  planting  in  1997.  Plant  weight  and  leaf  area  index  were  determined 
at  anthesis  (DOY  21 1 in  1997  and  207  in  1998).  Grain  yield  and  its  components,  namely 
plant  population  (PP),  grain  number  per  ear  (G#EAR),  grain  number  per  unit  area 
(G#AREA)  and  unit  grain  weight  (UGW)  were  determined  at  maturity  (DOY  289  in  1997 
and  288  in  1998). 

The  field  average  yield  was  different  in  the  two  years:  In  1997  the  average  was 
9571  kg/ha  (range:  6.6  - 12.0  t/ha)  and  in  1998  the  average  was  11331  kg/ha  (range:  9.8 
- 12.6  t/ha).  Plant  population  varied  between  5.4  and  10.3  pl./m2  with  an  average  of  8.4 
pl./nr  in  1 997  and  between  7.4  and  1 0.3  pl/m2  with  an  average  of  9.0  pl/m2  in  1 998.  Soil- 
water  content  of  each  15-cm  increment  interval  (down  to  effective  soil  depth)  for  each 
site  was  determined  regularly  (12  times  between  DOY  170  and  290)  during  the  season 
using  a neutron  moisture  meter. 

Rainfall,  maximum  and  minimum  temperature  and  solar  radiation  were  recorded 
daily.  Rainfall  in  1 997  was  lower  than  in  1 998  during  the  crop  season  and  whole  year 
(Figure  2).  The  final  values  for  1997  and  1998  were  849.1  and  529.1  mm  for  the  year 
and  225.9  mm  and  436.1  mm  (difference  was  210.2  mm)  for  the  crop  season.  From  a 
weather  series  (1979-1988)  of  the  same  location,  the  average  seasonal  rainfall  amount 
was  459  mm  with  minimum  and  maximum  values  of  313  and  701  mm.  Considering 
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these  statistics,  1998  can  be  considered  an  average  year  in  terms  of  seasonal  rainfall 
whereas  1 997  was  a dry  one  (the  driest  in  the  series). 

The  soils  in  the  field  have  a profile  drainage  constraint  caused  by  a clay  layer 
located  at  variable  depths.  It  is  usually  shallower  in  zones  of  higher  elevation  and  deeper 
in  depressions.  Typically,  this  layer  causes  a perched  water  table  that  is  present  in  the 
profile  at  the  beginning  of  spring,  due  to  snow  melt,  and  at  the  end  of  the  summer  in  high 
rainfall  years  (summer  storms).  In  the  middle  of  the  season,  the  water  table  disappears 
once  depleted  by  roots.  The  water  table  plays  a very  importance  role  in  supplying  water 
during  critical  yield  forming  phases:  flowering  and  grain  filling.  In  years  with  higher 
rainfall,  the  water  table  is  present  throughout  the  season  at  higher  and  more  variable 
water  levels  causing  excess  water  stress  to  plants. 

Crop  Simulation  Model 

The  crop  simulation  model  used  for  this  study  was  CERES-Maize  V 3.5 
(Hoogenboom  et  al.,  1994).  A clear  description  of  the  model  principles  is  provided  by 
Ritchie  et  al.  (1998).  It  is  a process-based  model  that  uses  a comprehensive  approach  to 
crop  growth  and  development  and  their  environmental  determinant  factors  and 
interactions.  The  model  simulates  phenological  development  and  crop  growth  as  affected 
by  genetics,  management  and  environment  using  a daily  step.  The  clear  distinction 
between  crop  development  and  growth  in  the  model  is  a key  to  the  simulation  approach. 
Crop  development  is  simulated  using  the  "thermal  time"  concept  and  is  affected  by 
temperature,  genotype  and,  to  some  extent,  the  photoperiod  during  floral  induction.  Crop 
growth  is  simulated  using  the  "radiation  use  efficiency"  concept  and  is  affected  by 
genotype  and  environment.  The  environmental  factors  affecting  crop  growth  are 
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temperature,  solar  radiation,  soil  water,  nitrogen  and  aeration  status  and  root  system 
distribution  (described  by  the  root  distribution  factor,  SRGF). 

The  genetic  characteristics  of  each  cultivar  are  described  by  six  coefficients: 
thermal  time  from  seedling  emergence  to  the  end  of  the  juvenile  phase  (PI);  extent  to 
which  development  (days)  is  delayed  for  each  hour  increase  in  photoperiod  above  the 
longest  photoperiod  at  which  development  proceeds  at  a maximum  rate  (P2);  thermal 
time  from  silking  to  physiological  maturity  (P5);  maximum  possible  number  of  kernels 
per  plant  (G2);  kernel  filling  rate  (mg/day)  during  the  linear  grain  filling  stage  and  under 
optimum  conditions  (G3)  and  the  interval  in  thermal  time  between  successive  leaf  tip 
appearances  (PHINT). 

In  order  to  account  for  the  effects  of  soil  water  and  nitrogen  stresses,  the  model 
simulates  water  and  nitrogen  dynamics  in  the  soil.  Clear  descriptions  of  these  sub- 
models can  be  found  in  Ritchie  (1998)  and  Godwin  and  Singh  (1998). 

The  soil  water  redistribution  during  infiltration  is  simulated  using  the  "cascading" 
approach.  In  such  an  approach,  the  soil  is  divided  into  representative  layers  (usually 
more  than  five).  The  lower  soil  water  holding  limit  (LL),  the  drained  upper  limit  (DUL) 
and  saturation  (SAT)  parameters  describes  the  behavior  of  each  layer.  The  three  water- 
holding limits  correspond  to  the  classical  permanent  wilting  point,  field  capacity  and 
saturation.  Water  is  held  in  each  layer  if  it  is  below  the  drained  upper  limit  value.  If  soil 
water  content  surpasses  the  drained  upper  limit,  water  is  lost  due  to  drainage  to  the  next 
layer.  Drainage  of  each  layer  can  be  calculated  in  the  model  through  a saturated 
hydraulic  conductivity  (SSKS)  or  through  a daily  percent  of  the  difference  between 
saturation  and  drained  upper  limit  (SLDR).  By  setting  the  saturated  hydraulic 
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conductivity  of  a specified  layer  to  low  values,  it  is  possible  to  simulate  a perched  water 
table  above  that  layer. 

Infiltration  is  simulated  by  the  difference  between  rainfall  and  runoff.  Runoff  is 
simulated  using  an  adaptation  of  the  USDA-Soil  Conservation  Service  curve  number 
method  (Williams  et  al.,  1991).  Potential  evapotranspiration  is  simulated  based  on 
Ritchie  (1972)  and  Priestley  and  Taylor  (1972). 

The  model  provides  the  possibility  of  simulating  crop  growth  and  yield  without 
nitrogen  and/or  water  stresses.  This  is  done  simply  turning  off  the  corresponding  sub- 
models. This  capability  allows  the  user  to  use  only  those  sub-models  considered 
important  (ones  limiting  growth)  for  their  specific  conditions  avoiding  the  need  for 
unnecessary  inputs. 

Model  Inputs 

The  model  requires  inputs  relative  to  cultivar,  crop  management,  soil  and  weather. 
Crop  management  (day  of  planting,  plant  population  and  row  spacing)  and  weather  inputs 
(solar  radiation,  air  maximum  and  minimum  temperatures  and  rainfall)  were  obtained 
from  measured  data  described  previously.  Cultivar  genetic  inputs  were  obtained  from  an 
independent  dataset  (using  the  same  cultivar)  described  by  Cora  et  al.  (1999). 

Effective  soil  depth  for  all  sites  was  obtained  from  measured  data  described  in  the 
field  experiment  section.  SLDR  was  considered  to  be  0.5  day'1.  The  curve  number 
(SLRO)  was  assumed  to  be  homogenous  among  sites  in  the  field  based  on  direct 
observation  since  almost  no  runoff  occurred  in  either  year.  We  used  the  value  of  75 
based  on  Ritchie  et  al.  (1989). 

The  values  of  LL,  DUL,  SAT  and  SRGF  for  each  layer  down  to  effective  soil 
depth  and  SSKS  of  the  deepest  layer  were  estimated  using  an  estimation  procedure  fully 
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described  in  chapter  6.  This  procedure  is  based  on  the  reverse  modeling  process  and 
parameter  optimization  and  uses  the  observed  time  series  of  soil-water  content  as  the 
objective  function  variable. 

Simulation  Analysis 

In  order  to  test  the  ability  of  CERES-Maize  to  predict  the  spatial  variability  of 
corn  growth  and  yield  the  model  was  run  for  both  years  with  the  most  accurate  inputs  for 
each  site.  Predictions  of  corn  growth  and  yield  were  thus  independent  of  observed  values 
in  this  field  for  both  years.  The  agreement  between  predicted  and  observed  crop  growth 
(LAI  and  plant  weight  at  anthesis),  grain  yield  and  components  (G#EAR,  UGW  and 
G#AREA)  was  evaluated  by  the  root  mean  square  error  of  prediction  (RMSEP)  among 
sites.  A field-wide  evaluation  was  also  carried  out  through  visual  inspection  of  the  spatial 
pattern  of  yield  variation  obtained  through  kriging.  In  addition,  the  spatial  patterns  of  the 
differences  between  predicted  and  observed  yield  are  presented. 

To  study  the  relative  contribution  of  each  spatially  variable  input  (plant-available 
soil-water  and  plant  population)  the  model  was  run  for  both  years  with  site-specific  plant- 
available  soil-water  (soil-water  holding  limits  and  effective  soil  depth)  and  homogenous 
plant  population.  The  comparison  between  both  runs  was  made  by  comparing  the 
RMSEP  of  grain  yield  among  sites.  The  homogenous  plant  population  was  chosen  to  be 
the  spatial  average  of  the  measured  values  (8.4  pl./m2  in  1997  and  9.0  pl/m2  in  1998). 

Results  and  Discussion 
Model  Predictability  of  Growth  and  Yield 

Figure  37  shows  the  predicted  versus  observed  values  of  grain  yield.  The  RMSEP 
values  for  the  respective  relation  are  displayed  in  each  graph.  The  average  RMSEP  of 
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grain  yield  for  both  years  was  501  kg/ha,  which  is  an  exceptionally  low  value.  The 
proportions  of  yield  variability  explained  by  the  model  were  88%  and  76%  for  1997  and 
1998,  respectively.  The  better  performance  of  the  model  in  1997  compared  to  1998 
makes  sense  if  we  consider  that  1997  was  a dry  year.  The  model  did  not  explain  the 
remaining  variability  because  either  the  model  is  not  perfect  or  we  did  not  take  all  the 
spatially  variable  inputs  into  account  or  both. 

Although  the  RMSEP  values  for  grain  yield  were  quite  low,  in  1997  the  model 
showed  a tendency  to  under  estimate  higher  yields  (>10000  kg/ha)  and  in  1998  the  model 
showed  a tendency  to  over  estimate  yield  for  the  whole  range  of  observed  values. 

Figure  38  through  Figure  42  display  the  predicted  versus  observed  values  of  the 
yield  components  and  growth  variables.  The  model  performed  very  well  in  predicting 
these  variables  as  well.  The  model  only  performed  poorly  in  predicting  LAI  at  anthesis 
for  1998.  There  was  a clear  lack  of  sensitivity  of  the  model  because  lower  values  were 
over  estimated  and  higher  values  under  estimated.  Nevertheless,  this  error  did  not  affect 
the  overall  performance  of  the  model  because  all  observed  LAI  were  above  3.5,  which 
intercepts  most  of  the  incoming  solar  radiation. 


Predicted  (kg/ha)  Predicted  (kg/ha) 


85 


Figure  37.  Predicted  versus  observed  grain  yield  (kg/ha)  for  1997  (above)  and  1998 
(below). 


Predicted  Predicted 
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Figure  38.  Predicted  versus  observed  grain  number  per  ear  for  1997  (above)  and  1998 
(below). 
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Figure  39.  Predicted  versus  observed  unit  grain  weight  (mg)  for  1997  (above)  and  1998 
(below). 
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Figure  40.  Predicted  versus  observed  grain  number  per  unit  area  (#/m2)  for  1997  (above) 
and  1998  (below). 
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Figure  41 . Predicted  versus  observed  LAI  at  anthesis  for  1997  (above)  and  1998  (below). 
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Figure  42.  Predicted  versus  observed  plant  weight  at  anthesis  (kg/ha)  for  1997  (above) 
and  1998  (below). 
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The  fact  that  the  model  was  able  to  predict  grain  yield  as  well  as  yield 
components  and  growth  in  a very  satisfactory  way  increases  the  confidence  in  its 
performance.  A model  that  is  able  to  predict  yield  but  does  a poor  job  in  all  other 
variables  may  have  questionable  validity.  Since  a crop  simulation  model  is  a mechanistic 
tool,  all  intermediate  variables  should  be  correctly  predicted.  Only  in  this  case  can  we  be 
sure  that  the  model  works  appropriately. 

The  maps  of  predicted  grain  yield  for  1997  and  1998  are  presented  in  Figure  43 
and  can  be  compared  with  the  observed  yield  map  displayed  in  Figure  7.  The  respective 
maps  of  the  difference  between  predicted  and  observed  are  presented  in  Figure  44.  The 
maximum  possible  error  had  similar  magnitude  in  both  years  and  was  below  1000  kg/ha. 
As  for  the  variability  among  sites,  the  model  predicted  the  spatial  pattern  of  grain  yield  in 
a very  satisfactory  way  since  those  errors  are  perfectly  compatible  with  most  model 
applications. 

Factor  Contribution  to  Predictability  of  the  Model 

Figure  45  shows  the  simulated  versus  observed  values  of  grain  yield  for  1997  and 
1998  assuming  homogenous  plant  population  across  sites.  The  RMSE  values  for  the 
respective  relation  are  displayed  in  each  graph.  The  average  RMSE  of  grain  yield  over 
both  years  was  565  kg/ha.  The  proportion  of  yield  variability  explained  by  the  model 
were  83%  and  71%  for  1997  and  1998,  respectively.  These  percentages  represent  a 
decrease  of  only  5%  points  of  total  observed  variability  when  site-specific  plant 
population  was  not  considered. 


92 


Figure  43.  Map  of  predicted  grain  yield  (kg/ha)  for  1997  and  1998. 


Figure  44.  Map  of  the  difference  between  predicted  and  observed  grain  yield  (kg/ha)  of 
1997  (above)  and  of  1998  (below). 
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These  results  mean  that,  on  average,  93%  of  the  yield  spatial  variability  explained 
by  the  model  (82%)  was  due  to  spatially  variable  plant-available  soil-water  (soil-water 
holding  limits  and  effective  soil  depth)  and  7%  due  to  spatially  variable  plant  population. 
The  poor  contribution  of  plant  population  to  the  spatial  variability  of  grain  yield  is 
probably  due  to  the  attenuation  effect  introduced  by  the  compensation  among  yield 
components.  The  decreases  in  RMSEP  by  including  spatially  variable  plant  population  in 
addition  to  spatially  variable  plant-available  soil-water  were  only  256  kg/ha  and  1 12 
kg/ha  for  1997  and  1998,  respectively.  These  results  are  consistent  with  the  fact  that 
plant  population  was  lower  on  average  and  had  higher  spatial  variability  in  1 997  than  in 
1998. 


Conclusions 

CERES-Maize  V3.5  accurately  predicted  the  spatial  variability  of  com  growth 
and  yield  for  two  years  of  independent  data  when  accurate  inputs  were  provided  for  soil 
properties  and  plant  population.  The  average  root  mean  square  error  of  prediction  of  the 
model  was  501  kg/ha.  The  model  accounted  for  more  than  80%  of  total  observed  yield 
variability.  The  model  also  performed  very  well  in  predicting  most  of  the  yield 
components  and  growth  variables  as  well. 

The  model-based  approach  allowed  us  to  study  the  contribution  of  each  spatially 
variable  input  to  the  spatial  variability  of  yield.  For  the  studied  field,  the  main  factor 
responsible  for  spatial  yield  variability  was  plant-available  soil-water  (as  affected  by 
effective  soil  depth  and  soil-water  holding  capacity),  which  alone  explained  an  average  of 
93%  of  the  spatial  grain  yield  variability  accounted  by  the  model  (total  average  of  82%). 
Only  the  remaining  7%  were  accounted  by  plant  population. 


Simulated  Yield  (kg/ha)  Simulated  Yield  (kg/ha) 
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Figure  45.  Predicted  versus  observed  grain  yield  (kg/ha)  for  1997  (above)  and  1998 
(below)  assuming  homogenous  plant  population  across  sites. 


95 


These  excellent  results  indicate  that  the  model  can  be  used  for  further  SSM 
studies  when  inputs  are  accurately  known.  Nevertheless,  prudence  is  advised  in  the 
generalization  of  such  results  because  of  the  specificity  of  the  studied  field  situation:  an 
intensively  observed  field  experiment  with  no  stresses  other  than  those  related  to  soil- 
water.  In  addition,  the  two  years  of  experiment  did  not  include  a wet  year  where  excess 
water  stress  plays  an  important  role  in  shaping  final  grain  yield  spatial  pattern.  In  more 
common  situations  the  spatial  pattern  of  grain  yield  is  an  integration  of  multiple 
interacting  dynamic  stresses  caused  by  nitrogen  deficiency,  high  weed  densities,  disease 
incidence,  etc.  The  wide  use  of  crop  simulation  models  in  Site-specific  Management  is 
not  compatible  with  the  sampling  intensity  carried  in  this  field  to  obtain  inputs  and 
parameters  essential  for  the  accuracy  of  the  predictions.  Methods  to  obtain  site-specific 
model  inputs  and  parameters  other  than  direct  measurement  have  to  be  investigated  to 
allow  the  generalization  of  crop  simulation  models  use  in  Site-Specific  Management. 
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CHAPTER  6 

USING  OPTIMIZATION  TO  ESTIMATE  SOIL  INPUTS  OF  CROP  MODELS  FOR 
USE  IN  SITE-SPECIFIC  MANAGEMENT 

Introduction 

Crop  simulation  models  have  shown  great  potential  for  accurately  predicting 
within-field  spatial  variability  of  grain  yield  (Chapter  5 in  this  volume;  Cora  et  al.,  1999; 
Fraisse  et  al.,  1999;  Paz  et  al.,  1997).  Predicting  the  spatial  variability  of  grain  yield  is  of 
crucial  importance  for  Site-Specific  Management  (SSM)  because  it  allows  for  testing  of 
management  prescriptions  and  for  correct  assessment  of  agronomic  and  economic 
outcomes  (Sudduth  et  al.,  1996).  Nevertheless,  one  common  limitation  of  crop 
simulation  model  use  in  SSM  is  the  need  for  accurate  values  of  many  inputs  for 
numerous  sites  in  a field.  This  limitation  has  led,  in  part,  to  some  unsuccessful 
applications  of  crop  simulation  models  to  predict  observed  yield  variability  (e.g.  Sadler  et 
al.,  1995;  Sudduth  et  al.,  1998). 

Some  model  inputs,  such  as  soil-water  holding  limits  of  each  soil  layer  or 
saturated  hydraulic  conductivity  of  any  drainage  restrictive  layer,  are  difficult  or 
expensive  to  determine,  but  are  of  crucial  importance  to  model  performance  in  non- 
irrigated  conditions.  It  is,  therefore,  extremely  important  that  we  measure  or  estimate 
them  correctly.  Laboratory  procedures  do  not  always  give  the  best  estimates  and  are  not 
as  inexpensive  and  practical  as  needed  for  site-specific  management.  Field  procedures 
are  extremely  laborious  and  time  consuming.  A third  way  to  obtain  model  inputs 
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involves  optimization  for  the  estimation  of  parameters  using  more  easily  measured 
variables. 

Site-Specific  Parameters  Estimation  Using  Optimization 

In  the  optimization  approach,  a framework  is  established  where  all  model  inputs, 
with  the  exception  of  the  ones  to  be  estimated,  are  considered  to  be  known.  For 
estimating  soil  water  holding  limits  in  a spatially  variable  field,  cultivar  inputs,  weather 
inputs,  and  soil  parameters  other  than  those  being  estimated  are  assumed  to  be  accurately 
known.  Also,  it  is  assumed  that  the  present  model  accurately  simulates  yield  when  inputs 
are  known  and  accurate.  With  these  assumptions,  the  lack  of  fit  between  observed  and 
simulated  data  is  attributed  to  non-optimal  estimation  of  the  spatially  variable  inputs. 

The  goodness  of  fit  is  measured  by  the  difference  between  simulated  and  observed 
variable  or  variables  considered  being  most  important  for  the  particular  study. 
Subsequently,  new  input  estimations  are  obtained  in  an  attempt  to  improve  the  fit. 

The  estimated  inputs  are  not  totally  independent  of  the  other  inputs  and  of  the 
model  parameters  used  in  the  procedure.  Therefore,  to  evaluate  the  input  estimates 
obtained  from  this  method,  a completely  independent  dataset  is  still  needed  for 
validation. 

The  procedure  to  estimate  the  new  parameters  can  be  of  various  forms  with 
different  computational  power  and  speed  requirements.  These  search  techniques  probe 
through  a multidimensional  space  of  feasible  values  based  upon  a certain  criterion. 
Generally,  the  criterion  is  a goodness-of-fit  measure  such  as  the  root  mean  square  error 
(RMSE).  The  major  difference  among  different  input-search  techniques  lies  in  the 
method  to  select  the  new  estimates.  They  can  vary  from  simple  systematic  search 
processes  (usually  designated  as  grid  search)  to  more  elaborate  methods,  such  as  genetic 
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algorithms  and  simulated  annealing-based  algorithms.  Press  et  al.  (1992)  provides  an 
introduction  to  all  these  methods  and  specific  references.  There  are  few  recipes  to  decide 
which  procedure  works  better  for  each  case  since  “better”  means  a trade-off  between 
computational  time,  problem  complexity  and  guarantee  of  the  best  or  optimal  solution. 

Although  the  search  space  of  each  individual  parameter  being  estimated  is  pre- 
defined according  to  feasible  values,  there  is  the  possibility  that  the  combination  of 
parameters  that  provides  the  best  goodness  of  fit  is  totally  unrealistic.  For  this  reason,  it 
is  a good  practice  to  perform  a validation  of  the  estimated  parameters. 

In  addition  to  providing  site-specific  crop  model  inputs,  this  procedure  allows  one 
to  test  hypotheses  relating  the  most  important  causes  of  spatial  variability.  The  rationale 
for  this  use  is  “if  the  estimated  soil-water  holding  limits  provide  an  acceptable  goodness 
of  fit  between  simulated  and  observed  variables  then  soil-water  related  stresses  are  the 
most  limiting  processes  for  that  environment  and,  therefore,  the  causes  of  variability”. 
This  reasoning  assumes  that  other  model  inputs  are  spatially  constant.  Paz  et  al.  (1997) 
provides  a good  example  of  this  application  of  crop  model  input  optimization.  They  used 
CROPGRO-Soybean  (Hoogenboom  et  al.,  1994)  to  test  the  hypothesis  that  soil-water 
related  stresses  were  responsible  for  grain  yield  spatial  variability  of  a soybean  field  (16 
ha)  in  Iowa.  After  estimating  profile  drainage  rate  and  rooting  depth  through 
optimization  (using  final  grain  yield  as  the  objective  function  variable),  the  model 
explained  69%  of  the  observed  yield  spatial  variability. 

Once  site-specific  crop  model  parameters  are  obtained  in  this  way  for  a specific 
field  and  validated,  they  can  be  used  to  run  the  crop  model  to  test  management 
prescriptions  and  to  correctly  assess  their  agronomic  and  economic  outcomes.  Such 
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procedures  could  be  used  by  industry  in  their  efforts  to  use  crop  simulation  models  in 
site-specific  management  programs  because  databases  of  multiple-year  data  (yield  map, 
weather,  management,  etc)  from  specific  customers'  fields  are  being  built.  Each  new  year 
of  data  added  to  the  database  at  the  end  of  the  crop  season  could  be  used  to  increase  the 
robustness  of  the  optimized  crop  model  inputs,  further  improving  the  accuracy  of  its 
prediction. 

There  is  one  serious  question  regarding  the  validity  of  site-specific  parameter 
estimation  concerning  the  model  variable  used  to  judge  the  goodness  of  fit.  Final  grain 
yield  may  be  the  ultimate  variable  of  interest,  but  the  objective  function  of  the  search 
technique  should,  as  much  as  possible,  use  a model  variable(s)  that  is  (are)  directly 
related  to  the  parameters  being  estimated.  For  example,  the  time  series  of  measured  soil- 
water  content  should  be  used  to  estimate  soil-water  holding  limits  in  the  optimization 
procedure,  and  not  the  final  grain  yield  values.  If  final  grain  yield  is  used  as  the  objective 
function  variable,  the  result  could  be  a simulation  that  merely  reproduces  grain  yield 
values  in  each  location  in  a field  at  the  cost  of  unrealistic  parameter  estimates. 

Conversely,  when  the  time  series  of  measured  soil-water  content  is  used  as  the  objective 
function  variable  the  result  is  a simulation  that  correctly  simulates  final  grain  yield  “for 
the  correct  reasons”.  Nevertheless,  the  time  series  of  measured  soil-water  content,  or 
most  of  the  intermediate  simulated  variables,  is  rarely  available,  especially  for  industry 
application.  Therefore,  the  ability  a crop  simulation  model  to  predict  the  spatial 
variability  of  grain  yield  using  site-specific  soil  parameters  estimated  from  optimization 
based  on  final  grain  yield  must  be  tested  and  compared  with  the  same  ability  resultant 
from  using  other  objective  function  variables. 
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Objective 

The  objective  of  this  paper  is  to  compare  the  ability  of  a crop  simulation  model  to 
predict  the  spatial  variability  of  grain  yield  using  site-specific  soil  parameters  estimated 
from  optimization  based  on  two  distinct  objective  function  variables.  The  objective 
function  variables  were  the  final  grain  yield  and  the  time  series  of  soil-water  content. 

The  comparison  aims  to  answer  the  following  questions:  Do  estimated  soil-water  holding 
limits  based  on  grain  yield  as  the  objective  function  variable  provide  an  accurate 
simulation  of  soil-water  content?  Do  estimated  soil-water  holding  limits  based  on  time 
series  of  soil-water  content  as  the  objective  function  variable  provide  good  simulation  of 
grain  yield? 

A negative  answer  to  the  second  question  independently  of  the  answer  to  the  first 
question  means  that  crop  model  should  be  revised.  The  optimal  scenario  is  a positive 
answer  to  both  questions  because  it  implies  that,  if  input  estimations  are  valid,  an  easily 
measured  variable  such  as  grain  yield  could  be  used  to  estimate  site-specific  soil 
parameters  for  crop  simulation.  If  however,  only  the  second  question  has  an  affirmative 
answer,  and  input  estimations  are  valid,  then  only  the  time  series  of  soil-water  content 
could  be  used  to  estimate  soil-water  holding  limits  and  the  extensive  application  of  the 
parameters  estimation  procedure  is  endangered. 

Material  and  Methods 

Field  Experiment 

This  study  was  based  on  the  field  experiment  fully  described  in  chapter  2.  The 
experiment  was  carried  out  in  a com  farm  field  located  in  Michigan,  USA  during  the  crop 
seasons  of  1997  and  1998.  It  consisted  of  observing  variables  of  interest  in  a uniformly 
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managed  field.  The  sites  for  observing  soil  and  crop  variables  (total  of  43)  were  spatially 
distributed  in  a grid  of  approximately  30.4-m  (Figure  1).  The  crop  was  planted  on  April 
29  (DOY  1 1 9)  in  1 997  and  on  May  1 1 (DOY  1 3 1 ) in  1 998. 

Effective  soil  depth  (simply  designated  as  soil  depth)  and  particle  size  analysis 
(carried  out  for  each  15-cm  increment  layer  up  to  the  constraining  clay  layer)  were 
determined  prior  to  planting  in  1997.  Plant  weight  and  leaf  area  index  were  determined 
at  anthesis  (DOY  21 1 in  1997  and  207  in  1998).  Grain  yield  and  its  components,  namely 
plant  population  (PP),  grain  number  per  ear  (G#EAR),  grain  number  per  unit  area 
(G#AREA)  and  unit  grain  weight  (UGW)  were  determined  at  maturity  (DOY  289  in 
1997  and  288  in  1998). 

The  field  average  yield  was  different  in  the  two  years:  In  1997  the  average  was 
9571  kg/ha  (range:  6.6  - 12.0  t/ha)  and  in  1998  the  average  was  11331  kg/ha  (range:  9.8 
- 12.6  t/ha).  Plant  population  varied  between  5.4  and  10.3  pl./m2  with  an  average  of  8.4 
pl./m2  in  1 997  and  between  7.4  and  1 0.3  pl/m2  with  an  average  of  9.0  pl/m2  in  1 998.  Soil- 
water  content  of  each  15-cm  increment  interval  (down  to  effective  soil  depth)  for  each 
site  was  determined  regularly  (12  times  between  DOY  170  and  290)  during  the  season 
using  a neutron  moisture  meter. 

Rainfall,  maximum  and  minimum  temperature  and  solar  radiation  were  recorded 
daily.  Rainfall  in  1997  was  lower  than  in  1998  during  the  crop  season  and  whole  year 
(Figure  2).  The  final  values  for  1997  and  1998  were  849.1  and  529.1  mm  for  the  year 
and  225.9  mm  and  436.1  mm  (difference  was  210.2  mm)  for  the  crop  season.  From  a 
weather  series  (1979-1988)  of  the  same  location,  the  average  seasonal  rainfall  amount 
was  459  mm  with  minimum  and  maximum  values  of  313  and  701  mm.  Considering 
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these  statistics,  1998  can  be  considered  an  average  year  in  terms  of  seasonal  rainfall 
whereas  1997  was  a dry  one  (the  driest  in  the  series). 

The  soils  in  the  field  have  a profile  drainage  constraint  caused  by  a clay  layer 
located  at  variable  depths.  It  is  usually  shallower  in  zones  of  higher  elevation  and  deeper 
in  depressions.  Typically,  this  layer  causes  a perched  water  table  that  is  present  in  the 
profile  at  the  beginning  of  spring,  due  to  snow  melt,  and  at  the  end  of  the  summer  in  high 
rainfall  years  (summer  storms).  In  the  middle  of  the  season,  the  water  table  disappears 
once  depleted  by  roots.  The  water  table  plays  a very  importance  role  in  supplying  water 
during  critical  yield  forming  phases:  flowering  and  grain  filling.  In  years  with  higher 
rainfall,  the  water  table  is  present  throughout  the  season  at  higher  and  more  variable 
water  levels  causing  excess  water  stress  to  plants. 

Crop  Simulation  Model  and  Inputs 

This  study  was  based  on  the  crop  simulation  model  CERES-Maize  V 3.5 
(Hoogenboom  et  al.,  1994)  used  in  chapter  5.  A clear  description  of  the  model  principles 
is  provided  by  Ritchie  et  al.  (1998).  It  is  a process-based  model  that  uses  a 
comprehensive  approach  to  crop  growth  and  development  and  their  environmental 
determinant  factors  and  interactions. 

The  model  requires  inputs  relative  to  cultivar,  crop  management,  soil  and  weather. 
Crop  management  (day  of  planting,  plant  population  and  row  spacing)  and  weather  inputs 
(solar  radiation,  air  maximum  and  minimum  temperatures  and  rainfall)  were  obtained 
from  measured  data  described  previously.  Cultivar  genetic  inputs  were  obtained  from  an 
independent  dataset  (using  the  same  cultivar)  described  by  Cora  et  al.  (1999). 

Effective  soil  depth  for  all  sites  was  obtained  from  measured  data  described 
previously.  SLDR  was  considered  to  be  0.5  day'1.  The  curve  number  (SLRO)  was 
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assumed  to  be  homogenous  among  sites  in  the  field  based  on  direct  observation  since 
almost  no  runoff  occurred  in  both  years.  We  used  the  value  of  75  based  on  Ritchie  et  al. 
(1989). 

The  values  of  LL,  DUL,  SAT  and  SRGF  for  each  layer  down  to  effective  soil 
depth  and  SSKS  of  the  deepest  layer  were  estimated  using  the  optimization-based  input 
estimation  procedure  that  is  the  focus  of  this  chapter.  The  details  about  the  input  search 
technique  and  the  objective  function  variables  are  given  in  next  section. 

Input  Search  Technique  and  Objective  Function  Variables 

The  choice  of  soil  parameters  to  estimate  for  each  particular  site  in  the  field  (LL, 
DUL,  SAT  and  SRGF  for  each  layer  and  SSKS  for  bottom  layer)  was  based  on  previous 
studies  showing  plant-available  soil-water  to  be  the  main  factor  responsible  for  spatial 
variability  of  yield  (previous  chapters).  The  chosen  soil  parameters  of  the  crop  model 
were  the  ones  that  mainly  affect  plant-available  soil-water  (with  exception  to  effective 
soil  depth,  which  was  measured). 

A simulated  annealing  algorithm  (Goffe  et  al.  1994)  was  used  as  the  parameter 
search  technique  for  minimizing  the  objective  function  for  each  of  the  variables  of 
interest.  The  objective  function  was  the  root  mean  square  error  when  the  variable  of 
interest  was  either  final  yield  or  time-series  of  soil-water  content.  Error  was  calculated  as 
the  arithmetic  difference  between  observed  and  simulated  values. 

The  range  of  search  for  the  soil-water  holding  limits  were  from  0.030  to  0.160 
cm3cm  '3  for  lower  limit,  0.120  to  0.350  cm3cnf3  for  drained  upper  limit  and  0.320  to 
0.500  cm3cm'3  for  saturation.  These  ranges  were  considered  reasonable  for  the  textures 
at  all  sites  and  at  all  layers  based  on  the  pedotransfer  function  developed  by  Saxton  et  al. 
(1986)  and  validated  by  Kern  (1995).  Impossible  combinations  of  estimates  were 
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eliminated  by  attributing  to  them  an  extraordinary  high  objective  function  value.  The 
range  of  values  for  the  bottom  layer  SSKS  was  set  from  0.000  to  0. 1 50  cm  h'1 . This  gave 
a full  range  of  options  from  a total  impervious  layer  to  a moderately  drained  layer.  The 
initial  values  for  the  SRGF  were  obtained  as  a function  of  layer  depth  according  to 
Ritchie  et  al.  (1989). 

Initial  values  for  the  soil-water  holding  limits  were  obtained  from  the  particle  size 
analysis  results  using  the  pedotransfer  function  existing  in  DSSAT  V3.5  (Ritchie  et  al., 
1989).  Since  simulated  annealing  is  a global  search  technique,  there  was  no  need  to  use 
multiple  initial  values  for  the  soil-water  holding  limits.  The  ability  to  find  the  global 
solution  was  tested  in  a representative  example  with  the  highest  possible  complexity 
(higher  number  of  search  values).  The  final  value  vector  obtained  from  the  algorithm 
was  compared  with  the  one  from  a systematic  search.  The  results  were  the  same  at  a 
three  significant  digit  precision  and  independent  of  the  initial  values.  This  indicated  that 
simulated  annealing  successfully  found  the  global  optimal  solution.  Because  of  this 
result  we  opted  to  use  the  soil-water  holding  limits  from  DSSAT  V3.5  as  initial  values 
since  these  would  be  the  values  used  as  default  by  most  model  users. 

The  validation  of  the  soil  parameters  generated  by  the  site-specific  input 
estimation  procedure  was  done  by  comparing  the  estimates  generated  by  the  procedure 
with  the  ones  estimated  from  the  pedotransfer  function  developed  by  Saxton  et  al.  (1986). 

For  both  objective  function  variables,  1997  and  1998  experimental  years  were 
used  as  both  estimation  and  validation  data  sets  alternately.  This  setup  provided  a 
factorial  computational  experiment  with  two  objective  function  variables  and  two  years 
targeted  as  estimation  and  validation  data  sets.  Site-specific  soil  parameters  were 
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estimated  for  all  of  the  43  sites  described  in  the  experimental  design.  The  estimation  and 
validation  data  sets  were  evaluated  through  the  RMSE  (RMSEF  and  RMSEP, 
respectively)  among  sites.  A field-wide  evaluation  was  also  carried  out  through  visual 
inspection  of  the  spatial  pattern  of  yield  variation  obtained  through  kriging.  In  addition, 
the  spatial  patterns  of  the  difference  between  predicted  and  observed  yield  are  presented. 

Results  and  Discussion 
Grain  Yield  as  Objective  Function  Variable 

Estimating  soil  parameters  using  grain  yield  as  the  objective  function  variable  led 
to  minimal  differences  between  simulated  and  observed  values  in  the  estimation  data  set 
(Table  14  and  Table  15).  This  confirmed  that  the  ranges  set  for  the  soil  parameters 
accounted  for  the  total  grain  yield  spatial  variability  existing  in  the  estimated  sets.  When 
these  soil  parameters  were  used  to  predict  yield  in  the  independent  year  (validation  data 
set),  the  average  RMSE  was  765  kg/ha  (683.8  kg/ha  for  1997  and  845.2  kg/ha  for  1998, 
Table  13).  The  predicted  versus  observed  yield  graphs  for  the  validation  data  set  are 
shown  in  Figure  47.  Although  the  RMSE  was  low  for  both  years,  a slight  tendency  to 
under  estimate  higher  yields  can  be  observed  for  1997.  The  resultant  maps  of  grain  yield 
are  displayed  in  Figure  48  and  can  be  compared  with  the  observed  yield  map  presented  in 
Figure  7.  The  maps  of  RMSEP  for  both  years  are  displayed  in  Figure  49.  Zones  of 
higher  error  (>  1000  kg/ha)  are  clearly  more  abundant  in  1998.  In  this  year  errors  higher 
than  1000  kg/ha  occurred  in  three  zones  whereas  in  1997  that  almost  never  happened. 
Despite  these  differences  the  errors  in  both  years  were  small  and  comparable. 
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Table  1 1 . Parameter  estimates  for  lower  limit  (LL,  cm3/cm3),  drained  upper  limit 
(DUL,  cm3/ cm3)  and  saturation  (SAT,  cm3/cm3)  for  site  2 resultant  from  the 
DSSAT  and  Saxton  pedotransfer  functions  and  from  the  optimization  based 
procedure  using  yield  (yield-based)  and  soil  water  content  (SWC-based)  as 
objective  functions. 


DSSAT 

Saxton 

Depth 

LL 

DUL 

SAT 

LL 

DUL 

SAT 

30  cm 

0.130 

0.261 

0.411 

0.129 

0.279 

0.478 

45  cm 

0.134 

0.270 

0.430 

0.130 

0.298 

0.488 

60  cm 

0.131 

0.268 

0.432 

0.127 

0.297 

0.487 

75  cm 

0.132 

0.265 

0.417 

0.129 

0.286 

0.481 

90  cm 

0.103 

0.226 

0.373 

0.109 

0.227 

0.443 

105  cm 

0.118 

0.239 

0.378 

0.123 

0.235 

0.454 

120  cm 

0.093 

0.213 

0.360 

0.100 

0.209 

0.428 

135  cm 

0.112 

0.233 

0.374 

0.117 

0.229 

0.448 

150  cm 

0.077 

0.194 

0.342 

0.085 

0.179 

0.402 

Yield- 

Yield- 

Based  97 

Based  98 

Depth 

LL 

DUL 

SAT 

LL 

DUL 

SAT 

30  cm 

0.113 

0.281 

0.436 

0.114 

0.276 

0.425 

45  cm 

0.117 

0.317 

0.376 

0.129 

0.279 

0.436 

60  cm 

0.117 

0.305 

0.480 

0.118 

0.271 

0.421 

75  cm 

0.119 

0.303 

0.420 

0.132 

0.286 

0.424 

90  cm 

0.093 

0.264 

0.461 

0.093 

0.237 

0.379 

105  cm 

0.094 

0.249 

0.440 

0.109 

0.239 

0.368 

120  cm 

0.098 

0.213 

0.354 

0.074 

0.218 

0.359 

135  cm 

0.069 

0.270 

0.364 

0.112 

0.241 

0.380 

150  cm 

0.102 

0.272 

0.448 

0.070 

0.212 

0.336 

SWC- 

SWC- 

Based  97 

based  98 

Depth 

LL 

DUL 

SAT 

LL 

DUL 

SAT 

30  cm 

0.121 

0.279 

0.479 

0.126 

0.284 

0.479 

45  cm 

0.130 

0.318 

0.425 

0.131 

0.312 

0.436 

60  cm 

0.124 

0.308 

0.424 

0.124 

0.301 

0.420 

75  cm 

0.129 

0.291 

0.426 

0.128 

0.290 

0.412 

90  cm 

0.100 

0.233 

0.378 

0.104 

0.233 

0.373 

105  cm 

0.112 

0.239 

0.394 

0.123 

0.240 

0.382 

120  cm 

0.101 

0.211 

0.374 

0.103 

0.212 

0.370 

135  cm 

0.101 

0.231 

0.366 

0.110 

0.230 

0.377 

150  cm 

0.066 

0.205 

0.414 

0.081 

0.212 

0.406 

107 


Table  12.  Parameter  estimates  for  lower  limit  (LL,  cm3/cm3),  drained 
upper  limit  (DUL,  cm3/cm3)  and  saturation  (SAT,  cm3/cm3)  for  site  42 
resultant  from  the  DSSAT  and  Saxton  pedotransfer  functions  and  from 
the  optimization  based  procedure  using  yield  (yield-based)  and  soil  water 


content  (SWC-based)  as  objective  functions. 


DSSAT 

Saxton 

LL 

DUL 

SAT 

LL 

DUL 

SAT 

30  cm 

0.094 

0.214 

0.361 

0.102 

0.209 

0.429 

45  cm 

0.109 

0.233 

0.379 

0.114 

0.236 

0.450 

60  cm 

0.100 

0.218 

0.363 

0.107 

0.210 

0.433 

75  cm 

0.980 

0.214 

0.361 

0.105 

0.200 

0.427 

90  cm 

0.097 

0.217 

0.362 

0.105 

0.211 

0.432 

105  cm 

0.049 

0.157 

0.331 

0.076 

0.160 

0.387 

Yield- 

Yield- 

Based  97 

Based  98 

LL 

DUL 

SAT 

LL 

DUL 

SAT 

30  cm 

0.090 

0.215 

0.352 

0.115 

0.231 

0.407 

45  cm 

0.109 

0.231 

0.368 

0.097 

0.274 

0.392 

60  cm 

0.103 

0.219 

0.340 

0.064 

0.240 

0.371 

75  cm 

0.099 

0.215 

0.359 

0.068 

0.236 

0.326 

90  cm 

0.097 

0.217 

0.470 

0.087 

0.201 

0.344 

105  cm 

0.031 

0.166 

0.322 

0.044 

0.173 

0.442 

SWC- 

SWC- 

Based  97 

based  98 

LL 

DUL 

SAT 

LL 

DUL 

SAT 

30  cm 

0.100 

0.211 

0.430 

0.098 

0.210 

0.433 

45  cm 

0.113 

0.239 

0.383 

0.111 

0.235 

0.384 

60  cm 

0.108 

0.211 

0.371 

0.107 

0.209 

0.370 

75  cm 

0.107 

0.204 

0.351 

0.106 

0.203 

0.348 

90  cm 

0.108 

0.213 

0.362 

0.108 

0.210 

0.370 

105  cm 

0.074 

0.158 

0.378 

0.061 

0.159 

0.386 
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Figure  46.  Map  of  total  profile  plant  available  soil-water  (mm)  resultant  from  soil 
parameters  estimated  from  yield  of  1997  (above)  and  1998  (below). 


Table  13.  Root  mean  square  error  of  grain  yield  (kg/ha)  for  the  soil  parameter  estimation 
strategies.  The  average  relative  error  (%)  is  presented  between  parentheses. 


Grain  yield  based  estimation  Soil-water  based  estimation 

Estimation  set  Validation  set  Estimation  set  Validation  set 

Estimation  for  97 

6.1 

845.2 

472.1 

444.6 

Validation  for  98 

(0.1) 

(8.1) 

(4.5) 

(4-2) 

Estimation  for  98 

6.5 

683.8 

412.5 

555.6 

Validation  for  97 

(0.1) 

(6.5) 

(4.0) 

(5.3) 

Average 

6.3 

764.5 

442.3 

500.1 

(0.1) 

(7.3) 

(4.2) 

(4.8) 
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Table  14.  Observed  and  simulated  grain  yield  (kg/ha)  for  each  of  the  estimation 
procedures  for  1 997. 

site 

Observed  97 

Yield-based  97 

Yield-based  98 

SWC-Based  97 

SWC-Based  98 

1 

8877 

8877 

8138 

8457 

8555 

2 

11671 

11632 

10658 

11092 

10714 

3 

11963 

11963 

10835 

11301 

11301 

4 

11415 

11415 

10921 

10715 

10543 

5 

9458 

9458 

9805 

9583 

8756 

6 

11946 

11946 

11006 

10996 

10996 

7 

10946 

10946 

9638 

10253 

10026 

8 

11423 

11423 

10561 

10696 

10696 

9 

10114 

10114 

9970 

9686 

9578 

10 

10480 

10480 

9968 

9741 

9741 

1 1 

10123 

10123 

10193 

9647 

9647 

12 

9547 

9547 

8483 

9279 

9341 

13 

10431 

10431 

10264 

10123 

9993 

14 

10512 

10512 

9283 

9943 

9764 

15 

8748 

8740 

9330 

8537 

9238 

16 

11345 

11345 

10536 

10700 

10363 

17 

10543 

10543 

9502 

10574 

10487 

18 

9716 

9716 

9783 

8978 

9522 

19 

11345 

11345 

10483 

10728 

10517 

20 

8282 

8282 

8138 

8935 

8769 

21 

9246 

9246 

8549 

9080 

8965 

22 

9189 

9189 

10237 

9805 

9798 

23 

10548 

10548 

10368 

9919 

10492 

24 

6620 

6620 

7370 

7058 

6294 

25 

8735 

8735 

8381 

8968 

8758 

26 

8937 

8937 

8321 

8985 

8949 

27 

10632 

10632 

10697 

10696 

10310 

28 

9774 

9774 

10271 

9910 

9716 

29 

8747 

8747 

8707 

9104 

9101 

30 

9904 

9904 

9299 

9895 

9388 

31 

8986 

8986 

9633 

9453 

9397 

32 

8368 

8368 

8340 

8109 

7918 

33 

8347 

8347 

7876 

8063 

7929 

34 

8764 

8764 

8975 

8779 

7856 

35 

8286 

8286 

8532 

8755 

8988 

36 

7345 

7345 

8397 

7736 

7799 

37 

8647 

8647 

8320 

8264 

8518 

38 

8579 

8579 

8649 

8744 

8643 

39 

8062 

8062 

8960 

8122 

7698 

40 

9560 

9560 

9695 

9382 

9498 

41 

8756 

8756 

9439 

8837 

8900 

42 

6800 

6800 

7328 

7512 

7173 

43 

9843 

9843 

8883 

9299 

8938 
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Table  15.  Observed  and  simulated  grain  yield  (kg/ha)  for  each  of  the  estimation 
procedures  for  1998. 


site 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 


Observed  98  Yield-based  98  Yield-based  97  SWC-Based  98  SWC-Based  97 


10795 

10795 

10824 

10802 

10367 

12366 

12366 

13451 

12811 

13155 

12629 

12629 

13475 

13108 

13108 

12531 

12531 

12825 

12178 

12445 

11593 

11593 

10912 

11346 

11833 

12626 

12626 

13342 

12962 

12962 

11114 

11114 

12430 

10972 

11627 

12140 

12140 

12948 

12339 

12339 

11433 

11433 

11335 

10811 

10801 

11605 

11605 

12672 

11729 

11729 

12205 

12205 

11390 

11973 

11973 

10334 

10334 

11379 

10542 

10688 

11982 

11982 

11445 

12198 

12244 

11170 

11170 

12253 

11073 

11481 

11433 

11425 

9961 

11768 

11058 

12094 

12094 

12815 

11757 

12238 

11037 

11037 

11841 

11411 

12149 

11763 

11793 

11123 

12036 

11749 

11798 

11798 

12645 

12268 

12329 

10216 

10216 

10447 

10413 

10900 

10939 

10939 

12206 

11363 

11645 

12152 

12152 

11180 

12575 

12093 

12155 

12155 

12100 

12498 

12375 

10231 

10231 

8766 

9547 

10210 

10829 

10829 

10234 

10764 

10931 

10871 

10871 

11135 

10850 

10877 

12416 

12416 

11954 

12538 

12711 

12242 

12242 

12687 

12481 

12578 

10961 

10961 

10391 

11255 

11011 

11145 

11145 

12018 

12146 

12290 

11396 

11396 

11291 

11420 

11409 

10639 

10639 

9691 

10219 

10517 

9940 

9940 

10551 

10209 

10559 

11154 

11125 

9779 

10218 

10444 

10769 

10769 

9806 

10989 

11032 

10518 

10518 

9313 

10512 

10643 

10312 

10312 

11263 

10810 

10736 

10686 

10686 

10568 

10585 

11191 

11089 

11089 

9950 

10158 

10663 

11520 

11520 

12339 

12156 

11715 

11562 

11562 

11668 

11143 

11239 

9834 

9834 

9000 

10209 

10340 

11045 

11045 

11901 

10661 

11033 

Predicted  Yield  (kg/ha)  Predicted  Yield  (kg/ha) 
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Figure  47.  Predicted  versus  observed  grain  yield  for  1997  using  soil  parameters 
estimated  from  yield  of  1998  (above)  and  for  1998  using  soil  parameters  estimated  from 
yield  of  1997  (below). 
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Figure  48.  Map  of  predicted  com  yield  (kg/ha)  for  1997  using  soil  parameters  estimated 
from  yield  of  1998  (above)  and  for  1998  using  soil  parameters  estimated  from  yield  of 
1997  (below). 


Figure  49.  Map  of  the  difference  between  predicted  and  observed  grain  yield  (kg/ha)  for 
1997  using  soil  parameters  estimated  from  yield  of  1998  (above)  and  for  1998  using  soil 
parameters  estimated  from  yield  of  1997  (below). 
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Table  1 1 and  Table  12  display  the  soil  parameter  estimates  for  LL,  DUL  and  SAT 
for  sites  2 and  42  resultant  from  the  DSSAT  and  Saxton  pedotransfer  functions  and  from 
the  optimization  procedure.  These  two  sites  had  the  lowest  and  the  highest  values  of 
estimated  profile  plant  available  soil  water,  representing  the  two  extreme  situations 
(Figure  46). 

Figure  50  and  Figure  51  show  the  estimated  values  of  LL,  DUL,  and  SAT  and  of 
the  differences  between  DUL  and  LL  and  between  SAT  and  DUL  against  the  estimations 
of  the  same  quantities  by  the  Saxton’s  pedotransfer  function.  The  estimated  LL  tended  to 
be  consistently  below  the  expected  values  obtained  from  Saxton’s  pedotransfer  function 
for  LL  lower  than  0. 100  cm3/cm3.  The  same  occurred  for  SAT  but  for  the  whole  range 
of  values.  Estimated  DUL  values  agreed  well  with  Saxton’s  estimates.  Consequently, 
the  grain  yield  based  estimation  resulted  in  larger  plant-available  soil-water  values  for 
sites  with  lower  LL.  The  differences  between  SAT  and  DUL  were  lower  for  estimated 
values  than  for  Saxton’s  estimates.  Only  in  a few  profiles,  the  estimated  saturated 
hydraulic  conductivity  allowed  the  creation  of  a perched  water  table  (9  sites  for 
calibration  in  1997  and  5 sites  for  calibration  in  1998).  Instead,  the  plant  availability  of 
soil-water  to  plant  was  made  possible  by  unrealistic  low  values  of  LL. 

Table  16  displays  the  root  mean  square  error  and  the  relative  error  of  soil-water 
content  by  soil  layer  depth.  The  magnitude  of  the  error  was  consistent  for  the  two  years 
and  varied  according  to  layer.  Up  to  90  cm,  the  average  relative  error  was  always  below 
20%  and  increased  to  45%  for  deeper  depths.  This  was  due  to  the  fact  that  105-cm  is  the 
first  depth  that  was  the  bottom  most  layer  for  some  sites.  Since  the  estimated  soil 
parameters  did  not  allow  for  saturation  in  the  lower  profile  layers,  the  soil-water  contents 
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were  largely  underestimated.  For  the  non-bottom  layer,  SAT  is  not  as  critical  because 
saturation  is  only  reached  sporadically  during  the  season,  if  ever. 

The  errors  in  the  yield-based  estimation  of  soil-water  holding  limits  led  to 
unacceptable  errors  in  the  estimation  of  soil-water  content  throughout  the  crop  season, 
especially  for  the  bottom  layer  of  the  soil. 

The  RMSE  of  yield  for  the  grain  yield-based  estimation  procedure  is  good  enough 
for  most  SSM  applications  (less  than  10%  error).  However,  the  fact  that  soil-water 
content  presented  such  high  errors  suggests  that  the  model  might  have  worked  well  for 
the  two  years  of  experimental  data  but  gives  no  confidence  about  the  behavior  in  less 
moderate  conditions  such  as  a wet  year. 


Estimated  SAT  (cm3/cm3)  Estimated  DUL  (cm3/cm3)  Estimated  LL  (cm3/cm3) 
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Figure  50.  Validation  of  the  soil  parameter  estimates  based  on  grain  yield:  LL  estimates 
vs.  Saxton's  estimates  for  1997  (above  left)  and  1998  (above  right);  DUL  estimates  vs. 
Saxton's  estimates  for  1997  (middle  left)  and  1998  (middle  right);  SAT  estimates  vs. 
Saxton's  estimates  for  1997  (below  left)  and  1998  (below  right). 


Estimated  SAT-DUL  (cm 3/cm 3)  Estimated  DUL-LL  (cm 3/cm 3) 
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Saxton  DUL-LL  (cm 3/cm  3) 


Saxton  SAT-DUL  (cm  3/cm  3) 
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Figure  5 1 . Validation  of  the  soil  parameter  estimates  based  on  grain  yield:  DUL  minus 
LL  estimates  vs.  Saxton's  estimates  for  1997  (above  left)  and  1998  (above  right);  SAT 
minus  DUL  estimates  vs.  Saxton's  estimates  for  1997  (below  left)  and  1998  (below  right). 
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Table  16.  Root  mean  square  error  of  soil-water  content  (cm3  cm'3)  and  relative  error  by 
soil  layer  depth  for  the  yield  based  and  soil-water  based  estimation  procedures. 

Yield  based  estimation  Soil-water  based  estimation 


Estimation  set 
Abs.  % 

Validation  set 
Abs.  % 

Estimation  set 
Abs.  % 

Validation  set 
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8.4 
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38.4 

0.018 
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6.7 

150  cm 

0.142 

41.2 
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44.9 

0.022 

6.4 

0.028 

7.5 

30  cm 

0.021 
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0.019 
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7.0 

0.009 
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7.1 
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6.9 
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60  cm 

0.021 

13.5 

0.020 

12.8 

0.011 
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Soil-water  Content  as  Objective  Function  Variable 

Estimating  soil  parameters  using  soil-water  content  as  objective  function  variable 
led  to  larger  ditferences  between  simulated  and  observed  yield  values  in  the  estimation 
data  set  when  compared  to  the  grain  yield-based  estimation.  The  average  RMSEF  for 
both  years  was  442.3  kg/ha  (472.1  kg/ha  for  1997  and  412.5  for  1998,  Table  13).  The 
simulated  versus  observed  yield  graphs  for  the  estimation  data  set  are  shown  in  Figure 
53.  In  1997,  the  model  showed  a tendency  to  under  estimate  higher  yields.  When  soil 
parameters  were  used  to  predict  yield  in  the  independent  year,  the  average  RMSEP  was 
500.1  kg/ha  (555.6  kg/ha  for  1997  and  444.6  kg/ha  for  1998,  Table  13).  The  predicted 
versus  observed  yield  graphs  for  the  validation  set  are  shown  in  Figure  54.  In  1997,  the 
model  showed  a tendency  to  under  estimate  higher  yields  and  in  1998  the  model  showed 
tendency  to  over  estimate  yield  for  the  whole  range  of  observed  values. 

The  resulting  maps  of  simulated  (estimation  data  set)  and  predicted  (validation 
data  set)  grain  yield  for  both  years  are  displayed  in  Figure  55  and  Figure  56,  respectively, 
and  can  be  compared  with  the  observed  yield  map  presented  in  Figure  7.  The  maps  of  the 
difference  between  predicted  and  observed  grain  yield  for  both  years  are  displayed  in 
Figure  57  and  Figure  58.  Only  few  zones  had  RMSE  larger  than  1000  kg/ha  with  the 
majority  of  the  area  with  values  smaller  than  500  kg/ha  for  both  the  estimation  and  the 


validation  data  sets. 
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Figure  52.  Map  of  total  profile  plant  available  soil-water  (mm)  resultant  from  soil 
parameters  estimated  from  soil-water  content  of  1997  (above)  and  1998  (below). 


Simulated  Yield  (kg/ha)  Simulated  Yield  (kg/ha) 
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Figure  53.  Simulated  versus  observed  grain  yield  for  1997  using  soil  parameters 
estimated  from  soil-water  content  of  1 997  (above)  and  for  1 998  using  parameters 
estimated  from  soil-water  content  of  1998. 


Predicted  Yield  (kg/ha) 
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Figure  54.  Predicted  versus  observed  grain  yield  for  1997  using  soil  parameters 
estimated  from  soil-water  content  of  1998  (above)  and  for  1998  using  parameters 
estimated  from  soil-water  content  of  1 997. 
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Figure  55.  Map  of  simulated  com  yield  map  (kg/ha)  for  1997  using  soil  parameters 
estimated  from  soil-water  content  of  1997  (above)  and  for  1998  using  soil  parameters 
estimated  from  soil-water  content  of  1998  (below). 


Figure  56.  Map  of  predicted  com  yield  map  (kg/ha)  for  1997  using  soil  parameters 
estimated  from  soil-water  content  of  1998  (above)  and  for  1998  using  soil  parameters 
estimated  from  soil-water  content  of  1997  (below). 
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Figure  57.  Map  of  the  difference  between  predicted  and  observed  grain  yield  (kg/ha)  for 

1997  using  soil  parameters  estimated  from  soil-water  content  of  1997  (above)  and  for 

1998  using  soil  parameters  estimated  from  soil-water  content  of  1998  (below). 


■■I 


Figure  58.  Map  of  the  difference  between  predicted  and  observed  grain  yield  (kg/ha)  for 

1997  using  soil  parameters  estimated  from  soil-water  content  of  1998  (above)  and  for 

1998  using  soil  parameters  estimated  from  soil-water  content  of  1997  (below). 
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Table  1 1 and  Table  12  display  the  soil  parameter  estimates  for  LL,  DUP  and  SAT 
for  site  2 and  42  resultant  from  the  DSSAT  and  Saxton  pedotransfer  functions  and  from 
the  optimization  procedure.  The  respective  maps  of  profile  plant  available  soil  water  are 
displayed  in  Figure  52. 

Figure  59  and  Figure  60  show  the  values  for  the  optimized  LL,  DUL,  and  SAT 
and  for  the  differences  between  DUL  and  LL  and  between  SAT  and  DUL  against  the 
estimations  of  the  same  quantities  by  the  Saxton’s  pedotransfer  function.  The  estimated 
LL  and  DUL  agreed  well  with  the  expected  values  obtained  from  Saxton’s  pedotransfer 
function.  Similarly  to  the  yield-based  optimization,  SAT  values  continued  to  show  poor 
agreement  with  Saxton’s  pedotransfer  functions.  The  difference  is  that  there  was  an 
exception  of  the  values  for  the  two  layers  closer  to  the  soil  surface  and  the  bottom  layer. 
For  the  rest  of  the  layers,  the  agreement  was  poor.  Plant-available  soil-water  showed  a 
good  agreement  for  the  entire  range  of  optimized  values  whereas  the  difference  between 
SAT  and  DUL  only  showed  reasonable  agreement,  again,  for  the  two  top  layers  and  the 
bottom  layer.  The  agreement  for  the  bottom  layer  makes  sense  accounting  for  the  fact 
that  the  layer  is  saturated  through  the  first  part  of  the  season.  For  the  top  two  layers,  the 
agreement  is  possibly  due  the  temporary  soil-water  contents  above  DUL  after  certain  rain 
events.  The  intermediate  layers  do  not  experience  soil-water  content  above  DUL 
throughout  the  season  and  therefore  are  not  sensitive  to  SAT. 

Table  16  shows  the  root  mean  square  error  and  the  relative  error  of  soil-water 
content  by  soil  layer  depth.  The  magnitude  of  the  error  was  consistent  for  the  two  years 
and  all  the  soil  layers  and  always  smaller  than  10%  of  the  absolute  value  on  average. 
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Figure  59.  Validation  of  the  soil  parameter  estimates  based  on  soil-water  content:  LL 
estimates  vs.  Saxton's  estimates  for  1997  (above  left)  and  1998  (above  right);  DUL 
estimates  vs.  Saxton's  estimates  for  1997  (middle  left)  and  1998  (middle  right);  SAT 
estimates  vs.  Saxton's  estimates  for  1997  (below  left)  and  1998  (below  right). 
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Saxton  DUL-LL  (cm  3/cm  3)  Saxton  DUL-LL  (cm  3/cm 3) 
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Figure  60.  Validation  of  the  soil  parameter  estimates  based  on  soil-water  content:  DUL 
minus  LL  estimates  vs.  Saxton's  estimates  for  1997  (above  left)  and  1998  (above  right); 
SAT  minus  DUL  estimates  vs.  Saxton's  estimates  for  1997  (below  left)  and  1998  (below 
right). 
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The  results  for  the  soil-water-content-based  estimation  are  highly  satisfactory 
since  both  the  soil-water  content  and  yield  were  correctly  simulated  providing  a much 
better  scenario  for  model  use  in  other  conditions  different  than  the  ones  used  in  the  study. 

However,  caution  should  be  taken  in  using  the  model  for  wet  years  since  SAT 
values  were  not  apparently  well  estimated  for  intermediate  soil  layers.  The  estimated 
values  of  SAT  for  intermediate  soil  layers  would  be  certainly  improved  if  experimental 
data  were  available  for  a wet  year. 


Conclusions 

CERES-Maize  ability  to  predict  the  spatial  variability  of  grain  yield  using  site- 
specific  soil  parameters  estimated  from  optimization  based  on  two  distinct  objective 
function  variables  was  studied.  The  objective  function  variables  were  final  grain  yield 
and  the  time  series  of  soil-water  content. 

Estimating  site-specific  soil  parameters  from  grain  yield  led  to  acceptable  errors 
in  grain  simulation  and  prediction.  However,  soil-water  content  was  not  accurately 
predicted  making  the  strategy  unreliable  when  the  goal  is  to  predict  soil-water  content. 
The  errors  in  soil-water  content  were  particularly  high  in  the  bottom  soil  layer.  In 
addition,  most  of  the  soil-water  holding  limits  were  not  valid,  especially  lower  limit  and 
saturation. 

Estimating  site-specific  soil  parameters  from  soil-water  content  measurements  led 
to  acceptable  errors  in  grain  and  soil-water  content  simulation  and  prediction.  The 
estimated  soil-water  holding  limits  were  valid  with  an  exception  of  saturation  for  the 
intermediate  soil  layers.  If  a wet  year  had  made  part  of  the  experimental  set  of  years  then 
saturation  for  this  layer  may  have  been  well  estimated  as  well. 
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The  results  suggest  that  using  time-series  of  soil-water  content  was  most  accurate 
tor  estimating  site-specific  soil  parameters.  This  constitutes  positive  results  in  terms  of 
confidence  in  the  ability  of  crop  model  to  deal  with  spatial  variability  of  water-related 
stresses.  Nevertheless,  in  terms  of  industry  application  of  the  technology  the  results  are 
not  very  encouraging  since  time  series  of  soil-water  content  are  seldom  available. 
Research  on  variations  of  the  grain  yield  based  estimation  procedure  should  then  be 
undertaken.  These  might  include  narrowing  the  range  of  feasible  values  of  each  of  the 
parameters  being  estimated  and/or  the  inclusion  of  more  years  of  observed  data  in  the 
procedure.  Narrowing  the  range  of  feasible  values  can  be  accomplished  by  inference 
techniques  based  on  other  soil  properties  and/or  topography  and/or  remotely  sensed  data. 
Examples  of  promising  inference  techniques  for  this  purpose  are  neural  networks  and 
expert  systems. 


CHAPTER  7 

SUMMARY  AND  CONCLUSIONS 

The  technological  ability  to  apply  agricultural  inputs  to  field  crops  in  a spatially 
variable  and  economically  viable  manner  has  triggered  new  opportunities  in  crop 
management.  The  essence  of  this  new  area  of  endeavor,  designated  as  Site-Specific 
Management  (SSM),  is  the  increase  in  precision  with  which  usual  crop  and  soil 
management  practices  are  performed  in  order  to  apply  optimal  levels  of  inputs  in  all  field 
locations.  This  implies  spatially  variable  nutrient  application,  pesticide  application,  plant 
density,  irrigation,  tillage,  etc,  according  to  specific  characteristics  of  each  point  in  the 
field. 

The  success  of  SSM  depends  on  the  ability  to  relate  the  current  status  of  a crop 
system  and  its  final  physical  performance  (grain  yield  in  most  cases)  for  every  location  of 
a field.  One  of  the  current  challenges  of  SSM  is  the  implementation  of  a tool  for  that 
purpose.  The  desirable  tool  should  (1)  predict  the  spatial  distribution  of  yield  as  early  in 
the  season  as  possible  so  that  possible  repairing  agronomic  actions  could  take  place;  (2) 
use  the  minimal  number  of  variables  as  possible  to  make  it  easy  to  use;  (3)  use  easily 
available  data  to  reduce  cost;  (4)  be  portable  in  space  and  time  to  be  usable  in  new  fields 
and  in  different  years;  (5)  be  sensitive  to  most  important  yield  limiting  factors  to  be  able 
to  explain  variability;  and,  (6)  incorporate  uncertainty  in  inputs  to  estimate  confidence 
limits  for  the  variables  of  interest. 

Procedures  to  predict  the  spatial  pattern  need  to  be  developed  and  validated  over 
space  and  time.  The  focus  on  the  water  limiting  conditions  resulted  from  the  fact  that 
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most  intensively  fanned  fields  are  managed  to  limit  yield  reductions  due  to  nutrients, 
diseases,  pests  and  weeds.  Two  procedures  to  describe  and  predict  the  spatial  pattern  of 
grain  yield  deserved  special  attention:  artificial  neural  networks  and  crop  simulation 
models.  The  two  procedures  are  inherently  different  in  nature.  An  artificial  neural 
network  is  a totally  empirical  tool  whereas  crop  simulation  models  are  based  on  the 
simulation  of  physiological  and  physical  processes.  Prior  to  the  development  of  any 
procedure,  the  mechanisms  that  lead  to  spatial  yield  variation  should  be  understood. 

The  objectives  of  the  dissertation  were  accomplished  in  five  distinct  pieces  of 
research,  each  of  which  was  presented  in  a separate  chapter.  The  first  two  studies 
addressed  the  mechanisms  that  led  to  spatial  yield  variation  by  identifying,  describing 
and  quantifying  the  spatial  variability  of  the  relevant  variables  causing  spatial  yield 
variability.  In  the  last  three  studies  the  procedures  to  predict  the  spatial  pattern  were 
developed  and  validated  over  space  and  time. 

Field  Experiment 

A field  experiment  was  carried  out  for  com  in  a 3.6  ha  farm  field  in  Michigan, 
USA  during  the  crop  seasons  of  1997  and  1998.  The  field  experiment  had  no  imposed 
levels  of  any  kind  of  inputs.  The  experiment  consisted  of  observing  variables  of  interest 
in  a uniformly  managed  field.  The  farmer  managed  the  field  exactly  in  the  same  way  as 
he  had  been  doing  for  over  10  years  (no-tillage,  non-irrigated).  The  field  was  chosen 
from  among  others  for  its  characteristic  spatial  variability  of  final  grain  yield  according  to 
the  farmer.  The  sites  for  observing  soil  and  crop  variables  (total  of  43)  were  spatially 
distributed  in  a grid  of  approximately  30.4  m considered  to  give  a good  representation  of 
all  combinations  of  soil  and  crop  conditions  important  in  the  study. 
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Effective  soil  depth  was  determined  and  soil  particle  analysis  was  performed  for 
each  15 -cm  increment  layer  in  43  all  sites.  Average  clay  and  sand  contents  of  the  top 
seven  layers  (0-105  cm)  were  used  as  direct  indicator  variables  of  the  soil-water  holding 
limits  determining  soil-water  holding  capacity.  In  addition  crop  growth,  soil-water 
content  and  weather  were  monitored  during  the  crop  season.  Leaf  area  index  and  each 
individual  organ  weight  (leaves,  stem  and  grain)  were  determined  three  times  in  the  crop 
season.  Soil-water  content  of  each  15-cm  increment  interval  (down  to  effective  soil 
depth)  for  each  site  was  determined  using  a neutron  moisture  meter.  At  crop  maturity, 
grain  yield  and  its  components  (plant  population,  grain  number  per  ear,  grain  number  per 
unit  area  and  unit  grain  weight)  were  also  determined. 

Study  1 - Relating  Corn  Growth  and  Yield  to  Soil-Water  and  Plant  Population 
It  was  hypothesized  that  plant-available  soil-water  and  plant  population  varied 
considerably  in  space  and  accounted  for  most  of  the  observed  within-field  spatial 
variability  found  in  the  yield  map.  This  hypothesis  was  tested  by  describing  and 
quantifying  the  effects  of  plant  population,  effective  soil  depth,  soil-water  holding  limits 
(soil-water  holding  capacity)  and  seasonal  pattern  of  rainfall  on  the  spatial  pattern  of  corn 
growth,  yield  and  yield  components  using  correlation  analysis  for  simple  effects  and 
forward-stepwise  multiple  regression  analysis  for  multiple  combined  effects.  In  addition, 
the  significance  of  the  effects  of  the  independent  variables,  as  well  as  their  interactions, 
on  the  dependent  variables,  were  tested  using  analysis  of  variance  assuming  an 
unbalanced  fixed  factor  factorial  design.  For  this  purpose,  each  independent  variable  was 
categorized  into  two  classes  based  on  the  median. 
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The  weather  during  the  crop  season  in  1997  was  clearly  drier  than  in  1998.  The 
respective  seasonal  rainfall  amounts  were  226  and  436  mm.  Considering  the  statistics 
from  a daily  weather  data  series  (1979-1988)  of  the  same  location,  1998  was  considered 
an  average  year  in  terms  of  seasonal  rainfall  whereas  1997  was  a dry  one  (the  driest  in  the 
series).  This  difference  made  these  two  years  satisfactory  to  study  the  effects  of  seasonal 
pattern  of  rainfall  on  grain  yield  and  its  spatial  pattern.  The  longest  period  of  rainfall 
absence  in  1 997  followed  planting  and  was  accompanied  by  low  air  temperature  causing 
crop  emergence  to  be  delayed  (it  occurred  30  days  after  planting  versus  only  8 days  in 
1998).  Consequently,  plant  population  at  maturity  was  more  variable  in  1997  than  in 
1998. 

The  average  of  sand  and  clay  content  of  the  top  seven  layers  varied  between  34.3- 
67.7%  and  12.8-20.9%,  respectively.  Effective  soil  depth  varied  between  105  and  150 
cm.  For  the  textures  at  all  sites  and  all  depths  and  according  to  a pedotransfer  function, 
sand  (clay)  was  negatively  (positively)  correlated  with  soil-water  holding  capacity 
varying  between  0.080  and  0.170  cm3  cm  3.  The  total  profile  soil-water  holding  capacity 
varied  between  1 05  and  200  mm. 

The  effective  soil  depth  was  constrained  by  a clay  layer  located  at  variable  depths 
for  each  site.  Typically,  this  layer  causes  a perched  water  table  that  is  present  in  the 
profile  at  the  beginning  of  spring,  due  to  snow  melt,  and  at  the  end  of  the  summer  in  high 
rainfall  years  (summer  storms).  In  the  middle  of  the  season,  the  water  table  disappears 
once  depleted  by  roots.  Therefore,  the  water  table  plays  a very  important  role  in 
supplying  water  during  critical  yield  forming  phases:  flowering  and  grain  filling.  In  years 
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with  higher  rainfall,  the  water  table  is  present  throughout  the  season  at  higher  and  more 
variable  water  levels  causing  excess  water  stress  to  plants. 

The  maximum  plant  population  values  were  approximately  the  same  for  both 
years  (10.3  #/m'2)  whereas  minimum  values  were  much  smaller  for  1997  (5.4  #/m'2  in 
1997  and  7.4  #/m  2 in  1998).  The  spatial  pattern  of  plant  population  showed  little 
similarity  between  years.  The  zones  of  greater  plant  population  reduction  in  1997 
corresponded  to  an  area  of  higher  sand  and  shallower  profiles  suggesting  that  the  delayed 
emergence  was  not  the  only  cause  for  lower  plant  population  but  a combined  effect  with 
low  soil-water  content. 

Grain  yield  was  lower  in  1997  than  in  1998  in  the  entire  field  with  mean  of  9571 
vs.  11331  kg/ha.  The  spatial  pattern  was  similar  between  the  years.  The  difference 
between  years  varied  from  around  3000  kg/ha  in  lower  yielding  areas  to  about  500  kg/ha 
in  higher  yielding  areas.  The  range  of  variability  in  1 997  was  6600  kg/ha  to  12000  kg/ha 
and  9800  kg/ha  to  12600  kg/ha  in  1998.  The  grain  number  per  ear  had  a lower  mean  in 
1997  than  in  1998  (365  vs.  463)  and  its  spatial  pattern  was  very  distinct  for  both  years. 
The  unit  grain  weight  had  a higher  mean  in  1997  than  in  1998  (303.4  vs.  271 .2  Mg)  and 
there  were  similarities  in  the  spatial  pattern  for  both  years.  Grain  number  per  unit  area 
had  a higher  mean  in  1998  than  in  1997  with  no  similarity  between  the  spatial  patterns  of 
variation. 

These  results  emphasized  the  complexity  of  the  yield  forming  process,  its 
interactions  with  the  environment,  and  the  resulting  spatial  yield  variability.  In  1 997,  a 
dry  year,  plant  population  and  plant-available  soil-water  as  affected  by  effective  soil 
depth,  soil-water  holding  limits  and  seasonal  amount  of  rainfall  accounted  for  as  much  as 
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82%  of  the  total  final  grain  yield  spatial  variability.  In  1998,  with  average  seasonal 
rainfall  amounts  and  a more  uniform  stand,  plant-available  soil-water  alone  was  able  to 
explain  68%  of  final  grain  yield  spatial  variability.  The  multiple  linear  regression  models 
used  to  describe  each  year  of  data  performed  poorly  when  used  to  predict  the  grain  yield 
of  the  other  year.  This  result,  again,  stressed  the  peculiarity  of  each  individual  year's  data 
set  in  terms  of  limiting  factors  and  the  inability  of  these  multiple  linear  regression  models 
to  account  for  the  complex  interactions  present  in  a yield  formation  process. 

Plant-available  soil-water  as  affected  by  effective  soil  depth,  soil-water  holding 
limits  and  seasonal  amount  of  rainfall,  was  the  single  most  important  factor  responsible 
for  the  spatial  pattern  of  corn  growth,  yield  and  yield  components  in  the  studied  non- 
irrigated  field.  In  both  situations  effective  soil  depth  was  more  important  than  sand 
content,  or  specific  plant-available  soil-water,  in  causing  spatial  variability  of  grain  yield. 
Unit  grain  weight  was  more  important  in  defining  the  spatial  pattern  of  yield  variability 
whereas  the  number  of  grains  per  unit  area  was  more  important  in  defining  the  magnitude 
of  total  field  yield. 

Study  2 - Relating  Corn  Yield  and  Soil  Properties  to  Topography 

In  an  effort  to  characterize  and  understand  existing  spatial  variability  of  crop  and 
soil  variables  we  looked  at  topography.  It  was  hypothesized  that  the  spatial  pattern  of 
soil  properties  related  to  plant-available  soil-water  and  crop  growth,  yield  and  yield 
components  of  corn  grown  under  non-irrigated  conditions  were  related  to  certain 
topographic  terrain  attributes.  In  addition,  we  aimed  to  understand  how  topography 
affects  crop  growth  and  yield  for  this  particular  field.  It  was  hypothesized  that  the  effects 
could  be  through  the  short-term  effect  on  soil-water  content  (lateral  water  movement)  or 
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through  its  long-term  (geological  time)  effects  on  site-specific  soil  properties  that  in  turn 
determine  differences  in  plant-available  soil-water  patterns.  This  issue  was  considered 
relevant  because  of  its  implication  in  terms  of  modeling  efforts. 

These  hypotheses  were  tested  by  describing  and  quantifying  the  effects  of 
topography  on  effective  soil  depth,  clay  and  sand  contents,  corn  growth  variables,  and 
yield  and  yield  components  using  correlation  analysis  for  simple  effects  and  forward- 
stepwise  multiple  regression  analysis  for  multiple  combined  effects.  The  rationale  for 
understanding  how  topography  affected  crop  growth  and  yield  for  this  particular  field 
was  also  based  on  a correlation  analysis. 

The  topography  of  the  field  was  characterized  by  the  following  topographic 
terrain  attribute,  altitude,  slope,  plan  curvature,  profile  curvature,  specific  catchment  area, 
distance  to  network,  downstream  and  upstream  flow  length,  wetness  index  and  stream 
power  index. 

The  elevation  relief  of  the  field  was  2.65  m with  the  gradient  well  distributed  in 
the  area  with  an  average  slope  of  1 .3%.  Only  three  small  areas  (about  0.7%  of  the  field 
area)  had  slopes  higher  than  4.5%  that  corresponded  to  field  outlet  locations.  The  field 
had  four  drainage  basins  bigger  than  600  m2  the  largest  with  an  area  of  1 .25  ha.  Due  to 
the  predominantly  low  slopes,  specific  catchment  area  dominated  the  calculation  of  both 
wetness  index  and  stream  power  index. 

Field  areas  with  higher  elevation,  farther  from  the  stream  network,  with  lengthier 
downstream  flow  and  smaller  wetness  index  had  higher  (lower)  sand  (clay)  content  and 
vice-versa.  Field  areas  with  lower  elevation,  shorter  downstream  flow,  lower  slope  and 
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closer  to  the  flow  network  had  higher  effective  soil  depths.  The  strongest  relation  in  the 
study  was  between  effective  soil  depth  and  downstream  flow  length. 

The  number  of  significant  relationships  between  topographic  attributes  and  crop 
variables  in  1997  was  higher  and  they  were  larger  in  value  when  compared  to  1998.  This 
fact  made  sense  considering  that  1997  was  a drier  year  than  1998,  and  therefore  plant- 
available  soil-water  was  more  limiting  to  crop  growth  in  1997.  The  topographic 
attributes  that  were  more  significantly  correlated  with  crop  variables  were  the  same  for 
both  years:  Elevation,  distance  to  network,  wetness  index,  plan  curvature  and 
downstream  flow  length.  Furthermore,  these  were  the  topographic  attributes  that  showed 
the  highest  correlation  with  soil  properties.  Downstream  flow  length  was  the  topographic 
attribute  that  had  the  highest  correlation  with  yield  and  most  yield  components.  This 
attribute  was  also  more  highly  correlated  to  soil  variables,  especially  soil  depth. 

These  observations  were  consistent  with  the  relationships  described  previously 
between  the  soil  properties  and  yield  and  yield  components.  Effective  soil  depth  and  clay 
and  sand  contents  affected  yield  formation  and  were  affected  by  topography  (especially 
related  to  altitude,  distance  to  network  and  downstream  flow  length).  These  relationships 
were  affected  by  the  interactions  with  weather,  in  particular,  rainfall. 

The  results  seemed  to  show  that  for  dry  and  normal  rainfall  years,  the  effect  of 
topography  on  the  spatial  variability  of  corn  growth  and  yield,  in  this  particular 
landscape,  was  relatively  more  important  through  the  long  term  effect  (geological  time) 
on  soil  properties  (particularly  effective  soil  depth)  than  the  short  term  effect  on  seasonal 
plant-available  soil-water.  However,  in  wetter  years,  the  relative  importance  of  the  two 
effects  might  change  since  the  magnitude  of  the  seasonal  lateral  surface  and  sub-surface 
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soil-water  movement  changes.  Topographic  variables  were  clearly  shown  to  have  a great 
potential  in  describing  the  spatial  variability  of  soil  properties  and.  therefore,  should  be 
considered  in  efforts  to  predict  crop  performance. 

Study  3 - Prediction  of  Spatial  Yield  Variability  With  Neural  Network  Models 
In  an  effort  to  develop  a tool  to  predict  the  spatial  pattern  of  yield  we  first 
concentrated  on  neural  network  models.  Neural  network  models  are  especially  well 
suited  for  complex  systems  where  the  underlying  processes  are  not  well  understood  or  are 
difficult  to  translate  into  mathematical  equations  and  when  a myriad  of  variables  are 
assumed  to  have  an  effect  on  the  variable  of  interest.  In  this  empirical  approach  the 
inclusion  of  new  variables  is  simplified.  Three  distinct  neural  network  models  were 
trained  and  tested:  using  only  topographic  variables,  using  only  agronomic  variables  and 
using  both  types  of  variables. 

All  three  neural  networks  models  had  a similar  structure:  feedforward 
architecture,  backpropagation  weight  adjustment  algorithm,  one  hidden  layer  with 
sigmoid  transfer  function  and  linear  transfer  functions  for  input  and  output  layers.  The 
major  difference  between  the  neural  network  models  were  the  number  of  nodes  in  the 
hidden  layer. 

Since  the  training  and  validation  sets  were  obtained  from  a two-year  experiment 
in  one  field,  the  data  set  had  to  be  split.  The  field  was  divided  into  two  sections  along  the 
transversal  direction.  The  training  and  testing  data  sets  were  evaluated  through  the  root 
mean  square  error  of  fit  and  of  prediction  among  sites,  respectively.  A field-wide 
evaluation  was  carried  out  through  visual  inspection  of  the  spatial  pattern  of  yield 
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variation  obtained  through  kriging  and  also  by  a local  measure  of  the  root  mean  square 
error. 

As  expected,  the  root  mean  square  error  of  fit  for  all  training  sets  were  smaller 
than  the  ones  for  testing  sets.  The  topographic  attributes  neural  network  had  the  worst 
predictive  ability  (RMSEP  of  1264  kg/ha)  followed  by  the  agronomic  variables  network 
(RMSEP  of  1 034  kg/ha).  The  neural  network  that  account  for  all  the  variables  in  the 
study  resulted  in  the  lowest  RMSEP  (9 1 8 kg/ha).  The  gain  in  predictability  by  using  all 
the  variables  did  not  increase  linearly.  This  was  possibly  due  to  the  correlation  between 
variables. 

Neural  network  models  allowed  for  new  sources  of  data  (topography)  to  be  used 
in  predicting  crop  performance.  The  model  that  included  all  soil,  crop  and  topographic 
variables  presented  a reasonable  predictability  of  the  spatial  variability  of  final  grain  yield 
for  practical  applications.  However,  one  serious  limitation  of  this  approach  was  the 
limited  gain  in  the  meaningful  system  relationships  that  could  be  translated  into  'real 
system  knowledge'  and  used  to  predict  the  spatial  variability  of  final  grain  yield  in  other 
environments  (different  weather,  soils,  fields,  etc).  Nevertheless,  if  the  future  of  the 
application  of  prediction  tools  in  Site-Specific  Management  is  to  involve  specific  field 
calibrations,  neural  network  models  are  promising.  Their  role  can  either  be  as 
independent  tools  in  predicting  grain  yield  from  environmental  variables  or  as  a support 
to  other,  more  mechanistic  simulation  tools,  in  providing  initial  conditions  values  or  site- 
specific  soil  parameters. 
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Study  4 - Prediction  of  Spatial  Yield  Variability  With  Crop  Simulation  Models 
In  an  effort  to  develop  a tool  to  predict  the  spatial  pattern  of  yield  the  second 
procedure  that  deserved  our  attention  was  crop  simulation  modeling.  Crop  simulation 
models  are  process-based  algorithms  that  attempt  to  predict  final  grain  yield  by 
describing  the  most  important  physical/chemical/biological  processes  responsible  for 
yield  formation.  Since  crop  simulation  models  provide  a comprehensive  approach  to 
predicting  crop  growth  and  development  and  their  environmental  determinant  factors  and 
interactions,  they  are  especially  suited  for  testing  management  prescription  and  being 
able  to  be  used  in  a wide  range  of  environment. 

CERES-Maize  V 3.5,  a worldwide  used  corn  simulation  model,  was  used  to 
predict  the  spatial  variability  of  corn  growth  and  yield  and  components.  It  was  also  used 
to  study  the  relative  contribution  of  each  spatially  variable  input  to  the  yield  variability 
and  its  predictability.  Based  on  previous  studies,  it  was  hypothesized  that  factors 
affecting  plant-available  soil-water  (soil-water  holding  limits  and  effective  soil  depth) 
and  plant  population  were  the  major  causes  for  spatial  variability  of  crop  growth  and 
yield.  Accordingly  the  model  was  run  for  each  site  varying  plant  population,  soil  depth, 
lower  limit,  drained  upper  limit,  saturation,  root  distribution  factor  and  saturated 
hydraulic  conductivity  of  the  bottom  most  soil  layer. 

In  order  to  test  the  ability  of  the  model  to  predict  the  spatial  variability  of  corn 
growth  and  yield  it  was  run  for  both  years  with  the  most  accurate  inputs  for  each  site. 
Predictions  made  of  corn  growth  and  yield  were  thus  independent  of  observed  values  in 
this  field  for  both  years.  The  agreement  between  predicted  and  observed  crop  growth, 
grain  yield  and  components  was  evaluated  by  the  root  mean  square  error  of  prediction 
among  sites.  A field-wide  evaluation  was  carried  out  through  visual  inspection  of  the 
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spatial  pattern  ot  yield  variation  obtained  through  kriging  and  also  by  a local  measure  of 
the  root  mean  square  error.  To  study  the  relative  contribution  of  each  spatially  variable 
input  the  model  was  run  for  both  years  with  site-specific  plant-available  soil-water  and 
homogenous  plant  population  (average  among  sites). 

The  average  root  mean  square  error  of  prediction  of  grain  yield  for  both  years  was 
501  kg/ha,  which  is  an  exceptionally  low  value.  The  proportions  of  yield  variability 
explained  by  the  model  were  88%  and  76%  for  1997  and  1998,  respectively.  The  better 
performance  of  the  model  in  1997  compared  to  1998  made  sense  considering  that  1997 
was  a dry  year.  The  model  also  performed  very  well  in  predicting  most  of  the  yield 
components  and  growth  variables  as  well.  The  model  only  performed  poorly  in 
predicting  LAI  at  anthesis  for  1998.  The  fact  that  the  model  was  able  to  predict  grain 
yield  as  well  as  yield  components  and  growth  in  a very  satisfactory  way  increased  the 
confidence  in  its  performance  since  it  is  a mechanistic  tool.  On  average,  93%  of  the  yield 
spatial  variability  explained  by  the  model  (average  of  82%)  was  due  to  spatially  variable 
plant-available  soil-water  and  7%  is  due  to  spatially  variable  plant  population. 

CERES-Maize  V3.5  accurately  predicted  the  spatial  variability  of  corn  growth 
and  yield  tor  two  years  of  independent  data  when  accurate  inputs  were  provided  for  soil 
properties  and  plant  population.  These  excellent  results  indicate  that  the  model  can  be 
used  for  further  SSM  studies  when  inputs  are  accurately  known.  Nevertheless,  prudence 
is  advised  in  the  generalization  of  such  results  because  of  the  specificity  of  the  field 
studied  situation:  intensively  observed  field  experiment  with  no  stresses  other  than  ones 
related  to  soil-water.  In  addition,  the  two  years  of  experiment  did  not  include  a wet  year 
where  excess  water  stress  plays  an  important  in  shaping  final  grain  yield  spatial  pattern. 
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In  more  common  situations  the  spatial  pattern  of  grain  yield  is  an  integration  of  multiple 
interacting  dynamic  stresses  caused  by  nitrogen  deficiency,  high  weed  densities,  disease 
incidence,  etc.  The  wide  use  of  crop  simulation  models  in  SSM  is  not  compatible  with 
the  sampling  intensity  carried  out  in  this  field  to  obtain  inputs  and  parameters  essential 
for  the  accuracy  of  the  predictions.  Methods  to  obtain  site-specific  model  inputs  and 
parameters  other  than  direct  measurement  have  to  be  investigated  in  order  to  allow  the 
generalization  of  crop  simulation  models  use  in  SSM. 

Study  5 - Estimation  of  Site-Specific  Crop  Model  Soil  Inputs 

Finally,  a method  to  obtain  site-specific  soil  inputs  for  crop  models  (CERES- 
Maize  V3.5  in  this  case)  was  studied.  Despite  their  accuracy  and  value,  the 
generalization  of  crop  simulation  models  use  in  SSM  can  become  impractical  by  the  need 
for  accurate  values  of  inputs  for  numerous  sites  in  a field.  Some  model  inputs,  such  as 
soil-water  holding  limits  of  each  soil  layer  or  saturated  hydraulic  conductivity  of  any 
drainage  restrictive  layer,  are  difficult  or  expensive  to  determine,  but  are  of  crucial 
importance  to  model  performance  in  non-irrigated  conditions.  It  is,  therefore,  extremely 
important  that  we  measure  or  estimate  them  correctly.  Laboratory  and  field  procedures 
are  not  as  inexpensive  and  practical  as  needed  for  site-specific  management.  Another 
way  to  obtain  model  inputs  involves  optimization  for  the  estimation  of  parameters  using 
more  easily  measured  variables. 

In  the  optimization  approach,  a framework  is  established  where  all  model  inputs, 
with  the  exception  of  the  ones  to  be  estimated,  are  considered  to  be  known.  It  is  also 
assumed  that  the  present  model  accurately  simulates  yield  when  inputs  are  known  and 
accurate.  With  these  assumptions,  the  lack  of  fit  between  observed  and  simulated  data  is 
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attributed  to  non-optimal  estimation  of  the  spatially  variable  inputs.  Subsequently,  new 
input  estimations  are  obtained  in  an  attempt  to  improve  the  fit. 

Once  site-specific  crop  model  parameters  are  obtained  in  this  way  for  a specific 
field  and  validated,  they  can  be  used  to  run  the  crop  model  to  test  management 
prescriptions  and  to  correctly  assess  their  agronomic  and  economic  outcomes.  Such 
procedures  could  be  used  by  industry  in  their  efforts  to  use  crop  simulation  models  in 
SSM  programs  because  databases  of  multiple-year  data  (yield  map,  weather, 
management,  etc)  from  specific  customers'  fields  are  being  built.  Each  new  year  of  data 
added  to  the  database  at  the  end  of  the  crop  season  could  be  used  to  increase  the 
robustness  of  the  optimized  crop  model's  inputs,  further  improving  the  accuracy  of  crop 
model's  prediction. 

One  serious  question  exists  regarding  the  validity  of  site-specific  parameter 
estimation  concerning  the  model  variable  used  to  judge  the  goodness  of  fit.  Final  grain 
yield  may  be  the  ultimately  variable  of  interest,  but  the  objective  function  of  the  search 
technique  should,  as  much  as  possible,  use  a model  variable(s)  that  is  (are)  directly 
related  to  the  parameters  being  estimating. 

The  ability  of  a crop  simulation  model  to  predict  the  spatial  variability  of  grain 
yield  using  site-specific  soil  parameters  estimated  from  optimization  based  on  two 
distinct  objective  function  variables  was  studied.  The  objective  function  variables  were 
the  final  grain  yield  and  the  time  series  of  soil-water  content.  The  comparison  aimed  to 
answer  the  following  questions:  Do  estimated  soil-water  holding  limits  based  on  grain 
yield  as  objective  function  variable  provide  an  accurate  simulation  of  soil-water  content? 
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Do  estimated  soil-water  holding  limits  based  on  time  series  of  soil-water  content  as 
objective  function  variable  provide  good  simulation  of  grain  yield? 

A negative  answer  to  the  second  question  independently  of  the  answer  to  the  first 
question  means  that  the  crop  model  should  be  revised.  The  optimal  scenario  is  a positive 
answer  to  both  questions  because  it  implies  that,  if  input  estimations  are  valid,  an  easily 
measured  variable  such  as  grain  yield  could  be  used  to  estimate  site-specific  soil 
parameters  for  crop  simulation.  If,  however,  only  the  second  question  has  a positive 
answer,  and  input  estimations  are  valid,  then  only  the  time  series  of  soil-water  content 
could  be  used  to  estimate  soil-water  holding  limits  and  the  extensive  application  of  the 
parameters  estimation  procedure  is  endangered. 

The  choice  of  soil  parameters  to  estimate  for  each  particular  site  in  the  field 
(lower  limit,  drained  upper  limit,  saturation,  root  distribution  factor  and  saturated 
hydraulic  conductivity  of  the  bottom-most  soil  layer)  was  based  on  previous  studies 
showing  plant-available  soil-water  to  be  the  main  factor  responsible  for  spatial  variability 
of  yield.  The  chosen  crop  model  soil  parameters  were  the  ones  that  mainly  affect  plant- 
available  soil-water  (with  exception  of  effective  soil  depth  considered  to  be  known). 

A simulated  annealing  algorithm  was  used  as  the  parameter  search  technique  to 
minimize  the  objective  function  for  each  of  the  variables  of  interest.  The  objective 
function  was  the  root  mean  square  error  when  the  variable  of  interest  was  either  final 
yield  or  time-series  of  soil-water  content.  Error  was  calculated  as  the  arithmetic 
difference  between  observed  and  simulated  values. 

The  ranges  of  search  for  the  soil-water  holding  limits  were  based  on  a 
pedotransfer  function.  Since  simulated  annealing  is  a global  search  technique  there  was 
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no  need  to  use  multiple  initial  values  for  the  parameters.  The  validation  of  the  soil 
parameters  generated  by  the  site-specific  input  estimation  procedure  was  also  done  by 
comparing  the  estimates  generated  by  the  procedure  with  the  ones  estimated  from  a 
pedotransfer  function. 

For  both  objective  function  variables,  1997  and  1998  experimental  years  were 
used  as  both  estimation  and  validation  data  sets  alternately.  The  estimation  and 
validation  data  sets  were  evaluated  by  the  root  mean  square  error  of  fit  and  prediction  of 
grain  yield  among  sites.  A field-wide  evaluation  was  carried  out  through  visual 
inspection  of  the  spatial  pattern  of  yield  variation  obtained  through  kriging  and  also  by  a 
local  measure  of  the  root  mean  square  error. 

Estimating  soil  parameters  using  grain  yield  as  the  objective  function  variable  led 
to  minimal  differences  between  simulated  and  observed  yield  values  in  the  estimation 
data  set.  When  these  soil  parameters  were  used  to  predict  yield  in  the  independent  year 
(validation  data  set),  the  average  RMSEP  was  765  kg/ha.  The  estimated  lower  limit 
tended  to  be  consistently  below  the  expected  values  obtained  from  the  pedotransfer 
function  for  lower  limits  lower  than  0.100  cm3/cm3.  Saturation  was  below  the  expected 
value  for  the  whole  range  of  values.  On  the  contrary,  drained  upper  limit  agreed  well 
with  the  pedotransfer  estimates.  As  a result,  the  grain  yield  based  estimation  resulted  in 
larger  plant-available  soil-water  values  for  sites  with  lower  lower  limit.  The  magnitude 
of  the  error  for  soil-water  content  was  consistent  for  the  two  years  and  varied  according 
to  layer.  Down  to  90  cm,  the  average  relative  error  was  always  below  20%  and  increased 
to  45%  for  deeper  depths.  These  errors  were  considered  too  high  even  though  grain  yield 
was  acceptably  well  predicted.  However,  the  fact  that  soil-water  content  simulations 
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presented  such  high  errors  suggests  that  the  model  might  have  worked  well  for  the  two 
years  of  experimental  data  but  gives  no  confidence  about  the  behavior  in  less  moderate 
conditions  such  as  a wet  year. 

Estimating  soil  parameters  using  soil-water  content  as  objective  function  variable 
led  to  larger  differences  between  simulated  and  observed  yield  values  in  the  estimation 
data  set  when  compared  to  the  grain  yield-based  estimation.  The  average  root  mean 
square  error  ol  fit  for  both  years  was  442  kg/ha.  When  the  soil  parameters  were  used  to 
predict  yield  in  the  independent  year,  the  average  RMSEP  was  500  kg/ha.  The  estimated 
lower  limit  and  drained  upper  limit  agreed  well  with  the  expected  values  obtained  from 
the  pedotransfer  function.  Similarly  to  the  yield-based  optimization,  saturation  values 
continued  to  show  bad  agreement  with  pedotransfer  functions.  The  difference  was  that 
there  was  an  exception  of  the  values  for  the  two  layers  closer  to  the  soil  surface  and  the 
bottom  layer.  For  the  rest  of  the  layers,  the  agreement  was  poor.  Plant-available  soil- 
water  showed  a good  agreement  for  the  all  range  of  optimized  values.  The  agreement  for 
the  bottom  layer  makes  sense  accounting  for  the  fact  that  the  layer  is  saturated  through 
the  first  part  of  the  season.  For  the  top  two  layers,  the  agreement  is  possibly  due  the 
temporary  soil-water  contents  above  drained  upper  limit  after  certain  rain  events.  The 
intermediate  layers  do  not  experience  soil-water  content  above  drained  upper  limit 
throughout  the  season  and  therefore  are  not  sensitive  to  saturation.  The  magnitude  of  the 
error  for  soil-water  content  was  consistent  for  the  two  years  and  all  the  soil  layers  and 
always  smaller  than  10%  of  the  absolute  value  on  average. 

The  results  for  the  soil-water-content-based  estimation  were  highly  satisfactory 
since  both  the  soil-water  content  and  yield  were  correctly  simulated  providing  a much 
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better  scenario  tor  model  use  in  other  conditions  different  than  the  ones  used  for  in  the 
study.  However,  caution  should  be  used  in  applying  the  model  for  wet  years  since 
saturation  values  were  apparently  not  well  estimated  for  intermediate  soil  layers.  The 
estimated  values  of  saturation  for  intermediate  soil  layers  would  be  certainly  improved  if 
experimental  data  were  available  for  a wet  year. 

Estimating  site-specific  soil  parameters  from  grain  yield  led  to  acceptable  errors 
in  grain  simulation  and  prediction.  However,  soil-water  content  was  not  accurately 
predicted  making  the  strategy  less  reliable  when  the  goal  is  to  predict  soil-water  content. 
Estimating  site-specific  soil  parameters  from  soil-water  content  measurements  led  to 
acceptable  errors  in  grain  and  soil-water  content  simulation  and  prediction. 

These  results  suggest  that  using  time-series  of  soil-water  content  was  most 
accurate  for  estimating  site-specific  soil  parameters.  This  constitutes  positive  results  in 
terms  ot  confidence  in  the  ability  of  crop  model  to  deal  with  spatial  variability  of  water- 
related  stresses.  Nevertheless,  in  term  of  industry  application  of  the  technology  the 
results  are  not  very  encouraging  since  time  series  of  soil-water  content  are  seldom 
available.  Research  on  variations  of  the  grain  yield  based  estimation  procedure  should 
then  be  undertaken.  These  might  include  the  narrowing  the  range  of  feasible  values  of 
each  of  the  parameters  being  estimated  and/or  the  inclusion  of  more  years  of  observed 
data  in  the  procedure.  Narrowing  the  range  of  feasible  values  can  be  accomplished  by 
inference  techniques  based  on  other  soil  properties  and/or  topography  and/or  remotely 
sensed  data.  Examples  of  promising  inference  techniques  for  this  purpose  are  neural 
networks  and  expert  systems. 
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General  Conclusions.  Recommendations  and  Future  Research 
Site-specific  management  requires  precise  site-specific  rate  recommendations  of 
the  key  agricultural  inputs.  The  development  of  such  recommendations  involves  the 
consideration  of  multiple  variables  such  as  soil  characteristics,  crop  status,  cultivar 
specific  characteristics  and  weather.  Traditional  recommendation  procedures  are,  by 
definition,  not  applicable  to  site-specific  management.  Due  to  the  complexity  of  the 
interactions  among  critical  variables  important  to  site-specific  recommendations,  a 
system  is  needed  to  help  in  the  decision  making  process.  A key  component  in  such 
system  is  a procedure  that  inter-relates  the  variables  determining  final  economic  yield. 

We  have  tested  two  procedures  for  such  proposes  in  a soil-water  limiting 
environment:  Neural  network  models  and  crop  simulation  models.  Crop  simulation 
models  clearly  outperformed  neural  networks  models.  Crop  simulation  models  offered 
better  predictability  of  the  physical  performance  of  the  systems  as  well  a comprehensive 
approach  to  the  site-specific  limiting  factors.  Therefore,  this  would  be  the  indicated 
procedure  to  integrate  in  a decision  support  system  for  site-specific  management. 

The  success  of  crop  simulation  models  in  helping  decision  making  for  site- 
specific  management  will  obviously  depend  on  their  prediction  accuracy.  The  two  main 
factors  that  should  deserve  especial  attention  are,  then,  the  estimation  of  site-specific  crop 
model  inputs  and  parameters  (especially  for  soil  and  management)  and  the 
appropriateness  of  the  equations  used  in  the  crop  simulation  models.  Therefore,  future 
research  efforts  should  focus  on:  (1)  developing  estimation  procedures  of  site  specific 
crop  model  inputs  and  parameters  and  (2)  further  improving  the  model  equations  to 
adequately  account  for  all  yield  limiting  factors. 


APPENDIX  A 
CROP  AND  SOIL  DATA 

The  observed  values  of  all  variables  related  to  soil  and  crop  growth  and  yield  are 
displayed. 
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Table  17.  Observed  site  values  of  crop  variables  in  1997. 

Site 

Yield  (Kg/ha) 

Plant  population 
(#/m2) 

Unit  grain  weight 

(mg) 

Number  of  grain 
per  unit  area 
(#/m2) 

Number  of  grain 
per  ear 

l 

8877 

7.1 

266 

2826 

398 

11671 

7.7 

324 

2934 

381 

11963 

8.8 

320 

3335 

379 

5 

11415 

9.9 

322 

3594 

363 

9458 

9.1 

300 

3067 

337 

11946 

10.1 

314 

3151 

312 

7 

10946 

8.8 

327 

3054 

347 

11423 

9.7 

320 

3017 

311 

9 

10 
11 

10114 

9.8 

320 

2685 

274 

10480 

9.8 

306 

3342 

341 

10123 

8.1 

326 

2770 

342 

12 

13 

14 

9547 

9.5 

305 

3249 

342 

10431 

9.6 

326 

3158 

329 

10512 

8.0 

325 

3024 

378 

15  _ . 

16 

17 

8748 

8.1 

318 

2673 

330 

11345 

10.3 

321 

3214 

312 

10543 

9.8 

324 

2783 

284 

18 

9716 

8.6 

326 

2614 

304 

19 

11345 

10.1 

320 

2919 

289 

20 

8282 

7.5 

290 

3188 

425 

21 

9246 

7.8 

310 

2980 

382 

22 

9189 

8.2 

310 

2936 

358 

23 

10548 

8.9 

312 

3231 

363 

24 

6620 

5.9 

269 

2826 

479 

25 

8735 

7.2 

292 

3305 

459 

26 

8937 

7.4 

291 

2908 

393 

27 

10632 

9.7 

312 

3317 

342 

28 

9774 

9.5 

329 

3031 

319 

29 

8747 

8.9 

281 

2821 

317 

30 

9904 

7.9 

304 

2978 

377 

31 

8986 

8.2 

276 

3296 

463 

32 

8368 

7.5 

288 

2993 

399 

33 

8347 

5.8 

276 

2906 

501 

34 

8764 

8.5 

276 

2788 

328 

35 

8286 

8.4 

297 

2982 

355 

36 

7345 

5.7 

267 

2759 

484 

37 

8647 

7.1 

281 

2840 

400 

38 

8579 

7.6 

285 

3101 

408 

39 

8062 

8.1 

290 

2924 

361 

40 

9560 

8.8 

326 

2640 

300 

41 

8756 

8.2 

311 

3009 

367 

42 

6800 

7.3 

267 

2949 

404 

43 

9843 

8.9 

295 

3213 

361 

150 


Table  1 8.  Observed  site  values  of  crop  variables  in  1998 

Site 

Yield  (Kg/ha) 

Plant  population 
(#/m2) 

Unit  grain  weight 
(mg) 

Number  of  grain 
per  unit  area 
(#/m2) 

Number  of  grain 
per  ear 

1 

10795 

9.2 

256 

3882 

422 

12366 

8.9 

282 

4183 

470 

4 

5 

12629 

8.7 

283 

4376 

503 

12531 

9.7 

276 

4220 

435 

11593 

10.2 

282 

4406 

432 

12626 

9.6 

281 

4522 

471 

8 

11114 

8.9 

278 

3631 

408 

12140 

9.8 

277 

4253 

434 

11433 

9.9 

274 

3980 

402 

11605 

8.5 

276 

4174 

491 

12205 

8.9 

289 

4468 

502 

2 

10334 

8.1 

260 

4066 

502 

13  

14 

15 

11982 

8.8 

288 

4611 

524 

11170 

9.7 

279 

4084 

421 

11433 

8.3 

266 

4017 

484 

16 

17 

12094 

9.7 

281 

4443 

458 

11037 

8.3 

277 

4001 

482 

18 

11763 

8.9 

277 

4183 

470 

19 

11798 

10.3 

274 

4614 

448 

20 

10216 

9.0 

264 

4320 

480 

21 

10939 

9.4 

283 

4061 

432 

22 

12152 

8.3 

277 

4548 

548 

23 

12155 

8.3 

270 

4590 

553 

24 

10231 

8.6 

239 

4223 

491 

25 

10829 

9.4 

264 

3760 

400 

26 

10871 

8.6 

261 

3784 

440 

27 

12416 

9.1 

288 

4386 

482 

28 

12242 

9.2 

285 

4637 

504 

29 

10961 

10.2 

274 

3947 

387 

30 

11145 

9.8 

277 

4018 

410 

31 

11396 

9.0 

264 

4698 

522 

32 

10639 

8.2 

263 

3961 

483 

33 

9940 

7.4 

258 

3848 

520 

34 

11154 

8.6 

268 

3861 

449 

35 

10769 

8.7 

268 

4228 

486 

36 

10518 

8.2 

264 

4067 

496 

37 

10312 

8.9 

254 

4058 

456 

38 

10686 

9.6 

255 

3715 

387 

39 

11089 

9.6 

263 

4090 

426 

40 

11520 

9.8 

284 

4096 

418 

41 

11562 

8.7 

273 

4289 

493 

42 

9834 

9.3 

238 

4129 

444 

43 

11045 

8.8 

270 

3995 

454 

151 


Table  19.  Observed  site  values  of  soil  variables. 


Site 

Soil  depth 
(cm) 

Clay  Content  (%) 
0-105  cm 

Sand  Content  (%) 
0-105  cm 

1 

120 

15.8 

53.4 

2 

3 

150 

19.9 

37.2 

150 

20.3 

36.7 

150 

18.8 

37.8 

D 

135 

19.5 

43.1 

150 

20.3 

36.3 

8 

150 

18.0 

45.7 

150 

20.8 

36.7 

150 

18.8 

46.2 

10 

135 

17.3 

48.1 

11 

135 

18.1 

43.4 

12 

120 

17.1 

48.0 

13 

14 

135 

20.7 

38.6 

150 

18.3 

43.8 

15 

16 

135 

17.4 

45.2 

150 

18.1 

43.0 

17 

18 

150 

16.6 

44.4 

135 

18.5 

43.2 

19 

150 

19.7 

37.9 

20 

120 

17.9 

44.0 

21 

135 

16.4 

51.4 

22 

135 

20.9 

38.1 

23 

135 

20.9 

38.5 

24 

105 

14.7 

61.0 

25 

120 

18.5 

42.1 

26 

120 

17.4 

45.6 

27 

135 

20.5 

34.3 

28 

135 

19.1 

36.0 

29 

120 

19.0 

41.7 

30 

135 

18.0 

41.7 

31 

120 

16.9 

46.9 

32 

120 

15.0 

60.2 

33 

120 

15.7 

51.8 

34 

120 

15.6 

53.3 

35 

120 

14.3 

57.3 

36 

120 

13.4 

62.1 

37 

120 

16.7 

54.9 

38 

120 

16.7 

49.4 

39 

120 

14.8 

56.8 

40 

135 

16.2 

45.5 

41 

135 

15.4 

52.7 

42 

105 

12.8 

67.7 

43 

120 

16.7 

45.9 

152 


APPENDIX  B 

ADDITIONAL  GRAPHS  FOR  CHAPTERS  5 AND  6 
The  intent  of  this  appendix  is  to  display  graphs  of  the  time  course  of  observed  and 
simulated  soil  water  contents  for  the  different  estimation  procedures.  Due  to  the  extent  of 
the  dataset  used,  only  data  from  two  sites  are  displayed.  We  have  chosen  the  sites 
according  to  the  value  of  plant  available  soil  water.  These  sites  (2  and  42)  had  the  lowest 
and  the  highest  values  of  estimated  plant  available  soil  water,  respectively. 
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Figure  61 . Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  30  cm  depth  in  site  2 in  1997. 
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Figure  62.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  45  cm  depth  in  site  2 in  1997. 
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Figure  63.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  60  cm  depth  in  site  2 in  1997. 
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Figure  64.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  75  cm  depth  in  site  2 in  1997. 


SWC  (cm3/cm3) 


155 


0.500 
0.450 
0.400 
o'  0.350 


0.000 


y97 

W98 

- - - - W97 


150  165  180  195  210 


225 

DOY 


240  255  270  285  300 


Figure  65.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  90  cm  depth  in  site  2 in  1997. 


0.500 

0.450 

0.400 

0.350 

0.300 


0.000 


y97 

W98 

- - - -W97 
• OBS 


150  165  180  195  210 


225 

DOY 


240  255  270  285  300 


Figure  66.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  105  cm  depth  in  site  2 in  1997. 
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Figure  67.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  120  cm  depth  in  site  2 in  1997. 
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Figure  68.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  135  cm  depth  in  site  2 in  1997. 
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Figure  69.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  150  cm  depth  in  site  2 in  1997. 
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Figure  70.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  30  cm  depth  in  site  42  in  1997. 
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Figure  71 . Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  45  cm  depth  in  site  42  in  1997. 
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Figure  72.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  60  cm  depth  in  site  42  in  1997. 
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Figure  73.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  75  cm  depth  in  site  42  in  1997. 
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Figure  74.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  90  cm  depth  in  site  42  in  1997. 
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Figure  75.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  105  cm  depth  in  site  42  in  1997. 
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Figure  76.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  30  cm  depth  in  site  2 in  1998. 
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Figure  77.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  45  cm  depth  in  site  2 in  1998. 
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Figure  78.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  60  cm  depth  in  site  2 in  1998. 
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Figure  79.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  75  cm  depth  in  site  2 in  1998. 
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Figure  80.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  90  cm  depth  in  site  2 in  1998. 
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Figure  81 . Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  105  cm  depth  in  site  2 in  1998. 
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Figure  82.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  120  cm  depth  in  site  2 in  1998. 
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Figure  83.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  135  cm  depth  in  site  2 in  1998. 
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Figure  84.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  150  cm  depth  in  site  2 in  1998. 
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Figure  85.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  30  cm  depth  in  site  42  in  1998. 
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Figure  86.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  45  cm  depth  in  site  42  in  1998. 
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Figure  87.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  60  cm  depth  in  site  42  in  1998. 
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Figure  88.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  75  cm  depth  in  site  42  in  1998. 
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Figure  89.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  90  cm  depth  in  site  42  in  1998. 
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Figure  90.  Observed  (OBS)  and  simulated  soil-water  content  for  each  of  the  estimation 
procedures  (y97  yield-based  97;  y98  yield-based  98;  W97  SWC-based  97;  W98  SWC- 
based  98)  at  105  cm  depth  in  site  42  in  1998. 


LIST  OF  REFERENCES 


Acock,  B.  and  Y.  Pachepsky.  1997.  Holes  in  precision  farming:  mechanistic  crop  models, 
p.  397-404.  In  J.V.  Stafford  (ed.),  Proceedings  of  the  First  European  Conference 
on  Precision  Agriculture.  BIOS  Scientific  Publishers  Ltd.,  Oxford,  UK. 

Batchelor,  W.D.  1998.  Fundamentals  of  neural  networks,  p.  597-628.  In  R.M.  Peart  and 
R.B.  Curry  (ed.),  Agricultural  systems  modeling  and  simulation.  Marcel  Dekker, 
New  York. 

Bell,  J.C.,  C.A.  Butler,  and  J.A.  Thompson.  1995.  Soil-terrain  modeling  for  site-specific 
agricultural  management,  p.  209-227.  In  P.C.  Robert,  R.H.  Rust,  and  W.E.  Larson 
(ed.),  Site-specific  management  for  agricultural  systems.  ASA-CSSA-SSSA, 
Madison,  WI. 

Birrel,  S.  J.  and  Sudduth,  K.  A.  1995.  Corn  population  sensor  for  precision  farming. 
ASAE  paper  No.  95-1334.  ASAE,  St.  Joseph,  ML 

Booltink,  H.W.G.  and  J.  Verhagen.  1997.  Using  decision  support  systems  to  optimize 

barley  management  on  spatial  variable  soil,  p.  219-233.  In  M.J.  Kropff,  P.S.  Teng, 
P.K.  Aggarwal,  J.  Bouma,  B.A.M.  Bouma,  J.W.  Jones,  and  H.H.  van  Laar  (ed.). 
Application  of  systems  approaches  at  the  field  level.  Kluwer  Academic 
Publishers,  Dordrecht,  The  Netherlands. 

Boone,  M.Y.L.  and  J.A.  Landivar.  1999.  Performance  of  ICEMM:  A cotton  simulation 
model  in  a precision  farming  study,  p.  1291-1296.  In  P.C.  Robert,  R.H.  Rust,  and 
W.E.  Larson  (ed.),  Precision  agriculture.  ASA-CSSA-SSSA,  Madison,  WI. 

Bouma,  J.,  J.  Verhagen,  J.  Brouwer,  and  J.M.  Powell.  1997.  Using  systems  approaches 
for  targeting  site-specific  management  on  field  level,  p.  25-36.  In  M.J.  Kropff, 

P.S.  Teng,  P.K.  Aggarwal,  J.  Bouma,  B.A.M.  Bouma,  J.W.  Jones,  and  H.H.  van 
Laar  (ed.),  Application  of  systems  approaches  at  the  field  level.  Kluwer  Academic 
Publishers,  Dordrecht,  The  Netherlands. 

Braga,  R.P.,  J.W.  Jones,  and  B.  Basso.  1999.  Weather  Variability  in  Site-Specific 

Management  Profitability:  A Case  Study,  p.  1853-1863.  In  P.C.  Robert,  R.H. 

Rust,  and  W.E.  Larson  (ed.),  Precision  Agriculture.  ASA-CSSA-SSSA,  Madison, 
WI. 


169 


170 


Calmon,  M.A.,  W.D.  Batchelor,  J.W.  Jones,  J.T.  Ritchie,  K.J.  Boote,  and  L.C.  Hammond. 
1999a.  Simulating  soybean  root  growth  and  soil  water  extraction  using  a 
functional  crop  model.  Trans.  ASAE  42:1867-1877. 

Calmon,  M.A.,  J.W.  Jones,  D.  Shinde,  and  J.E.  Specht.  1999b.  Estimating  parameters  for 
soil  water  balance  models  using  adaptive  simulated  annealing.  Applied 
Engineering  in  Agriculture  15:703-713. 

Cora,  J.E.,  F.J.  Pierce,  B.  Basso,  and  J.T.  Ritchie.  1999.  Simulation  of  within  field 

variability  of  corn  yield  with  Ceres-Maize  model,  p.  1309-1319.  In  P.C.  Robert, 
R.H.  Rust,  and  W.E.  Larson  (ed.),  Precision  agriculture.  ASA-CSSA-SSSA, 
Madison,  WI. 

Crawford,  C.A.G.,  D.G.  Bullock,  F.J.  Pierce,  W.W.  Stroup,  G.W.  Hergert,  and  K.M. 

Eskridge.  1997.  Experimental  design  issues  and  statistical  evaluation  techniques 
for  site-specific  management,  p.  301-335.  In  F.J.  Pierce  and  E.J.  Sadler  (ed.),  The 
State  of  Site-Specific  Management  for  Agriculture.  ASA-CSSA-SSSA,  Madison, 
WI. 

Denholm,  K.A.,  J.LB.  Culley,  J.D.  Aspinall,  and  E.A.  Wilson.  1993.  Soil  landscape 

relations  and  their  influence  on  yield  variability  in  Kent  County,  Ontario,  p.  355- 
356.  In  P.C.  Robert,  R.H.  Rust,  and  W.E.  Larson  (ed.),  Soil  specific  crop 
management.  ASA-CSSA-SSSA,  Madison,  WI. 

Easton,  D.  1996.  Corn  population  and  plant  spacing  variability:  the  next  mapping  layer, 
p.  721-727.  In  P.C.  Robert,  R.H.  Rust,  and  W.E.  Larson  (ed.),  Precision 
agriculture.  ASA-CSSA-SSSA,  Madison,  WI. 

Everett,  M.W.  and  F.J.  Pierce.  1996.  Variability  of  corn  yield  and  soil  profile  nitrates  in 
relation  to  site-specific  N management,  p.  43-53.  In  P.C.  Robert,  R.H.  Rust,  and 
W.E.  Larson  (ed.),  Precision  agriculture.  ASA-CSSA-SSSA,  Madison,  WI. 

Fitzjohn,  C.,  J.L.  Ternan,  and  A.G.  Williams.  1998.  Soil  moisture  variability  in  a semi- 
arid  gully  catchment:  implications  for  runoff  and  erosion  control.  Catena  32:55- 
70. 

Fraisse,  C.W.,  K.A.  Sudduth,  and  N.R.  Kitchen.  1999.  Evaluation  of  crop  models  to 

simulate  site-specific  crop  development  and  yield,  p.  1297-1308.  In  P.C.  Robert, 
R.H.  Rust,  and  W.E.  Larson  (ed.),  Precision  agriculture.  ASA-CSSA-SSSA, 
Madison,  WI. 

Frazen,  D.W.,  A.D.  Halvorson,  J.  Krupinsky,  V.L.  Hofman,  and  L.  Cihacek.  1999. 
Directed  sampling  using  topography  as  a logical  basis,  p.  1559-1568.  In  P.C. 
Robert,  R.H.  Rust,  and  W.E.  Larson  (ed.),  Precision  agriculture.  ASA-CSSA- 
SSSA,  Madison,  WI. 


171 


Gee,  G.E.  and  J.W.  Bauder.  1986.  Particle-size  analysis.  P.383-410.  In  Klute  (ed.) 

Methods  of  soil  analysis.  Part  1 . 2nd  Ed.  Agron.  Monogr.  9.  ASA,  Madison,  WI. 

Goffe,  W.L.,  G.D.  Ferrier,  and  J.R.  Rogers.  1994.  Global  optimization  of  statistical 
function  with  simulated  annealing.  Journal  of  Econometrics  60:65-100. 

Godwin,  D.C.  and  U.  Singh.  1998.  Nitrogen  balance  and  crop  response  to  nitrogen  in 

upland  and  lowland  cropping  systems,  p.  55-77.  In  G.Y.  Tsuji,  G.  Hoogenboom, 
and  P.K.  Thornton  (ed.),  Understanding  Options  for  Agricultural  Production. 
Kluwer  in  cooperation  with  ICASA,  Dordrecht/Boston/London. 

Hall,  G.F.  and  C.G.  Olson.  1991 . Predicting  variability  of  soils  from  landscape  models,  p. 
9-24.  In  M.J.  Mausbach  and  L.P.  Wilding  (ed.),  Spatial  variabilities  of  soils  and 
landforms  - Special  Publication  Number  28.  SSSA,  Madison,  WI. 

Hanna,  A.Y.,  P.W.  Harlan,  and  D.T.  Lewis.  1982.  Soil  available  water  as  influenced  by 
landscape  position  and  aspect.  Agron.  J.  74:999-1004. 

Hashimoto,  Y.  1997.  Applications  of  artificial  neural  networks  and  genetic  algorithms  to 
agricultural  systems.  Comput.Electron.Agric.  18:71-72. 

Hoogenboom,  G.,  P.W.  Wilkens,  W.D.  Batchelor,  W.T.  Bowen,  L.A.  Hunt,  N.B. 

Pickering,  U.  Singh,  D.C.  Godwin,  B.  Baer,  K.J.  Boote,  J.T.  Ritchie,  and  J.W. 
White.  1994.  Crop  models,  p.  95-244.  In  G.Y.  Tsuji,  G.  Uehara,  and  S.  Balas 
(ed.),  DSSAT  version  3 volume  2.  University  of  Hawaii,  Honolulu,  HI. 

Issacks,  E.H.  and  R.M.  Srivastava.  1989.  An  introduction  to  applied  geostatistics.  Oxford 
Univ.  Press.  Oxford. 

Kern,  J.S.  1995.  Evaluation  of  soil  water  retention  models  based  on  basic  soil  physical 
properties.  Soil  Sci.  Soc.  Am.  J.  59:1 134-1 141. 

Khakural,  B.R.,  P.C.  Robert,  and  D.J.  Mulla.  1996.  Relating  corn/soybean  to  variability 
in  soil  and  landscape  characteristics,  p.  1 17-128.  In  P.C.  Robert,  R.H.  Rust,  and 
W.E.  Larson  (ed.),  Precision  agriculture.  ASA-CSSA-SSSA,  Madison,  WI. 

Logsdon,  S.D.,  J.H.  Prueger,  D.  Meek,  T.  Colvin,  and  D.  James.  1999.  Crop  yield 
variability  as  influenced  by  water  in  rain-fed  agriculture,  p.  453-465.  In  P.C. 
Robert,  R.H.  Rust,  and  W.E.  Larson  (ed.),  Precision  agriculture.  ASA-CSSA- 
SSSA,  Madison,  WI. 

MacMillan.  R.A.,  W.W.  Pettapiece,  L.D.  Watson,  and  T.W.  Goddard.  1999.  A landform 
segmentation  model  for  precision  farming,  p.  1335-1346.  In  P.C.  Robert,  R.H. 
Rust,  and  W.E.  Larson  (ed.),  Precision  agriculture.  ASA-CSSA-SSSA,  Madison, 
WI. 


172 


Mallarino,  A.P.,  P.N.  Hinz,  and  E.S.  Oyazabal.  1996.  Multivariate  analysis  as  a tool  for 
interpreting  relationships  between  site  variables  and  crop  yields,  p.  151-158.  In 
P.C.  Robert,  R.H.  Rust,  and  W.E.  Larson  (ed.),  Precision  agriculture.  ASA- 
CSSA-SSSA,  Madison,  WI. 

Moore,  G.A.  and  J.P.  Tyndale-Biscoe.  1999.  Estimation  of  the  importance  of  spatially 
variables  nitrogen  application  and  soil  moisture  holding  capacity  to  wheat 
production.  Prec.  Agric.  1:27-38. 

Moore,  I.D.,  G.J.  Burch,  and  D.H.  Mackenzie.  1988.  Topographic  effects  on  the 

distribution  of  surface  soil  water  and  the  location  of  ephemeral  gullies.  Trans. 
ASAE  31:1098-1 107. 

Moore,  I.D.,  P.E.  Gessler,  G.A.  Nielsen,  and  G.A.  Peterson.  1993.  Terrain  analysis  for 
soil  specific  crop  management,  p.  27-55.  In  P.C.  Robert,  R.H.  Rust,  and  W.E. 
Larson  (ed.),  Soil  specific  crop  management.  ASA-CSSA-SSSA,  Madison,  WI. 

Morton,  Ty.W.,  G.W.  Buchleiter,  and  D.F.  Heermann.  1999.  Quantifying  the  effects  of 

water  availability  on  corn  yield  under  a center-pivot  irrigation  system,  p.  3 1-41 . In 
P.C.  Robert,  R.H.  Rust,  and  W.E.  Larson  (ed.).  Precision  agriculture.  ASA- 
CSSA-SSSA,  Madison,  WI. 

Mulla,  D.J.  1998.  Estimating  spatial  patterns  in  water  content,  matric  suction,  and 
hydraulic  conductivity.  Soil  Sci.  Soc.  Am.  J.  52:1547-1553. 

Mulla,  D.J.  and  J.S.  Schepers.  1997.  Key  processes  and  properties  for  site-specific  soil 
and  crop  management,  p.  1-18.  In  F.J.  Pierce  and  E.J.  Sadler  (ed.),  The  State  of 
Site-Specific  Management  for  Agriculture.  ASA-CSSA-SSSA,  Madison,  WI. 

Nolin,  M.C.,  S.P.  Guertin,  and  C.  Wang.  1996.  Within-field  spatial  variability  of  soil 
nutrients  and  corn  yield  in  a Montreal  lowlands  clay  soil,  p.  255-270.  In  P.C. 
Robert,  R.H.  Rust,  and  W.E.  Larson  (ed.).  Precision  agriculture.  ASA-CSSA- 
SSSA,  Madison,  WI. 

Pannatier,  Y.  1996.  Variowin:  software  for  spatial  data  analysis  in  2D.  Springer-Verlag, 
New  York,  Inc.,  Edwards  Brothers,  Inc.,  Ann  Arbor,  MI. 

Paz,  J.O.,  W.D.  Batchelor,  T.S.  Colvin,  S.D.  Logsdon,  T.C.  Kaspar,  and  D.L.  Karlen. 
1997.  Analysis  of  Water  Stress  Effects  Causing  Spatial  Yield  Variability  in 
Soybeans.  Trans.  ASAE  41 : 1527-1 534. 

Paz,  J.O.,  W.D.  Batchelor,  B.A.  Babcock,  T.  Colvin,  S.D.  Logsdon,  T.C.  Kaspar,  and 

D.L.  Karlen.  1999.  Model-based  technique  to  determine  variable  rate  nitrogen  for 
corn.  Agric.  Syst.  61 :69-75. 


173 


Pierce,  F.J.  and  E.J.  Sadler.  1997.  Preface,  F.J.  Pierce  and  E.J.  Sadler  (ed.),  The  State  of 
Site-Specific  Management  for  Agriculture.  ASA-CSSA-SSSA,  Madison,  WI. 

Plattner,  C.E.  and  J.W.  Hummel.  1996.  Corn  plant  population  sensor  for  precision 
agriculture,  p.  785-794.  In  P.C.  Robert,  R.H.  Rust,  and  W.E.  Larson  (ed.). 
Precision  agriculture.  ASA-CSSA-SSSA,  Madison,  WI. 

Pommel,  B.  and  R.  Bonhomme.  1998.  Variations  in  the  vegetative  and  reproductive 

systems  in  individual  plants  of  an  heterogeneous  maize  crop.  European  Journal  of 
Agronomy  8:39-49. 

Press,  W.H.,  S.A.  Teukolsky,  W.T.  Vetterling,  and  B.P.  Flannery.  1992.  Numerical 
recipes  in  FORTRAN  77  and  FORTRAN  90:  The  art  of  scientific  and  parallel 
computing.  Cambridge  University  Press,  Cambridge,  UK. 

Priestley,  C.H.B.  and  R.J.  Taylor.  1972.  On  the  assessment  of  surface  heat  flux  and 

evaporation  using  large-scale  parameters.  Monthly  Weather  Review  100:81-92. 

Ritchie,  J.T.  1972.  Model  for  prediction  evaporation  from  a row  crop  with  incomplete 
cover.  Water  Resour.  Res.  8:1024-1213. 

Ritchie,  J.T.  1998.  Soil  water  balance  and  plant  water  stress,  p.  41-54.  In  G.Y.  Tsuji,  G. 
Hoogenboom,  and  P.K.  Thornton  (ed.),  Understanding  Options  for  Agricultural 
Production.  Kluwer  in  cooperation  with  ICASA,  Dordrecht/Boston/London. 

Ritchie,  J.T.,  D.C.  Godwin,  and  U.  Singh.  1989.  Soil  and  weather  inputs  for  the  IBSNAT 
crop  models.,  p.  31-45.  In  Proceedings  of  the  IBSNAT  Symposium:  Decision 
support  systems  for  agrotechnology  transfer.  University  of  Hawaii,  Honolulu,  HI. 

Ritchie,  J.T.,  U.  Singh,  D.C.  Godwin,  and  W.T.  Bowen.  1998.  Cereal  growth, 

development  and  yield,  p.  79-98.  In  G.Y.  Tsuji,  G.  Hoogenboom,  and  P.K. 
Thornton  (ed.),  Understanding  Options  for  Agricultural  Production. 
Kluwer/ICASA,  Dordrecht/Boston/London. 

Robert,  P.C.,  R.H.  Rust,  and  W.E.  Larson.  1995.  Preface,  P.C.  Robert,  R.H.  Rust,  and 
W.E.  Larson  (ed.),  Site-specific  management  for  agricultural  systems.  ASA- 
CSSA-SSSA,  Madison,  WI. 

Sadler,  E.J.,  W.J.  Busscher,  and  D.L.  Karlen.  1995.  Site-specific  yield  histories  on  a SE 
coastal  plain  field,  p.  153-166.  In  P.C.  Robert,  R.H.  Rust,  and  W.E.  Larson  (ed.), 
Site-specific  management  for  agricultural  systems.  ASA-CSSA-SSSA,  Madison, 
WI. 

Sadler,  E.J.,  D.E.  Evans,  W.J.  Busscher,  and  D.L.  Karlen.  1993.  Yield  variation  across 
coastal  plain  soil  mapping  units,  p.  373-374.  In  P.C.  Robert,  R.H.  Rust,  and  W.E. 
Larson  (ed.),  Soil  specific  crop  management.  ASA-CSSA-SSSA,  Madison,  WI. 


174 


Saxton,  K.E.,  W.J.  Rawls,  J.S.  Romberger,  and  R.I.  Papendick.  1986.  Estimating 

generalized  soil-water  characteristics  from  texture.  Soil  Sci.Soc.Am.J.  50:1031- 
1036. 

Schueller,  J.K.  1996.  Impediments  to  spatially-variable  field  operations.  Comput. 
Electron.  Agric.  14:249-253. 

Schultz,  A.  and  R.  Wieland.  1997.  The  use  of  neural  networks  in  agricultural  modelling. 
Comput.  and  Electron.  Agric.  18:73-90. 

Sudduth,  K.A.,  N.R.  Kitchen,  D.F.  Hughes,  and  S.T.  Drummond.  1995.  Electromagnetic 
induction  sensing  as  an  indicator  of  productivity  on  claypan  soils,  p.  671-681 . In 
P.C.  Robert,  R.H.  Rust,  and  W.E.  Larson  (ed.),  Site-specific  management  for 
agricultural  systems.  ASA-CSSA-SSSA,  Madison,  WI. 

Sudduth,  K.A.,  S.T.  Drummond,  S.J.  Birrel,  and  N.R.  Kitchen.  1996.  Analysis  of  spatial 
factors  influencing  crop  yield,  p.  129-139.  In  P.C.  Robert,  R.H.  Rust,  and  W.E. 
Larson  (ed.),  Precision  agriculture.  ASA-CSSA-SSSA,  Madison,  WI 

Sudduth,  K.A.,  C.W.  Fraisse,  S.T.  Drummond,  and  N.R.  Kitchen.  1998.  Integrating 
spatial  data  collection,  modeling  and  analysis  for  precision  agriculture,  p.  1 66- 
173.  In  Proc.of  the  First  International  Conference  on  Geospatial  Information  in 
Agriculture  and  Forestry.  ERIM  International,  Inc.,  Ann  Arbor,  MI. 

Threlkeld,  G.W.  and  J.E.  Feenstra.  1974.  Soil  survey  of  Shiawassee  County,  Michigan. 
USD  A,  SCS  and  Michigan  Agricultural  Experiment  Station,  Washington,  DC. 

Timlin,  D.J.,  Y.  Pachepsky,  V.A.  Snyder,  and  R.B.  Bryant.  1999.  Soil  surface  and  the 
spatial  and  temporal  variability  of  corn  grain  on  a hillslope,  p.  203-214.  In  P.C. 
Robert,  R.H.  Rust,  and  W.E.  Larson  (ed.),  Precision  agriculture.  ASA-CSSA- 
SSSA,  Madison,  WI 

Torner,  M.D.  and  J.L.  Anderson.  1995a.  Variation  of  soil  water  storage  across  a sand 
plain  hillslope.  Soil  Sci.  Soc.  Am.  J.  59:1091-1 100 

Torner,  M.D.,  J.L.  Anderson,  and  J.A.  Lamb.  1995b.  Landscape  analysis  of  soil  and  crop 
data  using  regression,  p.  273-284.  In  P.C.  Robert,  R.H.  Rust,  and  W.E.  Larson 
(ed.),  Site-specific  management  for  agricultural  systems.  ASA-CSSA-SSSA, 
Madison,  WI. 

Van  Uffelen  , C.G.R.,  J.  Verhagen,  and  J.  Bouma.  1997.  Comparison  of  simulated  crop 
yield  patterns  for  site-specific  management.  Agric.  Syst.  54:207-222. 


175 


Verhagen,  A.,  H.W.G.  Booltink,  and  J.  Bouma.  1995.  Site-specific  management: 

Balancing  production  and  environmental  requirements  at  farm  scale.  Agric.  Syst. 
49:369-384. 

Williams,  J.R.  1991.  Runoff  and  water  erosion,  p.  439-455.  In  R.J.  Hanks  and  J.T. 
Ritchie  (ed.),  Modeling  plant  and  soil  systems.  Agronomy  Monograph  #31. 
American  Society  of  Agronomy,  Madison,  WI. 


BIOGRAPHICAL  SKETCH 


Ricardo  Nuno  da  Fonseca  Garcia  Pereira  Braga  was  born  in  Lisbon,  Portugal,  on 
September  1970.  He  started  studying  agriculture  in  high  school  taking  a 3-year 
specialization  in  'agro-pecuaria'  (animal  and  crop  production).  In  September  1988  he 
enrolled  at  the  Technical  University  of  Lisbon,  Agronomy  Superior  Institute  to  pursue  the 
5-year  degree  of  "Engenharia  Agronomica"  (Agronomy  Engineering).  In  1993  he 
received  this  degree  and  won  an  award  from  "Ordem  dos  Engenheiros"  (National 
Engineering  Association)  for  being  the  graduating  student  with  the  highest  score.  In  July 
1996  he  received  the  degree  of  Master  of  Science  in  Tropical  Agricultural  Production 
from  the  Technical  University  of  Lisbon.  His  Master  of  Science  thesis  entitled 
'Evaluation  of  CERES-Rice  Model,  Sensitivity  Analysis,  Validation  and  Weather  Data 
Quality'  was  supervised  by  Dr.  Pedro  A.  Pinto.  A month  later,  he  enrolled  in  the 
Agricultural  and  Biological  Engineering  Department,  University  of  Florida,  to  pursue  a 
Ph.D.  degree  under  the  supervision  of  Dr.  J.W.  Jones.  He  was  awarded  two  scholarships 
for  his  Ph.D.  education:  The  Science  and  Technology  Foundation  Praxis  XXI  Scholarship 
(Co-Funded  by  the  European  Union  ) and  the  Fulbright  Scholarship.  In  August  2000, 
Ricardo  received  his  Ph.D.  degree  after  defending  the  dissertation  entitled  'Predicting  the 
Spatial  Pattern  of  Com  Yield  Under  Water  Limiting  Conditions'. 


176 


I certify  that  I have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a dissertation  for  the  degree  of  Doctor  of  Philosc 


J/me^  W.  Jones,  Chair 

unguished  Professor  of  Agricultural 
and  Biological  Engineering 

I certify  that  I have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Joe^T I.  Ritchie 
Prqjessor  of  Crop  and  Soil  Sciences 
Michigan  State  University 

I certify  that  I have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Kenneth  J.  Boote^-* 
Professor  of  Agronomy 


I certify  that  I have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a dissertation  for  the  degree  of  Doctor  of  Philosophy. 

l/wtjAu 


Wendy  D.Ajranam 
Professor  (of  Agricultural  arf 
Engineering 


lological 


I certify  that  I have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Pedro  A.  Pinto 

Professor  of  Plant  and  Animal  Production 
Technical  University  of  Lisbon 


This  dissertation  was  submitted  to  the  Graduate  Faculty  of  the  College  of 
Agriculture  and  Life  Sciences  and  to  the  Graduate  School  and  was  accepted  as  partial 
fulfillment  of  the  requirements  for  the  degree  of  Doctor  of  Philosophy. 


August  2000 


£ . 


Dean,  College  of  Agricultur 


Dean,  Graduate  School 


